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Die Aufgaben einer Versuchsanstalt für Luftfahrt 


in Deutschland vor 25 Jahren. 
Von Wilhelm Hoff in Berlin. 


or 25 Jahren war die Zeit der ersten Luftfahrtpioniere vorüber. Ohne sich um die 

Meinung der Umwelt zu kümmern, hatten sie sieh für den uralten Menschheitstraum, es 
den Vögeln gleichzutun und zu fliegen, mit Gut und Blut eingesetzt. Viel war schon er- 
reicht. Luftschiffe verschiedener Bauart kreisten über deutschem Land. Ihre Konstrukteure 
und Besatzungen sammelten reichlich Erfahrungen, die grundlegend für den großen Auf- 
schwung und die Vormachtstellung des deutschen Luftschiffbaues werden sollten. Auf Flug- 
wochen und groß angelegten Rundflügen wurden die ersten brauchbaren Flugzeuge gezeigt. 
Die Führer hatten selbst ihre Beherrschung für Flüge bei günstigem Wetter erlernt, und 
waren dabei, die Schwierigkeiten des Schleehtwetterfluges zu überwinden und die ersten 
Figuren des Kunstfluges zu üben. Zusammen mit Sportfliegern hatten sich Hugbegeisterte 
Oftiziere des Heeres und der Marine, oft unter erheblichen Widerständen ihrer Dienststellen. 
daran beteiligt. Die ersten militärischen Schulungskurse waren im Gange. Eine Luftfahrt- 
industrie war im Entstehen. Sie bemühte sieh, die Kenntnisse verwandter technischer Ge- 
biete aufzugreifen und neue, eigene Grundlagen für die Entwieklung des Flugzeugbaues zu 
schaffen. Eine große Zahl geschulter Ingenieure hatte ihre Berufsaufgabe in der Luftfahrt 
eefunden. 


Nachdem Anfang November 1911 sich Vertreter der Flugwissenschaft in Göttingen 
getroffen und kennengelernt hatten, war am 3. April 1912 die Wissenschaftliche Gesellschaft 
für Flugtechnik (später „für Luftfahrt“, kurz „WGL* genannt) gegründet worden. In Göttingen 
arbeitete seit 19085 die Modell-Versuchsanstalt. An den Technischen Hochschulen in Aachen 
und Stutteart hatten neue Institute ihre Arbeiten begonnen. "An vielen anderen Technischen 
Hochschulen und Universitäten wurden die neuen Probleme dureh Vorlesungen an interessierte 
Studenten herangebracht. Die ersten Verbindungen mit dem Ausland wurden geknüpft. 


In diese Zeit des allgemeinen Aufstieges fiel die am 20. April 1912 erfolgte Gründung 
der Deutschen Versuchsanstalt für Luftfahrt, später allgemein „DVL* benannt. Die Vorver- 
handlungen hatten einige Monate gedauert. Anregungen aus dem Reichstag und dem Verein 
Deutscher Motorflugzeugbau-Industrieller waren erfolgt. Eine unmittelbare Weisung des 
deutschen Kaisers hatte den entscheidenden Anstoß gegeben. Das Reichsamt des Innern, 
damals vertreten durch Ministerialdirektor Dr. Lewald und Geh. Oberregierungsrat Albert, 
übernahm die Führung, tatkräftig unterstützt von allen beteiligten Zivil- und Militärbehörden, 
der Industrie, sowie von Vereinen und Gönnern der Luftfahrt. Die DVL wurde nicht als 
„Reichsanstalt“ geschaffen, sondern ein eingetragener Verein übernahm die Pflicht, „das 
deutsche Flugwesen und die deutsche Luftschiffahrt durch Errichtung, Ausbau und Unter- 


> 
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haltung einer Versuchsanstalt zu gemeinsamem Nutzen zu fördern.“ Dem deutschen Reichs- 
kanzler waren wichtige Entscheidungen in der Zielsetzung und Verwaltung dieses eingetragenen 
Vereins vorbehalten. 


Trotz vieler Stürme um ihren Bestand und ihre Gestaltung, ist der Grundgedanke, aus 
der DVI, keine Reichsanstalt mit der damit verbundenen behördlichen Schwerfälligkeit zu 
machen, bis heute erhalten geblieben. Die 25 Lebensjahre der DVL brachten manchen Bei- 
trag zu der oft erörterten Frage, ob eine wissenschaftliche Tätigkeit geeigneter im Rahmen 
einer Reichsanstalt, wie er früher in der Physikalisch-Technischen Reichsanstalt gefunden 
wurde, oder in der Art der späteren Kaiser Wilhelm-Institute durchzuführen ist. In diesem 
Aufsatz wird nieht auf dieses sehr interessante Gebiet, jedoch darauf eingegangen werden, 
wie die ersten Aufgaben der DVL zu ihrer Gründungszeit gesehen worden sind, und wie wir 
heute zurückschauend zu diesen Fragen stehen. 


In ausführlichen Denkschriften hatten im Sommer 1911 einerseits Geh. Regierungsrat 
Prof. Dr. H. Hergesell und Prof. Dr. L. Prandtl, sowie andererseits Dr.-Ing. F. Bende- 
mann, der spätere erste Direktor der DVL auf Ersuchen des Staatssekretärs des Reichs- 
amtes des Innern die Aufgaben einer Versuchsanstalt umrissen. 


Wie schnell schon damals die Entwicklung ging, ist daran zu erkennen, daß nur ein 
Teil der skizzierten Aufgaben später tatsächlich übernommen worden ist, und andere wich- 
tigere und umfassendere an deren Stelle getreten sind. 


Um es vorwegzunehmen, das Gebiet „des Windes und der Atmosphäre“ fand in der 
DVI, keine Aufnahme. Nur eine allerdings nieht abgeschlossene Forschungsarbeit dieses 
Arbeitsfeldes sei hier erwähnt. Zur Erforschung von Windströmungen in der Nähe von 
Bodenhindernissen sollten Rauchringe durch eine besonders hierfür gebaute Maschine erzeugt 
und hochgeblasen werden. Meteorologische und aerologische Arbeiten sind in Deutschland 
von anderen Stellen in so weitem Umfang aufgegriffen worden, daß es kein Fehler war, wenn 
sich die DVL stets von diesem Wissensgebiet ferngehalten hat. 


Dafür traten andere Fachrichtungen in den Vordergrund und blieben es auch im wesent- 
lichen in den ersten Jahren der DVL-Lebensgeschichte. 


Damit der erste Wettbewerb um den Kaiser-Preis für den besten deutschen Flugzeugmotor 
durehgeführt und sein Ergebnis im Januar 1913 vorgelegt werden konnte, mußten zunächst 
alle verfügbaren Kräfte auf dies im Brennpunkt des Luftfahrtinteresses stehende Ereignis 
angesetzt werden. Viel Zeit war nicht gegeben. Es galt zuerst einwandfreie Prüfstände zu 
schaffen, auf denen Flugzeugmotoren auf Höchstleistung unter Einhaltung ihrer späteren Be- 
triebsbedingungen, d.h. mit laufender Luftschraube und im Prüfrahmen wie im Flugzeug- 
rumpf gebettet, durchgemessen und auf Dauerbetrieb geprüft werden konnten. 


Die meßtechnische Hauptschwierigkeit lag in einwandfreier Ermittelung des auf die 
laufende Schraube wirkenden Drehmoments, das von einem Pendelralmen aus gemessen 
werden sollte, die Rückwirkung des Schraubendralls und der Auspuffgase war zu beachten. 


Die DVI, hatte die Genugtuung, daß das von ihr geschaffene Hilfsmittel sofort einschlug 
und zu keinen wesentlichen Beanstandungen führte. Alle zur Motorenbeurteilung notwendigen 
(Größen (z. B. Betriebsstoffverbrauch) konnten befriedigend erfaßt werden. Manche wertvolle 
Beobachtung wurde gesammelt, später verarbeitet und an die beteiligten Kreise weitergegeben. 


Im Wettbewerb hatten 26 verschiedene Arten vier- und sechszylindriger wassergekühlter 
stehender Motoren, sowie sechs- und siebenzvlindriger luftgekühlter Umlaufmotoren gestanden. 
Die Nutzleistungen bewegten sich zwischen etwa 45 PS bei den kleineren und etwas über 
I PS bei den größeren Motoren. Zum Vergleich wurden die Einheitsgewichte (einschließ- 
lich Kühler, Behälter und Betriebsstoffe 7 h) herangezogen, die sich auf 3,8 kg/PS bei dem 
Sieger, und bei dem schwersten Umlaufmotor auf 5,85 kg/PS beliefen. Der siegreiche vier- 
zylindrige Motor von Benz & Co. hatte insgesamt 13,5 h, zuletzt unter erschwerten Be- 
dingungen zu laufen. Er hatte eine Einheitsleistung von 10,75 PS/1 (bei einer Drehzahl von 
1255 U/min) erbracht. Besonders zu erwähnen war außerdem ein vierzylindriger Motor der 
Daimler-Motoren-Ges. mit hängenden Zylindern und Stirnradgetriebe. 


Der Motor, der sich als „bester* aus für damalige Zeit sehr scharfer Prüfung hervor 
getan hatte, fand nieht den Beifall der „zünftigen* Fachwelt. Ein anderer Motor wurde als 
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‚moralischer“ Sieger angesehen. Damit mußte die DVL sich zum erstenmal, wie später auch bei 
anderen Gelegenheiten, heftiger Angriffe wegen Weltfremdheit und mangelnden Verständ- 
nisses für die praktischen Bedürfnisse erwehren. Die nachfolgende Motorenentwicklung hat 
dem Wettbewerbsergebnis recht gegeben. Aber nur wenige erinnerten sich dann an die der 
DVL zu Unrecht gemachten Vorwürfe. 


Durch diesen Motorenwettbewerb war die DVL mit Erzeugern und Verbrauchern von 
Flugzeugmotoren in enge Verbindung gekommen. Zur Weiterentwicklung neuer Baumuster 
und ihres Zubehörs wurde die DVL gern herangezogen. Damit war sie in dies wichtige 
Teilgebiet der Luftfahrttechnik eingeschaltet. 


Die von der DVL geleisteten Vorarbeiten wurden riehtunggebend für Untersuchungs- 
verfahren an Luftfahrtmotoren. Viele deutsche Prüffelder wurden mit Motorprüfständen nach 
DVI Vorbild und zugehörigem Meßgerät ausgerüstet. 


In der Natur dieser Arbeit lag es, daß, nachdem der Bericht über das Ergebnis des 
Kaiserpreis-Wettbewerbs erstattet worden war, wenig in Fachzeitschriften bekanntgegeben 
werden konnte. Auch wurde bald der zweite Kaiserpreis-Wettbewerb vorbereitet, der ım 
Winter 1914/15 auf weiterentwickelten Motorenprüffeldern und mit verschärften Bedingungen 
ausgetragen werden sollte. 


Heute wird man sich fragen, warum nieht schon früher die für die Motorenentwicklung 
so wichtigen Einzylinderversuche aufgegriffen worden sind. Dies lag wohl daran, daß die 
damaligen Motorenmuster im ganzen eine Leistung aufbrachten, wie sie heute von einem 
Zylinder abgegeben wird. Vor 25 Jahren stand der schnellaufende Verbrennungsmotor noch 
in seiner Kinderzeit. Neue Erfahrungen im Gesamttriebwerk waren zu sammeln. Die da- 
maligen Auschauungen über Festigkeit waren nicht anwendbar, da sie überhohe Baugewichte 
zur Folge gehabt hätten. Erst auf viele Versuchsergebnisse gestützt setzten hier später 
theoretische Betrachtungen ein. 


Als weiteres Fachgebiet wurden in unmittelbarem Anschluß an die schon 1908 begonnenen 
Luftschraubenuntersuchungen der Geschäftsstelle für Flugtechnik des Sonderausschusses der 
Jubiläumsstiftung der deutschen Industrie in Lindenberg die Prüfung und die Eichung von 
Luftschrauben aufgenommen. Die DVL hatte hierzu Veranlassung, da die Wettbewerbs- 
motoren mit Luftschrauben gebremst werden sollten. Die Lindenberger Arbeiten hatten 
wichtige grundsätzliche Ergebnisse über Hubschrauben erbracht, die es Deutschland ermög- 
lieht haben, in der Hubschraubenfrage zunächst eine abwartende Stellung einzunehmen. Es 
fehlten aber die wichtigen Fahrtversuche mit Luftschrauben. Vornehmlieh diesen sollten 
Schleppversuche mit Schienenfahrzeugen dienen. Die Meßmethoden wurden ausgearbeitet, 
die Versuchswagen vorbereitet. Während der ersten Kriegszeit ruhten diese Arbeiten und 
wurden erst 1916 auf der Streeke der Militäreisenbahn bei Dümde neu aufgenommen. Die 
ersten Entwürfe für Meßnaben zur Bestimmung des Schubs und Drehmoments an Luft- 
schrauben wurden begonnen. 


Mit den Eichungen der Luftschraube wurden regelmäßig Schleuderprüfungen verbunden, 
soweit der Leistungsbereich der ersten Prüfstände dies zuließ. Manche ungeeignete Schraube 
konnte wegen Flatterns oder wegen Anrisse am Kantenbeschlag oder anderen Stellen aus- 
geschieden werden. Somit war die DVL anfangs nur an die Standversuche mit Luftschrauben 
gebunden, denen damals und während des Krieges eine sehr hohe Bedeutung zukam. Luft- 
schrauben mit im Stand einstellbaren Blättern waren noch unbekannt. Zu erwähnen ist, daß 
während des Weltkrieges viele zehntausend Luftschrauben auf DVL-Prüfständen geeicht und 
auf ihr Verhalten geprüft worden sind. 


Die geplanten und immer für vordringlich gehaltenen Fahrtversuche mit Luftschrauben 
kamen nicht zur Durchführung. Dies ist heute gesehen sehr zu bedauern, denn dadurch 
haben die wertvollen Lindenberger Arbeiten nur eine einseitige Fortsetzung gefunden. 


Bei ihrer Gründung wurde der DVL eine eigenartige Vorrichtung zum Messen des 
Winddrucks auf Bauwerke, Dächer, Schornsteine usw. in Obhut gegeben. Dieser nach Plänen 
des Marineingenieurs Gießen gebaute Apparat war das Ergebnis eines Wettbewerbs gewesen, 
der unter Führung des Preußischen Ministeriums der öffentlichen Arbeiten ausgeschrieben 
worden war. Dieses an sich sorgfältig durchdachte Meßgerät kam nicht mehr in Gebrauch, 
da andere einfachere Untersuchungsverfahren, hauptsächlich in Windkanälen, die gestellte 
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Aufgabe schneller und vielseitiger zu lösen vermochten, und da die Versuche auf der soeben 
erwähnten Prüfbahn ergänzend zu Windkanalversuchen eingesetzt werden sollten. 


Während noch bei der Gründung der DVL es als selbstverständlich angesehen wurde, 
daß ein Windkanal oder ein größerer Rundlauf für aerodynamische Versuche verfügbar sein 
sollte, wurden Arbeiten auf diesem engeren Gebiet der Aerodynamik nicht aufgenommen. 
In der ersten Zeit ihres Bestehens wollte sich die DVL auf die in Modellversuchsanstalten. 
besonders in Göttingen gewonnenen Ergebnisse stützen. Sie setzte sich ihrerseits zum Ziel. 
den vielen am ausgeführten Flugzeug auftretenden Fragen näherzukommen. Zu erwähnen 
ist aber, daß damals ein kleiner, aber sehr brauehbarer, mit primitivsten Hilfsmitteln erstellter 
Windkanal Eiffelscher Bauart entstanden ist, der zu Instrumenten-Untersuchungen benutzt 
wurde, und der zwölf Jahre im Betrieb geblieben ist, bis ein Brand ihn vernichtet hat. 


Wenn man heute rückschauend die DVL-Versuchsarbeiten beurteilen will, die dem 
lugzeug galten, so wird man das Schrifttum der damaligen Zeit heranziehen. Wir finden 
es im wesentlichen in den ersten Jahrgängen der Zeitschrift für Flugtechnik und Motorluft- 
schiffahrt. Die ersten Versuehsberichte über gewölbte Tragflügel lagen in dem prachtvoll 
ausgestatteten französischen Werk von G. Eıiffel vor. Man war sich noch nieht im klaren. 
ob man einfache dünne Zylinderschalen wählen, wie man die Eintrittskante gestalten, wohin 
man die größte Durchwölbung legen sollte; man wußte nicht die Bedeutung des Seiten- 
verhältnisses richtig einzuschätzen, die beobachtete höhere Tragwirkung auf den Oberflügel 
eines Doppeldeekers war noch ungeklärt. Die Gesetze der „Druckpunktswanderung* wurden 
erörtert. Grundlegende Arbeiten über die Längsstabilität, die Seitensteuerung und Seiten- 
stabilität waren erschienen, die Eigenschaften der Luftschraube in Fahrt wurden theoretisch 
beleuchtet und dargestellt. 


Technische Kenntnisse und brauchbare Beobachtungen über das Verhalten am ausgeführten 
"luezeug fehlten. An bemannte Drachenversuche, die in Lindenberg zur Feststellung der 
l.uftkräfte gemacht worden waren, und an eine Arbeit, die in der Technischen Hochschule 
Aachen entstanden war, konnte angeknüpft werden. 


Als erste Aufgabe wurde von der DVL in Angriff genommen, brauchbare Unterlagen 
über die Belastung des Flugzeuges im Flugbetrieb zu erhalten. Hierfür waren neue Me&- 
verfahren zu entwickeln, die unter Verwendung einer von F. Bendemann angegebenen, 
gesteuerten Meßdose gefunden wurden. Durch viele Versuchstlüge auf Eindeckern und Doppel- 
deekern wurden mit diesem Hilfsmittel die in wichtigsten Trag- und Steuerorganen währen 
(des Betriebes auftretenden Kräfte gemessen und damit ein Bild über Belastung des Flugzeugs 
im ungestörten Flug, beim Abfangen und im Böenflug gewonnen. Die Ergebnisse dieser 
Messungen brachten geeignete Unterlagen für die in Ansatz zu bringenden Lasten bei den 
Festigkeitsrechnungen der späteren Flugzeuge. 


Bei den ersten Flugzeugmustern war man bestrebt, Tragflügel und sonstige wichtige 
Flugzeugbauteile dureh eine Mindestzahl von gedrückten oder gebogenen Stäben, die als 
entscheidend für die Gewichtshöhe angesehen wurden, und mehr durch viele Zugorgane 
zegenseitig miteinander zu verbinden. Über die anzuwendenden Berecehnungsverfahren herrschte 
eroße Unsicherheit. 

Als die preußische Fliegertruppe mehrere Besatzungen durch Flügelbruch während des 
“luges verloren hatte, wurden die ersten Flugzeugmuster auf ihre Festigkeit, zunächst in 
Döberitz, alsdann aber systematisch dureh die DVL untersucht. Mancher Baufehler konnte 
dabei aufgedeekt werden. Regelmäßig wurde dabei versucht, auch einen Beitrag für die 
Festigkeitsreehnung der Flugzeuge zu liefern. Das erste Ergebnis einer ausführlichen Flügel- 
rechnung unter Verwendung der Methoden der Baustatik konnte kurz vor dem Kriege mit- 
eeteilt werden. 


Bei den großen Flugwettbewerben arbeiteten DVL-Beauftragte im technischen Stabe der 
Veranstalter mit. So wurde ebenso wie es auf dem Fachgebiet der Motoren schon einige 
Monate früher der Fall gewesen war, eine enge Fühlung zwischen Herstellern und Verbrauchern 
von Flugzeugen und der DVL herbeigeführt. 


Die damals für das Flugzeug begonnenen Prüf- und Forschungsarbeiten blieben lange 
Jahre ein Hauptbestandteil der DVL-Tätigkeit. Die Notwendigkeit ihrer Durchführung wurde 
auch in den schweren Nachkriegszeiten immer anerkannt und gab damit der DVL einen 
kückhalt für ıhr Bestehen. 
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Eine Fülle von allgemeinen Aufgaben wurde an die DVL herangetragen. Zu Beginn 
des praktischen Fliegens wurde nur nach „Gefühl“ und nach „Sicht“ geflogen. Ganz allmählich 
erst wurden Meßgeräte zur Unterstützung des Flugzeugführers übernommen. Am wichtigsten 
waren die Motorüberwachungsgeräte. Der Umdrehungszähler wurde auch zur Beobachtung 
der Flugzeuggeschwindigkeit mit verwendet. Später übernahmen Windräder und Staudruck- 
meßgeräte diese Aufgabe. An der Entwicklung und Eiehung von zuverlässigen Motordreh- 
zählern. von Windmessern, Staudruckmeßgeräten und Höhenmessern war die DVL stark 
beteiligt. Dazu kamen viele Sonderfragen, die von der DVI als maßgebende Stelle beantwortet 
werden mußten. Größeren Umfang nahmen die ersten Baustoffuntersuchungen an. 


Der Mitarbeiterkreis der DVL bestand im ersten Jahr ausschließlich aus Ingenieuren, 
denen die Arbeit am Reißbrett und am Prüfstand und der Verkehr mit Herstellern und Ver- 
hrauchern vertraut war, aber wenig Neigung zeigten, gewonnene Erfahrungen zusammenzufassen 
und in Forsehungsberichten niederzulegen. Dabei wurden anfangs in der DVL unter Prüfungs- 
arbeiten solehe Arbeiten verstanden, die gegen Bezahlung in auswärtigem Auftrag ausgeführt 
wurden, und unter Forschungsarbeiten Arbeiten aus eigenem Antrieb und auf DVL-Kosten. 
In der ersten Zeit waren auch diese „Prüfungsarbeiten* stets als Forschungsarbeiten ım 
heutigen Sinne angelegt, und meist die weit anregenderen und interessanteren „Forschungs- 
arbeiten“ im alten Sinne füllten empfindliche Lücken in der Erkenntnis aus. Hier wurde 
versucht, den Wirkungsgrad der DVL-Tätigkeit im Innern zu steigern und für ihr Ansehen 
in der teehnischen Welt nach außen zu sorgen. Vor dem Kriege wurde das erste DVL-Jahrbuch 
herausgebracht. dem allerdings erst 1926 weitere folgen sollten. 


Die Geschieke der jungen DVL wurden von dem ersten Vereinspräsidenten General- 
leutnant z. D. Rieß von Scheuernschloß in wohlwollender Weise behütet. Der Sachverwalter 
im Reiehsamt des Innern wurde Geh. Oberregierungsrat Albert. In den weiteren Vorstand, 
den Verwaltungs- und Technischen Ausschuß wurden namhafte Männer gewählt oder vom 
Reiehskanzler berufen. Als erster Leiter wurde Prof. Dr.-Ing. F. Bendemann bestellt, der 
sieh mit einem Mitarbeiterkreis umgab, der in ihm den vorbildlichen Leiter und fähigen 
Ingenieur verehrte. 


Bei Ausbruch des Weltkrieges wurden die Anlagen‘ der DVL vorbehaltlos für den 
technischen Dienst der Fliegertruppe zur Verfügung gestellt. Ihre ersten Mitarbeiter, die sich 
bei ernster Berufsarbeit und froher Erholung (Kegelspiel) zusammengefunden hatten, wurden 
auseinander gerissen, um ihre Wehrpflicht zu erfüllen. Erst später wurden viele zurück- 
gerufen, um neue dringende technische Aufgaben zu übernehmen. 


Die DVL verlor Dr. Dipl.-Ing. Fuhrmann, Dipl.-Ing. Linek und Dipl.-Ing. Tütschulte 
durch Tod an der Front, Generalleutnant Rießk von Scheuernschloß dureh schwere im Feld 
erworbene Krankheit und Dr.-Ing. Heller dureh tödlichen Flugunfall, außerdem eine Reihe 
tüchtiger Angehöriger ihres Betriebs. An das Wirken dieser vorbildliehen und ausgezeichneten 
Männer sei an «dieser Stelle dankbar und ehrend gedacht. 


In diesem Aufsatz wurde versucht, die ersten Aufgaben einer Deutschen Versuchsanstalt 
für Luftfahrt zu umreißen. Der jüngere Leser wird sich die Frage stellen, warum diese oder 
jene Aufgabe, die heute Lebensbedeutung besitzt, nicht früher erkannt und bearbeitet worden 
ist. Diesem sei gesagt, daß trotz des schönen Aufschwunges vor 25 Jahren die Luftfahrt 
noch in den Anfängen gestanden hat, und daß die DVL im ersten Jahr nach ihrer Gründung 
dem Umfang nach etwa einem heutigen DVL-Institut entsprochen hat. 


Die DVL ist heute groß geworden und mit ausgezeichneten neuartigen Einriehtungen 
ausgerüstet. Sie bildet ein notwendiges Glied für die technische Arbeit in der deutschen 
Luftfahrt. Große und wichtige Aufgaben sind ihr vom Herrn Reiehsminister der Luftfahrt 
gestellt. Für die kommenden Jahre mögen ihrer Leitung und ihren Mitarbeitern schöne 
Erfolge beschieden sein, dies ist der aufriehtige Wunsch aller Luftfahrtingenieure in den 
Tagen des Rückblieks auf ihr 25 jähriges Bestehen. 2 
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Die Hubschraube in Bodennähe. 


Von A. Betz VDI in Göttingen. 


1. Einleitung. Wenn eine Hubschraube vom Durchmesser d=2r, bzw. der Kreisfläche 
F=d er, einen Schub 5 erzeugen soll, so muß sie Luft nach unten beschleunigen. 


Infolgedessen herrscht bereits in der Schraubenebene selbst eine Durchflußgeschwindigkeit m’. 
Zusammen mit dem Schub S ergibt sieh aus dieser Durchflußgeschwindigkeit die erforder- 
liche Mindestleistung 


Wenn sieh nun unterhalb der Hubschraube ein Boden befindet, so wird die nach unten ab- 
strömende Luft, der Schraubenstrahl, abgelenkt, und der Zustrom der Luft zum Propeller 
behindert. Dadurch wird bei festgehaltenem Schub die Durchflußgeschwindigkeit mw’ geändert, 
womit sich auch die erforderliche Leistung L,- S w’ entsprechend ändert. Insbesondere 
arbeitet die Hubsehraube in der Nähe des Bodens im allgemeinen günstiger als in größerer 
Höhe. Über die Größe dieses Bodeneinflusses soll die nachstehende Überlegung eine Abschätzung 
zeben. Für die beiden Fälle, wenn die Entfernung a des Propellers vom Boden klein oder 
eroß gegen den Propellerradius r, ist, lassen sich die Verhältnisse gut übersehen. 


Wir betrachten den Einfluß des Bodens zunächst unter der Voraussetzung, daß der 
Schub konstant gehalten wird, was sich dureh passende Änderung der Drehzahl sicher erreichen 
läßt. An sich würde sich sonst der Schub infolge der Änderung der Durehflußgeschwindigkeit 
bei Annäherung an den Boden ändern. Diese Anderung des Schubes bei konstanter Drehzahl 
wird im 4. Kapitel besonders behandelt. 


Für unsere Überlegungen wollen wir voraussetzen, daß der Propellerschub S gleichmäßig 
über die ganze Propellerkreistläche F verteilt ist und der Drall des Schraubenstrahles ver- 
nachlässigt werden kann. Dann erfährt jedes Teilchen beim Durchtritt durch den Propeller 
einen Druckzuwachs 

N 
Wenn das Teilchen dann in größerer Entfernung vom Propeller wieder in ein Gebiet ungestörten 


„w, um. Die 


"lüssigkeit, welehe den Propeller durehströmt hat, fließt daher mit einer Geschwindigkeit mw, ab, 
für welehe 


Druckes kommt, so setzt sich diese Drucksteigerung in kinetische Energie 


ist. Am Rande des Propellerstrahles herrscht überall der Druck der umgebenden Luft, und 
da wir diesen mit hinreichender Genauigkeit konstant gleich dem Druck in großer Entfernung 
vom Propeller setzen können, so ist die Geschwindigkeit auf der ganzen Oberfläche des 
Strahles »°,, und nur im Innern des Strahles herrscht in der Nähe des Propellers ein Über- 
und damit eine Geschwindigkeit < 


2. Propellerabstand vom Boden klein gegen Propellerradius. Für einen Propeller in der 
Nähe des Bodens ist diese Abströmgeschwindigkeit »w, waagerecht radial nach außen gerichtet 
(Abb. bi. Ist die Entfernung vom Boden sehr klein gegen den Propellerradius, so spielt sich 


Propelierkreisfläche 


h 


Abb. 


ie Umsetzung von Druck in Geschwindigkeit in unmittelbarer Nähe des Randes der Propeller- 
fläche ab, ın einem Gebiet, dessen radiale Erstreekung klein gegen den Propellerradius ist. 
Infolgedessen können wir bei der Betrachtung dieses Vorganges den Propellerrand ohne großen 
Fehler als geradlinig (anstatt kreisförmig) ansehen. Wir wenden nun auf die durch den 
Propeller strömende Flüssigkeit (es findet nur in der Nähe des Randes eine merkliche Durch- 


lech. 


et 
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trömung statt) eine Impulsüberlegung an. In waagerechter Richtung wirkt der zwischen dem 
Propeller und dem Boden herrschende Überdruck p,. Ist der Abstand des Propellers vom 
Boden a und betrachten wir eine senkrecht zum Propellerrand stehende Schicht von der 
Dieke Öd. so wirkt der Überdruck w auf eine Fläche ad. Er übt also eine Kraft 


aus. Die Zuströmung zum Propeller bzw. zum Rande des Propellers ist von allen Seiten 
vleichmäßig. Ein waagerechter Impuls ist daher nieht vorhanden. Die aus dem Überdruck 
sich ergebende waagerechte Kraft muß daher gleich sein der Reaktion der abströmenden 
Jüssigkeit. Ist die Durchflußmenge durch den Propeller 9, so trifft auf ein Randstück von 
der Länge ö die Menge 


Der Impuls dieser abströmenden Menge ist 


Durch Gleiehsetzen von J und P erhalten wir unter Verwendung von Gl. (3) die sekundliehe 
Durchflußmenge 


28 


Ist die Höhe des abfließenden Strahles im Abstand r von der Propellerachse h, so ergıbt 
sich aus 


durch Vergleich mit Gl. (9) | 


Um die Menge O0 sekundlich gegen den Überdruck p, durch den Propeller zu befördern, 


ist die Leistung 
IX 
0 Po d oF 


erforderlich. Für den ungestörten Propeller ist die Endgeschwindigkeit im Schraubenstrahl 
ebenso groß wie die Abflußgeschwindigkeit des durch den Boden gestörten Schraubenstrahles, 


nämlich 
/ 2 Ss 
oF 


(Gl. (3)). Ferner ist nach den bekannten Überlegungen der einfachen Schraubenstrahltheorie') 
die Durehflußgeschwindigkeit in der Schraubenkreisebene für den ungestörten Propeller 


und die Durchflußmenge 
Die erforderliche Minimalleistung ist 


!) Hütte, I. Bd., 26. Aufl., S 414 oder Handb. d. Phys., VI. Bd., S. 264. 
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Der Einfluß des Bodens bei unverändertem Schub läßt sich demnach durch das Verhältnis 


() 


in einfacher Weise zum Ausdruck bringen. 


3. Propellerabstand vom Boden groß gegen Propellerradius. Die wesentliche Wirkung 
eines Propellers auf die umgebende Flüssigkeit in größerer Entfernung besteht in dem An- 
saugen der Flüssigkeit. Sie läßt sich rechnerisch als Feld einer Senke darstellen, in welche 
die vom Propeller angesaugte Flüssigkeitsmenge verschwindet. Für eine Hubschraube, d. i. 
für einen Propeller ohne Fahrtgeschwindigkeit, ist die ausgesaugte Menge gleich der dureh 


den Propeller strömenden Menge 0, (Gl. 12). Demnach ist die im Abstand 7 vom Propeller 
| £ 


erzeugte Zusatzgeschwindigkeit 


(15), 


sıe ist auf den Propeller zu gerichtet. Befindet sich der Propeller in der Nähe einer unendlich 
ausgedehnten Begrenzungsebene, aber in soleher Entfernung, daß die Störung durch diese 
Ebene nur klein gegen die Durehströmgeschwindigkeit ist, so können wir die Wirkung dieser 
Ebene «durch Spiegelung des Propellers und Ersatz des gespiegelten Propellers durch eine 
Senke rechnerisch erfassen. 


Befindet sieh die Hubschraube im Abstand a über dem Boden, so ist ihr Spiegelbild 
bzw. «die es ersetzende Senke im Abstand 2a unter ihr (Abb. 2). Die Ergiebigkeit der Senke 


! 
d 
Abb. 2. Abb. 3. 


ist 9, (Gl. (12). Sie ruft am Ort des Propellers eine senkrecht nach unten gerichtete 


() F N 
(16) 


Iat2a® 16aaV 20F 


hervor, welehe sieh der ungestörten Durehflußgeschwindigkeit mw’. überlagert. Die resultierende 
ist demnach 


NS IF 


Da sieh bei festgehaltenem Schub die Leistungen wie die Durchflußgeschwindigkeiten ver- 
halten, so ıst 

| | | (18). 


64 


Während bei geringem Abstand des Propellers vom Boden die Antriebsleistung nach Gl. (14) 
verringert wird, tritt bei großem Abstand nach Gl. (15) eine, wenn auch sehr geringe, Er- 
höhung eın. 


= 
| 
| 
| 
| 
| 
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Fassen wir die Ergebnisse der beiden Überlegungen zusammen, so können wir den 
Einfluß des Bodens auf die erforderliche Propellerleistung bei konstantem Schub angenähert 
etwa dureh die in Abb. 3 angegebene diek ausgezogene gebrochene Linie darstellen. Die ım 
vorstehenden ermittelten Abhängigkeiten nach Gl. (14) und (185) sind dünner gezeichnet. Die 
wirkliche Einflußfunktion dürfte sich von der gekniekten Linie nur dureh eine Abryndung der 
Keke unterscheiden’). 


{. Einfluß des Bodens auf den Schub. Bei den vorstehenden Überlegungen hatten 
wir immer vorausgesetzt, daß der Schub bei Annäherung an den Boden konstant gehalten 
werden soll. Da sich dabei die Durehflußgeschwindigkeit durch den Propeller ändert, so 
kann der konstante Schub nur durch gleichzeitige Änderung der Drehzahl oder durch 
Änderung der Propellerform, z. B. des Steigungswinkels, aufrechterhalten werden. Von einer 
Formänderung wollen wir aber abschen. Ist für den Propeller außer dem Schub am Stan«dl 
auch die Abnahme des Schubes bei kleinem Fortschrittsgrad r/u bekannt Fortschritts- 
veschwindigkeit, « = Umfangsgeschwindigkeit), so können wir den Schub in Bodennähe bei 


N 

beliebiger Drehzahl angeben. Wir wollen die Abhängigkeit der Schubziffer kr, ” 

> Fu‘ 
(S  Propellersehub, Luftdiehte, Propellerkreistläche) linear annehmen 

wobei die Schubziffer in unbegrenzter Luft am Stand (4 0) bedeutet und = 
(A 


konstant sein soll, was bei den kleinen Steigungen der Hubsehrauben einigermaßen zutrifft. 


Ist die effektive Steigung der Schraube MH, so ergibt sieh für 4 = eg der Schub Null, so daß 


ok, _, da 


ist. Ist bei konstantem Schub die Durchtlußmenge 9 durch den Boden um 10.0, 0 
vermindert, so ist die mittlere Durchflußgeschwindigkeit um 


F 


Aw=no (21) 


verkleinert. Der Propeller verhält sich daher so, als wenn er eine Fortschrittsgeschwindigekeit 


habe. War zur Konstanthaltung des Schubes eine Herabsetzung der Umfangsgeschwindig- 
keit «, auf «’ nötig, so ist der Fortschrittsgrad 


ET 
A=7 25) 
und deshalb nach Gl. (19) 
Am x» Aw da 
Der Schub ist dann für irgendeine Drehzahl » bzw. Umfangsgeschwindigkeit u 


°) Eine Untersuchung, welehe den Verlauf genauer ermittelt, ist mit Aussicht auf Erfolg in Angriff genommen. 
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Durch Vergleich mit dem Schub 


5 daAw 
H 7) 


\ 


ergibt sich 


Wenn wir uns auf Einflüsse beschränken, bei denen Jw< w’ ist, was praktisch wohl immer 
der Fall ist, so können wir in dem dann kleinen Korrekturglied in der Klammer «’ gegen 
, vertauschen. Da ferner nach Gl. (21) und (11) 


ist, wobei | die ermittelte Einflußfunktion des Bodens bei konstantem Schub ist (Abb. 3), 


können wir den Schub bei beliebiger Umfangsgeschwindigkeit durch diese Funktion ausdrücken 


N da / S, I, 
S, I+ Hu, 20 (29). 
Insbesondere ist bei konstanter Drehzahl «= u, 
N dr S, L, 


5. Berücksichtigungen der zusätzlichen Verluste. Bei Anwendung der vorstehenden Er- 
gebnisse ist zu beachten, daß außer dem Energiebedarf infolge der .Durchflußgeschwindig- 
keit »w’ (Gl. 1) noch ein weiterer Energiebedarf zur Überwindung des Flügelwiderstandes 
und wegen störender Nebeneinflüsse vorliegt. Infolgedessen ist die Leistung bereits bei 
unendliehem Abstand vom Boden L, um einen Betrag Lr„ größer als die theoretische 
Leistung /,, nach Gl. (15) 


Bei fester Gleitzahl der Flügelprofile ändert sich diese zusätzliche Leistung in Bodennähe bei 


festgehaltener Drehzahl », proportional dem Schub. Darnach erhalten wir für diese zusätz- 
liche Leistung bei einer Drehzahl » 


Sn 
32 
(32), 
wobei aus Gl. (29) zu entnehmen und 
Io N, 


Da die verlustlose Leistung Z, = S mw’ also proportional S®? ist, wird die Gesamtleistung 
bei einer beliebigen Drehzahl n 


sn L,\ S/ı/ S 
Hierbei sind die Werte „ und Lp, aus Gl. (29) und (31) zu errechnen und N wieder der 


Kıntluß bei konstantem Schub. 


Zusammenfassung. Auf Grund einfacher Überlegungen wird der Einfluß des Bodens auf 
den Leistungsbedarf einer Hubschraube in kleiner und in großer Entfernung vom Boden 
ermittelt und durch eine einfache geknickte Linie dargestellt. Die gefundene Einflußfunktion 
wird dazu benützt, um auch die Schubänderung in Bodennähe bei beliebiger Drehzahl zu 
ermitteln, und schließlich wird noch gezeigt, wie man den Einfluß der zusätzlichen Propeller- 
verluste berücksichtigen kann. 
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Hodographenmethode der Gasdynamik. 


Von A. Busemann in Braunschweig. 


1. Einleitung. Im Jahre 1590 hat bereits der Holländer P. Molenbroek') die Hodo- 
„raphenmethode der Gasdynamik für das Gas mit dem Verhältnis der spezifischen Wärmen 
ee, 1 entwickelt und mit der Legendreschen Transformation verglichen. Der Russe 
A. Tschapligin?) hat diese Methode auf beliebige Gase erweitert (1904). Beide Arbeiten 
blieben zunächst wenig beachtet, bis sie durch D. Riaboucehinski’) und B. Demtschenko’) 
1932 ans Tageslicht gezogen und durch neuere Betrachtungen vermehrt wurden. 


Unter Hodographenmethode versteht man die Integration einer Strömung mit Hilfe der 
Abbildung ihrer Grenzbedingungen in den Geschwindigkeitsplan. Bei der ebenen Flüssigkeits- 
strömung läßt sich das komplexe Potential F=P+iY konform auf die Ebene des konjugiert 
komplexen Geschwindigkeitsvektors w=u — iv abbilden (vgl. Prandtl-Tietjens, Bd. 1, 
S, 169). Gelingt es, die Grenzbedingungen ohne Kenntnis der Lösung in die Geschwindigkeits- 
ebene zu übertragen, so kann man die Integration hier vollziehen. Freie Strahlgrenzen 
zegenüber einem Raum mit gegebenem konstantem Druck lassen sich auf diese Weise sogar 
leichter behandeln, weil nicht der Verlauf, sondern die konstante Geschwindigkeit längs der 
Grenzstromlinie bekannt ist. Die Abbildung in die Geschwindigkeitsebene verwandelt dem- 
nach die freien Strahlgrenzen in feste Grenzen. Außer diesen lassen sich aber auch feste 
veradlinige Begrenzungen mit Hilfe ihrer Richtungen in die Geschwindigkeitsebene übertragen. 
Aus geradlinigen festen Grenzen und aus freien Strahlgrenzen zusammengesetzte Rand- 
bedingungen ergeben daher feste Grenzen in der w-Ebene. Dies ist der Vorteil der Abbildung, 
denn die Differentialgleichung ist in beiden Ebenen die gleiche. 


2. Die Differentialgleichung der Gasströmung. In der Gasströmung muß man beachten, 
daß sich die Drehungsfreiheit auf den Geschwindigkeitsvektor iv, die Quellfreiheit aber auf 
den Stromdiehtevektor 8 bezieht, der durch Multiplikation mit der veränderlichen Diehte aus 
dem Geschwindigkeitsvektor hervorgeht. Die Diehte o soll sich adiabatisch mit dem Druck p 
verändern, sie wird dann ebenso wie der Druck eine Funktion des Geschwindigkeitsbetrages m: 

\ 2 09 


(N). 


Die angegebene Art der Funktion gilt für vollkommene Gase mit einem Verhältnis der spezi- 
fischen Wärmen e,/e,=x, einer Dichte o, und einer Schallgeschwindigkeit «a, im Ausgangs- 
zustand ==0, Um den Stromdiehtevektor bei kleinen Geschwindigkeiten in den Gesehwindig- 
keitsvektor übergehen zu lassen, also den Anschluß an die Flüssigkeitsströmung zu erhalten, 
wird man folgende Definition vorziehen: 


Die reibungslose stationäre Gasströmung hat dann die Differentialgleichungen: 


) 


In der ebenen Gasströmung mit den Lagekoordinaten x und y und den Geschwindigkeits- 
komponenten « und » folgt aus diesen Gleichungen die Existenz eines Geschwindigkeits- 
potentials ? und einer Stromfunktion 7 mit folgender Bedeutung: 


Zwischen beiden Funktionen bestehen daher die Beziehungen: 


') P.Molenbroek: Über einige Bewegungen eines Gases bei Annahme eines Gesehwindigkeitspotentials. 
Arch. d. Mathem. u. Phys., Grunert Hoppe 1890, Reihe ?, Bd. 9, S. 157. 

”) A. Tschapligin: Wiss. Ann. d. Univ. Moskau (Math. Phys., Bd. 21, 1 bis 121, 1904). 

») D. Riabouchinski: Comptes rendus de l’Ae. d. Se. Paris (194, 1215, 1932). 

*) B. Demtschenko: Comptes rendus de !’Ac. d. Se. Paris (194, 1218 u. 1720, 1932) und Publ. Math. Univ. 
Belgrade (2, 85, 1933). 


(26) 
27). | 
ner 
gen | 
23 | 
3), | 
en 
29). | 
0). 
| 
Ig- 
es 
)ei 
he 
1). 
am 
| | 


- Ztschr. f.angew. 
Busemann, Hodographenmethode der Gasdynamik Math. und 


Differentiiert man die erste Gleichung nach y und die zweite nach ., so ergibt ihre Differenz: 


a” | 


oa — | 


Hierin bedeutet a die Schallgeschwindigkeit, sie ergibt sich bei der Differentiation von Gl. (1) 
ın folgender Weise: f’ = do/dw= —- owf/a?’, Löst man die Gl. (5) nach Y auf, so erhält man 
(ureh Differentiation die entsprechende Differentialgleiehung für das Potential: 


a a? 


In dieser Form erscheinen beide Differentialgleichungen identisch. Da jedoch « und » selbst 
wieder Ableitungen und zwar verschiedene für ?P und Y nach Gl. (5) sind, so handelt es sich 
in Wahrheit um zwei verschiedene nicht lineare Differentialgleichungen zweiter Ordnung. 
Nach Riemann (1860) kann man diese Differentialgleichungen durch Legendresche Trans- 
[formation aus der in die bzw. die entsprechende yqr-Ebene der Strom- 
(liehte übertragen (vgl. Handbuch d. Physik, Bd. 7. S. 315). Man erhält dann in den trans- 
formierten Funktionen 


lineare Differentialgleichungen zweiter Ordnung 


nit Aa 0/0,. Die Grenzbedingungen lassen sich nieht mehr vor der Lösung in 
die neue Ebene übertragen. Nach der Lösung findet man die Orte der „Ebene dagegen 
sehr leicht dureh Differentiation: 


Nur wenn etwa in der „Ebene eine Lane mit festen Geschwindigkeitsvektoren gegeben 
wäre, könnte man diese Bedingung vor der Lösung in beide Ebenen übertragen (vgl. z. B. die 
\ever-Strömung dureh eine Lavaldüse, Forschung. d. Ing.-Wes. 1908). Die transformierten 
Funktionen y und y haben mit dem Potential und der Stromfunktion nicht mehr viel zu tun. 
Beide Möglichkeiten sind vielmehr völlig gleichwertige Darstellungsarten. Man wird jedoch 
das transformierte Potential in der wr-Ebene vorziehen, weil hier Unterschall- und Über- 
schallbereich auseinandergebreitet erscheinen, während sieh diese Bereiche in der Ebene der 
Stromdiehten überdecken. 


3. Die Umformung von Molenbroek und Tschapligin. Im Gegensatz zur Legendreschen 
Transformation werden jetzt das Geschwindigkeitspotential P und die Stromfunktion Y selbst 
auf die Gesehwindigkeitsebene abgebildet. Da die günstigsten Verhältnisse in verzerrten 
Geschwindigkeitsebenen erhalten werden, empfiehlt es sich, die logarıthmische konforme Ab- 
bildung der Gesehwindigekeitsebene mit den kartesischen Koordinaten In ı= und dem Stromlinien- 
winkel 5 zugrunde zu legen. Diese läßt sich leichter dureh Verzerrung der Abszisse beliebigen 
Funktionen des Gesehwindigkeitsbetrages anpassen als die polare Darstellung der Geschwindig- 
keit. Da der Logarithmus von »w nur beim Übergang zur Flüssigkeitsströmung nötig ist. 
sollen Abszisse und Ordinate zunächst mit ww und 5 bezeichnet werden. 


Die Umformung geht von Gl. (4) aus, die jetzt in folgender Form zu schreiben ist: 


dP w-cosp- de+m-sinp-dy, „m cos 
Oo Oo 


Nach dr und dy aufgelöst bedeutet dies: 


cos sin I0 
m 
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nz: Soweit sich die Gasströmungen nach den Geschwindigkeitsverteilungen umkehrbar eindeutig 
auf die Geschwindigkeitsebene abbilden lassen, müssen sich = und y als Funktionen von 
und 5 schreiben lassen, und dann sind auch ® und Y Funktionen von w und P. 
(6). 
Hieraus folgt für (12): 
Ibst Außerdem müssen folgende Bedingungen erfüllt sein: 
IST 
Ins- 
Diese Bedingungen enthalten nur erste Ableitungen für das Potential und die Stromfunktion: 
(S) 
Da in beiden Gleichungen P’.„ und Y’5 sowie P’5 und Y”,. mit proportionalen Faktoren vor- 
kommen, erhält man bei der Auflösung nach 2’. und 2’ jeweils nur eine erste Ableitung 2 
von 
2 


Beim Verzerren der Abszisse ww ändert sich das Differentialgleichungspaar so, daß der Faktor 
10) bei 7’5 mit einer Funktion von w multipliziert und zugleich der Faktor bei Y",. dureh die- 
selbe Funktion dividiert wird. Man erhält demnach bis auf das Vorzeichen den gleichen 


a: Faktor in beiden Gleichungen, wenn man folgende Verzerrung vornimmt: 
die 
ten 
‚dw diom) dw 
un. de | | (15). 
on A n 
Es lautet dann das Gleichungspaar: 
ler 

hst Diese eindeutige Wahl der Abszisse @ liefert ein Gleichungspaar, das in gleicher Form 
ten zwischen dem elektrischen Potential und der elektrischen Stromfunktion gilt, wenn man einen 
Ah. Elektrolyten in einem flachen Gefäß von einem elektrischen Strom durchfließen läßt und am 
| 0, m” 
ven Orte der Oberfläche eine Tiefe T umgekehrt proportional der Funktion | 
. 
lıo- . r . op .. 
ww vorsieht. Vertauscht man dagegen die elektrischen Größen gegenüber den entsprechenden 
st, 


Begriffen der Gasströmung, indem man die Stromfunktion durch das elektrische Potential 
und das Geschwindigkeitspotential durch die elektrische Stromfunktion darstellt, so wird die 
Tiefe direkt proportional dieser Funktion. Das Gefäß müßte also in jedem Falle einen 
zylindrisch gewölbten Boden erhalten, der bei direkter Analogie beim Erreichen der Schall- 
1) geschwindigkeit eine unendliche Tiefe ergibt, bei der inversen Analogie eine mit Annäherung 
an die Schallgeschwindigkeit abnehmensle Tiefe erzeugt. Diese elektrische Strömung bleibt 
bei einer konformen Abbildung der Oberfläche erhalten. Man kann daher auch in Polar- 
darstellung mit der verzerrten Geschwindigkeit e” +77 arbeiten. Die notwendige Tiefe T bleibt 
erhalten und ordnet sich auf konzentrischen Kreisen an, so daß der Boden eine achsen- 
2). symmetrische Fläche wird. Diese Polardarstellung ist für die Anwendung gut geeignet, denn 
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der ganze Unterschallbereich bildet dann nur eine begrenzte Fläche. Wie man aus Gl. (15 
erkennt, ist der logarithmische Zuwachs des Radius das geometrische Mittel zwischen den 
logarıthmischen Geschwindigkeits- und Stromdichtezunahmen. Der Radius bleibt demnaelı 
endlich, weil sowohl die Geschwindigkeit als auch die Stromdiehte begrenzt sind. Er liegt 
zwischen Geschwindigkeit und der dazugehörigen Stromdichte, die um den Faktor o/o, kleiner 
ist. Der Überschallbereich läßt diese elektrische Analogie leider nicht zu, die Tiefe der Schale 
würde imaginär. Dieses Ergebnis sagt nur, daß die elektrische Strömung mit ihrer elliptischen 
Differentialgleichung nicht geeignet ist, für die Integration einer hyperbolischen Differential- 
gleichung verwendet zu werden. 

Kliminiert man das Strömungspotential aus dem Gleichungspaar (14), so ergibt sich eine 
lineare Differentialgleichung zweiter Ordnung für die Stromfunktion. Man erhält nur zweite 
Ableitungen, wenn man die Abszisse vorher so verzerrt, daß der Faktor auf der rechten Seite 
der zweiten Gleichung von mw unabhängig wird. Auf jeden Fall erhält man eine wesentlich 
einfachere Differentialgleichung in der Geschwindigkeitsebene als in der »=y-Ebene. Der 
Vorteil der Hodographenmethode ist demnach größer als bei der Flüssigkeitsströmung; dort 
blieb die Differentialgleichung vom gleichen Typ, während sie hier bei der Abbildung line- 
arisiert wird. Grundsätzlich wird jeder für die Flüssigkeitsströmung lösbare Fall durch das 
elektrolytische Verfahren bei beliebigem Adiabatengesetz des Gases gelöst. 


4. Die „langsame“ Gasströmung. Tschapligin entdeckte auch ein Adiabatengesetz, 
für das die Differentialgleiehung in der Geschwindigkeitsebene die Differentialgleichung der 
Potentialtheorie wird. Bei diesem Gas bleibt die Tiefe des Elektrolyten konstant. Analytisch 
kann man hier Stromfunktion und Potential durch konforme Abbildung bestimmen. Es ist 
dies das Gas mit einer geradlinigen Adiabaten im Druck-Volumen-Diagramm. Da die wirk- 
lichen Gase eine gekrümmte Adıabate besitzen, so ist diese Vereinfachung nur dort von prak- 
tischem Wert, wo sich die Adiabate durch ihre Tangente annähern läßt. In der gewöhnlichen 
Aerodynamik benutzt man überhaupt nur einen Punkt der Adiabaten, nämlich o = o, = konst., 
ohne daß dadureh der brauchbare Geschwindigkeitsbereich verschwindend klein wäre. Wieviel 
erößer muß daher der zulässige Geschwindigkeitsbereich ausfallen, wenn man die Adiabate 
dureh ihre Tangente ersetzt und damit endliche Diehtedifferenzen berücksichtigen kann. Nur 
um einen Anhalt zu bieten, sei ausgesprochen, daß der Geschwindigkeitsbereich etwa zwischen 
der 0.2-fachen und der O,5-fachen Schallgeschwindigkeit liegt, bei dem nieht mehr die Annahme 
konstanter Dichte, wohl aber die Adiabatentangente ausreicht. Spricht man erst dann von 
einer Gasströmung, wenn die Volumenänderung Einfluß auf die Strömung bekommt, so kann 
man dies Gebiet mit Demtschenko als die „langsame* Gasströmung bezeichnen. 


Bei vollkommenen Gasen erhält der Faktor in der Gl. (16) folgenden Wert (vgl. auch (1)): 


— 
| ) 2 
Kr ist demnach veränderlich außer für 2° — 1. Unter Berücksichtigung der Tatsache, daß 


man der Strömung einen hydrostatischen Druck p, überlagern darf, kann man nämlich auch 
(die geradlinie Adiabate im Druck-Volumen- Diagramm in der Form (p, — = konst. schreiben. 
Für die Diehte gilt dann (p, pl=C-o°'. Die Konstante © läßt sich durch die Schall- 
geschwindigkeit berechnen, denn es gilt dp/do = a’. Hieraus folgt = a,” für die Adiabaten- 
tanzente beim Zustand o,. Die Strömung für ein Gas mit = — 1 läßt sich in der Ebene mit 
den Koordinaten » und 5 und jeder aus ihr dureh konforme Abbildung entstandenen Ebene 
dureh komplexe Potentiale wie für eine Flüssigkeitsströmung lösen. Wesentlich ist die 
Verzerrung dieser Hodographenebene gegenüber der Geschwindigkeitsebene, die dureh Gl. (15) 
gegeben ist. Für ein beliebiges vollkommenes Gas liefert die Integration von (15): 


W= a, 
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Für = — 1 erhält man also: 
2a, m 
ISa 
(15a) 
Die Umkehrung dieser Beziehung lautet: 
W 
9). 
4a, — W° ( 


Der für die Integration der Gl. (12) außerdem notwendige Faktor o,/o bekommt bei diesem 
Gas (vel. Gl. mit = —1) folgenden Ausdruck in W: 


w” 4a, 


dar 


9% 


Kennt ıman demnach das Potential ? und die Stromfunktion 7 in der W-Ebene mit den 
Koordinaten U= W-.cosß und V= W-sin 5, so ist die Übertragung der gefundenen Lösung 
dureh die Gl. (12) in Verbindung mit (19) und (20) in die Koordinaten = und y möglich: 


Die Funktionen P (U, V) und Y(U, V) sind durch die Grenzbedingungen und die Singularitäten 
in der W-Ebene bestimmt. Sie werden daher in den meisten Fällen mit Funktionen P (uw, ®) 
und Y(w,v) übereinstimmen, die für eine gleichartige, nur in den maßgebenden Bestimmungs- 
stücken etwas veränderte Flüssigkeitsströmung gelten. Auf der gesuchten freien Strahl- 
erenze ist nun Y’=konst., also d und es ist W=konst. Berücksichtigt man diese 
Umstände bei der Integration der Gl. (21), so ergeben sich in der betreffenden Flüssigkeits- 
strömung und in der gesuchten Gasströmung geometrisch ähnliche Strahlgrenzen. Der Ver- 
kleinerungsfaktor jeder einzelnen Begrenzungskurve beträgt (1 - W°’/4 a,’). Dieses Gas behilft 
sich bei den freien Strahlgrenzen mit den Formen der Flüssigkeitsstrahlen, indem diese für 
etwas andere Strahldicken, andere Wandabstände oder andere Geschwindigkeitsverhältnisse 
benutzt werden. Die Gase mit gekrümmten Adiabaten sind nur bei kleinen Geschwindigkeiten 
ebenso bescheiden, bei höheren Geschwindigkeiten macht sich die veränderte Differential- 
geleichung bemerkbar, das elektrolytische Verfahren zur Darstellung der Potentiale liefert 
entsprechend der Bodenwölbung des Gefäßes ganz neuartige Verteilungen der Potentiale und 
damit andere Strahlformen. Bei der Auswertung genügt es entweder, die Potentiale oder die 
Stromfunktionen zu kennen, weil die in Gl. (21) fehlenden Werte durch die Beziehungen (16) 
zu ermitteln sind. 


5. Strömung durch eine Blende’). Um das Hodographenverfahren auch in den Einzel- 
heiten übersehen zu können, soll als einfaches Beispiel die Strömung durch eine Blende 
behandelt werden. Ein rechteckiger Kanal von der Breite 24, ist am Ende durch eine Blende. 
mit der Schlitzbreite 2, verengt. Die zweidimensionale Ausflußströmung mit Strahlbildung 
durch diese Blende soll bestimmt werden (Abb. 1). Den Ausfluß aus einem Kessel erhält 
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j 
U 
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Abb. 1. Abb. 2. 


usemann: Die Expansionsberichtigung der Kontraktionsziffer von Blenden. Forschung Bd. 4, 
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ınan, wenn , gegen unendlich geht. Auf allen Stromlinien herrscht im Kanal in großer 
Entfernung vom Ende die Geschwindigkeit «,. im Strahl gelangen die Stromlinien in großer 
Entfernung von der Austlußöffnung wieder zu einer einheitlichen Geschwindigkeit «,, die 
dureh des Druckgefälle zwischen dem Kessel und dem Außenraum bestimmt ist. Aus dieser 
Tatsache ergibt sieh ohne Kenntnis der gesamten Lösung, daß im Hodographen (Abb. 2) der 
Bildpunkt 7, für die Geschwindigkeit «, die Quelle aller Stromlinien und der Bildpunkt 7, 
für die Gesehwindigkeit #, die Senke aller Stromlinien darstellt. Die Werte U, und U, sind 
nach Gl. (15) zu ermitteln: 


2 a,u, 2a,u, 55, 

Die beiden begrenzenden Stromlinien laufen in Abb. I zunächst entlang der Kanalwand bis 
zum Staupunkt in der Eeke zwischen Wand und Blende. Auf den begrenzenden Stromlinien 
zibt es «daher auf diesem Stück nur die "Komponente der Geschwindigkeit mit Werten 
zwischen #, und Null. Die Kanalwände bilden sich daher auf der T-Achse zwischen ÜU= TU, 
und 2=0 ab. Vom Staupunkt aus gelangen die äußeren Stromlinien entlang der Blende bis 
zur Blendenöffnung. Hier besteht nur eine »-Komponente der Geschwindigkeit. Daher wird 
(lie benetzte Blendenwand dureh die V-Achse innerhalb des Kreises U? + V?’= 1,’ wieder- 
vergeben. Die freien Strahlgrenzen sind dureh die Bedingung konstanter Geschwindigkeit 
auszezeiechnet, daher liegen sie auf dem Kreis, und zwar sind es die beiden Bögen von der 
V-Achse bis zum Punkt 7, der !- Achse. Der ganze Strömungsbereich zwischen diesen beiden 
äußeren -Stromlinien bildet im Hodographen die Halbkreisfläche rechts von der V-Achse, sie 
ist aber auf der U-Achse zwischen 7-1, und Ü=0 geschlitzt vorzustellen. Die beiden 
Ränder des Schlitzes bedeuten die zusammenfallenden Geschwindigkeitswerte entlang der 
oberen und unteren Kanalwand. Die „-Achse der Abb. 1 ist in ihrem ganzen Verlauf dureh 
das Stück der !-Achse zwischen 7, und 7, dargestellt, dies ist zugleich die mittlere Strom- 
linie 7==0. Diese Bedingungen lassen sich sämtlich erfüllen, wenn man die Quelle in U, 
und die Senke in 7, sowohl an der V-Achse als auch am Kreise spiegelt. Die Spiegelung 
am Kreise ergibt eine Verdoppelung der Senke in 7, und eine weitere Quelle n U, = U/U,. 
Die Spiegelung an der V-Achse liefert dieselben Singularitäten auch auf der negativen Ü-Achse. 
/ur Vereinfachung der Rechnung soll im Gegensatz zur früheren Festlegung mit W=U-+iV 
die komplexe Koordinate und mit WU @V die konjugiert komplexe Koordinate der Hodo- 
eraphenebene bezeichnet werden. sowie konjugiert komplex P- 
ist das komplexe Strömungspotential. Mit der Quellstärke 9 lautet dies: 


FIW)= U)—2ln(W- . (23). 


Die Übertragung dieser Lösung in die Koordinaten . und y läßt sich nach Gl. (21) auch komplex 
in 2. schreiben: 


W-dF 
d W | . . . . . . . (21 a). 


Die Integration läßt sich bei beiden Summanden leicht für jede einzelne Quelle oder noch 
besser für ein symmetrisch zur V-Achse gelegenes Quellenpaar durchführen. Für alle drei 
Quellenpaare erhält man schließlieh: 


9, W-D, Qu, W-U, 9, W-U 


W-T, 207, W--U, 
W+L, la; n 


1 


(24). 


Will man nur «die Hauptabmessungen und deren Verhältnis kennenlernen, so genügt es, die 
Koordinaten 7, und zu bereehnen. Die Koordinate heißt allgemein für einen Punkt 
unter Benutzung der dureh ihn festgelegten Winkel &,. +, und &, (vgl. Abb. 3) bei einer Quell- 
stärke U 
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Diese Winkel sind bereits so gezählt, daß sich die «-Achse (y=0) auf das Stück der U-Achse 
wischen U= U, und U=T, abbildet. Schiebt man den Punkt P in Abb. 3 auf die positive 
[’.Achse links von U,, so erhält man die Höhe „, der oberen Kanalwand. Hier ergibt sieh 
a und 28,0. Daraus folgt nach (24a): 


. (24b). 


[A6833] 


Abb. 3. 


Die Höhe 7, der Blendenkante erhält man dureh deren Bildpunkt W=i T,.auf der V-Achse. 


Dort ist &, =2 arctg (U — e, und —r/2. Demnach: 


Schließlich findet man die halbe Breite des ausströmenden Strahles ,, wenn man auf dem 


Kreise an den Punkt U, heranrückt. Dabei wird dann &,=r,-0 und »,: ni2: 
u? 
ta, 


Beim Ausfluß aus einem Kessel geht „, gegen unendlich, d. h. die Geschwindigkeit «, und 
damit U, geht gegen Null. In diesem Falle erhält man aus Gl. (24e) und (24d) mit einer 
Umformung nach Gl. (20) die folgende Kontraktionsziffer: 


4a, +U, 


Jede einzelne Strahlgrenze ist in diesem Fall genau die geometrisch ähnliche Verkleinerung 


der Strahlgrenze beim Ausfluß einer Flüssigkeit. 


6. Zusammenfassung. Die Übertragung der Hodographenmethode von der Flüssigkeits- 
strömung auf die Gasströmung ist von Molenbroek 18% angegeben. Sie erweist sich bei der 
(rasströmung mit Unterschallgeschwindigkeit noch wesentlich vorteilhafter im Vergleich zu 
den übrigen Berechnungsverfahren der Hydro- und Gasdynamik, weil mit der Übertragung 
der Strömung in die Hodographenebene zugleich eine Linearisierung der Differentialgleichung 
verbunden ist. Für ein Gas mit geradliniger Adiabate ist die Hodographenmethode nach 
Fschapligin sogar eine Zurück führung auf die Differentialgleichung der Flüssigkeitsströmung 
innerhalb der Hodographenebene. (Die von Tsehapligin für vollkommene Gase entwickelten 
eihen konnten in dieser Arbeit nicht berücksichtigt werden. Vgl.) oder bei *) die letzt- 
genannte Arbeit.) Nach Molenbroeks Worten wäre die Linearisierung durch Legendresche 
Transformation „wegen des natürlicheren und kürzeren Weges der Ableitung vorzuziehen, allein 
die Nebenbedingung an der Strahloberfläche nimmt bier eine verwickeltere Gestalt an. Es 
ist dies die Ursache, wenn doch die Hodographenmethode in den Anwendungen wahrscheinlich 
den Vorteil davonträgt.“ 65) 


| 
| W=U-iV 
| 

| 
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Die zweite Invariante des Verzerrungstensors. 


Von M. Lagally in Dresden. 
1. Bezeichnungen. Die Navier-Stokesschen Gleichungen für die Strömung einer zähen 
Flüssigkeit lassen sich in vektorieller Form!) in die eine Gleichung 


u 


dt 


zusammenfassen. Dabei ist 
die Geschwindigkeit; 
p der isotrope (hydrostatische) Druck; 
o die Diehte (sie kann vom Ort und der Zeit abhängen); 
\" das Potential der äußeren Kräfte, bezogen auf die Masseneinheit; 
« der Koeffizient der inneren Reibung. 


Für «0 geht (1) in die Eulersche Gleichung für reibungslose Flüssigkeiten über: 


dv 
dv. 
H ist «lie substantielle Beschleunigung; sie zerfällt nach 


in einen lokalen Teil \} und einen konvektiven Teil v-[ dv. 
( 


Neben (1) bzw. (la) besteht die Kontinuitätsbedingung 
do 
dt 


2. Der Verzerrungstensor und seine Invarianten. In (2) tritt eine Dyade | v auf, die als 
Verzerrungstensor (zenauer Tensor der Verzerrungsgeschwindigkeit) zu bezeichnen ist. In 
rechtwinkligzen Koordinaten ist 


dabei ist +d Jo die substantielle Anderung der Dichte. 


(4). 


Man erhält in folgender Weise ein vollständiges System von (algebraischen) Invarianten des 
Verzerrungstensors. Eine Invariante 3. Grades ist die Determinante von Y d 


dr dm 
ou dr dw 
D). 


ou dr dw 
02 02: dz 


Je eine Invariante 1. und 2. Grades wird von der Summe der Diagonalglieder und von der 
Summe der Diagonalminoren von J, gebildet: 


Or 


J.= + = =. . 6): 
dw ou Om ou Or 
04 dx 0x 
or dw ou ou 


', Wegen der Bezeichnung vgl. man: M. Lagally: Vorlesungen über Vektorrecehnung (1928). 
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Außerdem gibt es noch eine vektorielle Invariante 


| ‚(dw ov (du om (2° 


(lie man nötigenfalls durch 3 skalare Invarianten ersetzen kann, die, wenn [ d symmetrisch 
ist, gleichzeitig mit 1 X dv verschwinden; darunter befindet sich eine Invariante 2. Grades, nämlich 
rotd ? | 


0 3 02 0%, 


Damit ist ein vollständiges System von Invarianten gewonnen. Man kann natürlich die 
gewählten Fundamentalinvarianten durch andere ersetzen; so kann man an Stelle von J, die 
>. Verjüngung von (V vd)’ oder das doppelt skalare Produkt 


setzen, das mit J, dureh die Gleichung 
Lv - J,’ 2 . . . . . . . . . (10) 


zusammenhängt. 

Von diesen Invarianten sind J, = divd und rotd sehr bekannt und mechanisch leicht 
zu deuten: auch J, ist einer einfachen Deutung zugänglich: Das Maß der Verzerrung, die die 
strömende Flüssigkeit in der Zeiteinheit in der Umgebung U’ eines Punktes /’ erleidet (Ver- 
zerrungsgeschwindigkeit), ist der Unterschied der Geschwindigkeit in benachbarten Punkten 
vd. Diese Gleichung vermittelt eine affıne Abbildung der durch den Ortsvektor dr 
bestimmten Umgebung U auf ein Bildfeld dv (die Verzerrung von U). Dann ist J, der 
()uotient aus dem Volumen des Bildfeldes von U und dem Volumen von U; dieser (Quotient 
ist nieht mit der Volumdilatation zu verwechseln’), 

Besonders einfach ist die Bedeutung von J,--0; dann ist J d eine ausgeartete (planare) 
Dvade, und die Verzerrungsvektoren dv aller Ortsvektoren dr von U sind einer festen Ebene 
parallel. 

3. Die Dissipationsfunktion®). Eine Invariante 2, Grades spielt eine wichtige Rolle in 
der in der Theorie der reibenden Flüssigkeiten auftretenden Dissipationsfunktion; die Energie, 
welche infolge der Reibung in Wärme übergeführt wird, also mechanisch verloren geht, ist, 
bezogen auf die Masseneinheit 


Kine Umformung ergibt 


> u 


D= 


| 


-- 


2u 


| 
D |rotd? 


| 


Für reibungslose Flüssigkeiten verschwindet D; aber nicht etwa dadurch, daß eine Invariante 
2. Grades verschwindet, sondern dadurch, daß «0 wird. 


4. Der Satz von Hamel’) und seine Erweiterung. Vor kurzem hat Hamel einen Satz 
über den Druck einer strömenden Flüssigkeit angegeben, der eine Invariante 2. Grades ent- 
hält und für reibende und reibungslose Flüssigkeiten in gleicher Weise gilt, der aber aller- 
dings Volumbeständigkeit voraussetzt; der Satz besteht in der Gleichung 

Dieser Satz läßt sich für kompressible Flüssigkeiten in folgender Weise erweitern. 
Setzt man inden Navier-Stokesschen Gleichungen zur Abkürzung +V=P: =R, 

wobei P die Summe aus der Kompressionsenergie und der potentiellen Energie der äußeren 
Massenkräfte, R die Reibungskräfte, bezogen auf die Masseneinheit, bedeutet, so wird (1) 
dv 
At 
”) G.Hamel: Elementare Mechanik, (1912), S. 566. 
») H. Lamb: Hydrodynamik, deutsch von J. Friedel, (1907), 8. 668. 


ei G. Hamel: Ein allgemeiner Satz über den Druck bei der Bewegung volumbeständiger Flüssigkeiten. Monats 
hefte für Mathematik und Physik, 43 (1936), S. 349. 
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Bildet man in dieser Gleichung beiderseits die Divergenz, so ergibt sich 


dv 
»- 
Dabei ist Vo vo Hier tritt in der 
Summe aus dem 1. und 3. Summanden die konvektive Anderung von J, =] »vd auf: ArHaR 
1.J 
ersetzt man außerdem vd nach (10) durch 2,J,, so wird 
dv 
1.J 
dt 


Das ist die Erweiterung des Hamelschen Satzes. Eine unmittelbare mechanische 
Deutung scheint der erweiterte Satz ebensowenig zuzulassen wie der ursprüngliche; was ihn 
interessant macht, ist das Auftreten der Invariante J, in einer für Jede Art von Flüssigkeits- 


‚fl 
bewegung geltenden Gleichung. Mit I -R = ul (or) kommt ein fremdes Element in die 


Gleichung herein; FR ist, wie übrigens auch Iv, als Differential-Invariante von | d aufzu- 
fassen. Zu bemerken ist, daß I -R in zwei wesentlich verschiedenen Fällen verschwindet; 
nämlich, wenn «#«=0 ist (Fall der reibungslosen Flüssigkeit) und wenn 1 -v=,J,—=0 ist 
(Hamelscher Fall der inkompressiblen Flüssigkeit). 


5. Invarianten des Richtungsfeldes der Geschwindigkeit. Führt man an Stelle der Ge- 
schwindigkeit dvd den Einheitsvektor ihrer Richtung t=,,, ein, so kann man die Invarianten 


von insbes. (ft), und (J t), geometrisch deuten. 


Im allgemeinen besitzt eine Kurvenkongruenz mit dem Tangenteneinheitsvektor t keine 
Orthogonalllächen: das gleiche gilt von den beiden Kurvenkongruenzen, deren Richtungen 
durch die Hauptnormalen n und die Binormalen b der Kurvenschar t gegeben sind. 


Die Geometrie eines dreifach-orthogonalen Kurvensystems°) wird durch die Bewegung 
eines begleitenden Dreibeins beherrscht. Bezeichnet man mit t,, t,, t, die Einheitsvektoren 
der Kurventangenten, mit ds,, ds,, ds, die Linienelemente der Kurven, so wird sich ein 
begleitendes Dreibein t,. t,. t, beim Übergang von einer Ausgangslage in eine Nachbarlage, 
wobei sein Scheitel auf einer Kurve einer der drei Scharen fortrückt, in drei verschiedenen 
Weisen drehen. Die Änderungen eines mit dem Dreibein fest verbundenen Vektors r lassen 
sich also in folgende Form bringen: 


Or Or Or 


\ 
08, 29 


daber sind ıt,. u, die Vektoren der Drehung. 


Jeder der drei Drehvektoren läßt sich in Komponenten in Richtung der Tangenten der 
drei Kurven zerlegen. Man kann folgenden Ansatz machen: 


Die geometrische Bedeutung der Koeffizienten T und N ergibt sich in folgender Weise. Wenn 
man auch die Kurven der drei Kongruenzen nieht zu dreifach orthogonalen Flächen zusammen- 
fassen kann, so kann man doch jede Kurve in doppelter Weise mit einem Flächenstreifen 
verbinden, auf dem sie liegt, und dessen Normale in die Richtung einer der anderen Kurven 
fällt. Dann bedeutet N die Normalkrümmung, T die geodätische Torsion einer Kurve auf 
einem Flächenstreifen. Der untere Index bezeichnet die Tangentenriehtung der Kurve, der 
obere Index bei den N die Normalenrichtung des Flächenstreifens. 


= 5) M. Lagally: Dreifach-orthogonale Kurvenkongruenzen. Sitzungsberichte der Bayerischen Akademie der 
Wissenschaften (1932), S. 133. 
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Durch Einsetzen von (17) in (16) ergeben sich 9 Gleichungen, wenn man r durch t,,t,,t, 
ersetzt. Von ihnen sollen diejenigen 3 angeschrieben werden, welche die Anderungen von t, 
erkennen lassen; gleichzeitig sollen die Symbole für die Richtungsdifferentialquotienten 


durch die ausführlichere Schreibweise t,-V,t,-N,t,-V ersetzt werden. 


Aus diesen 3 Gleichungen kann man die Dyade Vt, ablesen: 


Ihre Invarianten sind: 

N’—T, 0 

N,'— N,?0 


Die Bedingung für die Existenz einer zur Kurvenkongruenz t, orthogonalen Flächen- 
schar ist bekanntlich 


t,- rott,—=0; oder T+T,=0; 
setzt man T, T, = T, so ergibt sich 


dabei bedeutet H die mittlere Krümmung, KA die Gaußsche Krümmung der Orthogonal- 
flächen. Von den Formeln (21) ist die erste sehr bekannt und auf anderem Weg leicht zu 
beweisen; auch die zweite findet sich schon (1912) bei Burali-Forti‘). 


6. Geometrische Invarianten des Geschwindigkeitsfeldes. Aus den Invarianten von Vt, 
(wofür jetzt wieder Yt geschrieben werden soll), solche von V d herzuleiten, ist in voller 
Allgemeinheit weitschweifig und mühsam; es soll deshalb nur auf den flächennormalen Fall 
eingegangen werden, in dem dv=yy oder mindestens proportional zu V y ist, wo y eine 


skalare Feldfunktion bedeutet. Setzt man |Y 9” M, so ist jedenfalls t M 2yg. Also ist 


(2), 


= 


Yt=M 


wenn man abkürzend die symmetrische Dyade P=|”ı setzt. 
Die erste Invariante von ® ist 


setzt man noch M- I so kann man die Invarianten von It, ins- 


besondere H und K, durch Simultaninvarianten von | g und J?y ausdrücken, die man auch als 
Differentialinvarianten von V y auffassen kann; einfache Zwischenrecehnungen sind im folgenden 
weggelassen. 


VYM?-LAg . (24e). 


‘), Burali-Forti: Fondamenti per la geometria differentiale eol metodo vettoriale generale. Rendieonti di 
Palermo 33, (1912). 
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Der Ausdruck für die Gaußsche Krümmung A der Niveauflächen läßt eine bemerkens- 
werte Umformung zu. Man überzeugt sich leieht durch Ausmultiplizieren, daß die reziproke 
Dvade von P = It; Yau durch folgende Determinante gegeben wird: 

5 


Setzt man die geränderte Determinante —- P, so ist 
P’ hat einen ähnlichen Aufbau wie 
Den Zusammenhang vermittelt die Cayley-Hamiltonsche Gleichung zwischen 4 aufein- 
anderfolgenden Potenzen einer Dyvade 
—- oder P®— PD, P'=0 .. (27), 
wo / die Einheitsdyade ist. Multipliziert man die letztere Gleichung zweimal skalar mit Vo, 
so ergibt sich die gesuchte Beziehung 


Damit ergibt sieh nach (24e) für das Krümmungsmaß der Niveauflächen der einfache Ausdruck, 


der übrigens auch durch direktes Ausrechnen zu bekommen ist‘). 
Bemerkenswert sind die Vereinfachungen von HM und K im Fall der Laplaceschen 


Felder; für ist: 
3 


Der letztere Ausdruck läßt eine weitere Vereinfachung zu. Berechnet man 

... Bl), 
so sieht man, daß für IQ -0 auf der rechten Seite die dritten Differentialquotienten hinaus- 
fallen. Für 19 =0 wird überdies auch PD: P= —2®P,; also ergibt sich für die Gaußsche 
Krümmung der Niveauflächen eines Laplaceschen Feldes 


| 


7. Zusammenhang mit dem Hamelschen Satz. Die letzten Umformungen geben, wenigstens 
für den Fall der Laplaceschen Felder, einen interessanten Zugang zum Hamelschen Satz. 
Bemerkt man, daß M = ist, so wird 


| | 


bei Wiedereinführung der alten Bezeichnung; nimmt man hierzu den Energiesatz (Bernoulli- 
sche Gleichung) als physikalische Tatsache: 


») 


r 
+/ +V = cosst, 
so geht aus (35) der Hamelseche Satz in seiner speziellen Fassung (12) hervor: 
686 


') Paseal-Timerding: Repertorium der Höheren Mathematik, II, 2. Kap. XLI. — Salkowski: Allgemeine 
Flächentheorie, S. 1074. 
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Über die Tragfähigkeit eines längsbelasteten Plattenstreifens 
nach Überschreiten der Beullast. 


Von K. Margwerre, Berlin-Adlershof, und E. Trefftz 7, Dresden '). 


I. Einleitung. Die folgende Untersuchung beschäftigt sich mit der Frage nach der Trag- 
fähigkeit eines Plattenstreifens, der in der Längsrichtung über die Knickgrenze, bei welcher 
der Streifen ausbeult, auf Druck beansprucht ist. Als Näherung für den Streifen, dessen 
Länge groß ist gegen seine Breite, wird der unendlich lange Streifen behandelt, der an seinen 
Längsseiten an starren Trägern drehbar, d. h. momentenfrei befestigt sei. 

Wir verwenden die energetische Methode. Den Ansatz für die Formänderungsenergie 
eewinnen wir, da es sich um eine dünne Platte handelt, in derselben Weise wie in der 
klassischen Biegungstheorie, die wir auf den Fall erweitern, daß die Durchbiegungen nicht mehr 
klein sind gegen die Plattendicke h. 

Koordinaten und Bezeichnungen. Wir legen die .«-Achse in die Plattenmitte, die y-Achse 
quer dazu, die z-Achse senkrecht zur Plattenebene. Die Plattenbreite sei b, die Plattendicke h; 
E. G und m seien die Elastizitätszahlen des Materials, 

Eh’ m? m 


N 12. (m? -1) 


(der sog. Plattenfaktor. 
Il. Differentialgeometrie der Plattenbiegung. Ein Punkt der Plattenmittelebene habe vor 
der Deformation die Koordinaten x und 9, bzw. den Ortsvektor 


wo &.€, €, Einheitsvektoren in den Achsenriehtungen sind. Bei der Deformation erfahren 
die Punkte der Mittelebene die Verschiebung 


und bilden nach der Deformation ‚die Fläche . 


Die wesentliche Voraussetzung der Plattenbiegungstheorie ist, daß die Normalen zur Mittel- 
ebene nach der Deformation geradlinig und normal zur deformierten Mittelfläche bleiben. 
Ist also W ein Einheitsvektor in der Normalen zur deformierten Mittelfläche, so geht ein 
beliebiger Punkt 


bei der Deformation in den Punkt 
über. Um den Normalenvektor W% zu bekommen, bilden wir die Vektoren 
0x 


die tangential zu den Kurven liegen, welche bei der Deformation der Mittelebene aus den 
Koordinatenlinien entstehen. Ihr vektorielles Produkt r.xX x, hat die Riehtung der Normalen. 
Diese selbst wird also 


Ir I, Ur x 


Die für die Energie der Biegung maßgebenden Krümmungsgrößen sind dureh die Ableitungen 
des Normalenvektors bestimmt. Um 


ON 
(y . . . . . . . . . . (‘ ) 


\ 
(x 


zu berechnen, beachten wir, daß x. und r, auf N senkrecht stehen: 


!, Fürdie Fassung der vorliegenden Arbeit zeichnet der oben zuerst genannte Verfasser allein verantwortlich; 
er hat sich bemüht, eine Darstellung zu geben, wie sie im Sinne seines verehrten Lehrers gelegen hat. 


| 

In- 
9), | 
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he 
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Differentiation ergibt 


bzw. (10). 
Kerner folet aus 1 


Mit Hilfe der Gl. (10) und (11) können X, und NW, durch r, und r, dargestellt werden. 
Setzen wir 

und multiplizieren diese Gleichung skalar mit Wr, und r,, so folgt aus (9) und (11) 


Für a und 5 ergeben sieh die linearen Gleichungen 


Daraus folgt mit der Abkürzung 


Wir erhalten also 


und entsprechend 
Ill. Die Formänderungsenergie. Die Komponenten des Verzerrungstensors, durch die wir 
die Formänderungsenergie (im folgenden abgekürzt: F. E.) ausdrücken werden, erhalten wir 


dureh die Betrachtung des Linienelementes. Für zwei Punkte mit den Koordinatendifferenzen 
d.r,dy.dz ist das (Juadrat ihrer Entfernung vor der Deformation 


nach der Deformation 
. . . . (0) 


bzw. unter Berücksichtigung der Gl. (9 und 


wo zur Abkürzung 
wesetzt Ist. 
Die Änderungen in den Koeffizienten des Linienelementes im deformierten Zustande 
werenüber «den Werten des undeformierten Zustandes 


\ L „2m 2 


. 2 20) 2 

Ya +221, 


sind die nieht-linearen Verzerrungskomponenten der elastischen Deformation. Ihren Zusammen- 
hang mit den Verzerrungsgrößen der gewöhnlichen Elastizitätstheorie erkennen wir, wenn wir 
(die Plattenmittelfläche betrachten, d.h. 2-0 setzen, und für r., und rt, ihre aus Gl. (3) 
folxzenden Werte einsetzen: wir erhalten dann 
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An den linearen Anteilen sehen wir, daß y,, bzw. y,, dem Doppelten der Dehnung e, bzw. &, 


! entspricht, und y,, der Winkeländerung yx,„. Die F. E. für die Volumeinheit beim zwei- 
. 
dimensionalen Spannungszustand 
0) m Y \ 
wird also, wenn wir das Hookesche Gesetz sinngemäß auf unseren Fall großer Deformationen 
übertragen ?), 
). 
(i m \ 
2) Setzen wir hier die Werte der y aus (23) ein und integrieren über die Plattendicke, d.h. 
von z=-— h/2 bis z= +h/2, so erhalten wir die F. E. für die Flächeneinheit der Platten- 
| mittelfläche 
ar 2 - ar 2 A 
Ye) dz eure Yo)dz . . . . . . . . (Zu) 
hi? — hi? 
4). 
| G m | 
2 2_9]? be 
4m 1 | + + 
>), 


)), 
| (28). 


S0 
‘) 
S). | 
Die hier vorkommenden Größen finden wir in der folgenden Weise. Aus (3) folgt: 
ip 
also 
und daraus 
)) 
(1). 
), 
Nach (7) ist 
(32); 
] 
| ferner nach (17/18): 
| \ 14% 2 
' worin die n,„ sich nach (10) ergeben zu: 
r 
?) Vgl. E. Trefftz: Über die Ableitung der Stabilitätskriterien...... Verhandl. des IIl. intern. Kongresses 
f. techn. Mech., Stockholm 1930. Bd. IIL, S. 44. 
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Vereinfachung. Bis zu dieser Stelle haben wir außer der Voraussetzung der Platten- 
biezungstheorie, daß die Normalen der Mittelebene bei der Deformation geradlinig und senkrecht 
zur deformierten Mittelfläche bleiben, keine Vereinfachungen vorgenommen. Um den Aus- 
druck für die F. E. jetzt in eine rechnerisch brauchbare Form zu bringen, entwickeln wir 
ihn nach Potenzen der Differentialquotienten von u,v,w. Durch eine Größenordnungs- 
betrachtung vergewissern wir uns, welche Glieder beibehalten werden müssen und welche 
zestrichen werden können. Nun weiß man, daß im Augenblick des Ausbeulens unter reiner 
Längskraft der gelenkig gelagerte Plattenstreifen von der Breite b die Form annimmt: 


m f cos co 
/ h 


(39), 
worin «er Biegepfeil f unbestimmt bleibt. Machen wir die Annahme, daß diese Beulform 
bei nieht zu großer Überschreitung der Beullast erhalten bleibt, so ist also w von der Größen- 
ordnung f, und und von der Größenordnung f/b. Aus Symmetriegründen (d.h. weil sie 
vom Vorzeichen von f unabhängig sein müssen) sind die Verschiebungen «, » mit f? propor- 
tional. ihre ersten Ableitungen also proportional (f/b)’: es sind daher z. B. «, und (w,)?’ von 
der gleichen Größenordnung, während (w,)’ gegen diese Größen klein ist. 


Bei der Integration der F. E. pro Volumeneinheit (Gl. 25) über die Plattendicke h treten 
Glieder auf, die mit h, proportional sind: 


a, hy, tur, + h’y,(u,r, ww). 


Nun ist A/b ebenfalls eine kleine Größe: und für das Verhalten einer Platte oberhalb der 
Beullast ist es kennzeichnend, daß f und k von der gleichen Größenordnung sind. Wenn 
man daher in y, z. B. die Glieder » +2-ter Ordnung in f/b mitnimmt, so muß man in y, die 


Glieder »-ter, in y, die »n —2-ter Ordnung beibehalten. — Die Zusammendrückung im Augen- 
blick des Ausbeulens ist bekanntlich proportional mit (h/b)’, ist also von.der gleichen 


Größenordnung wie 


Behält man in (25) nur die unter diesen Gesichtspunkten niedrigsten Glieder bei (h f* und 
so erhält man: 


Im 


(36). 


Auf dieser Stufe der Näherung setzt sich also die gesamte F. E. aus Dehnungsenergie da, und 
Biegungsenergie «A, additiv zusammen, was zuweilen als selbstverständlich betrachtet und 
unmittelbar zum Ausgangspunkte der Rechnung gemacht wird. 


IV. Der gedrückte Streifen. Wird der Plattenstreifen in der Längsriehtung gedrückt, so 
besteht vor dem Ausbeulen zwischen der gleiechförmig verteilten Spannung o, und der Dehnung # 
(z. B. bei 0, VO) der Zusammenhang: 0, Er, wofür wir, wenn wir mit p den mittleren Druck 


bi2 
— 
und mit &,: & die mittlere Stauchung bezeichnen, auch schreiben können 


Oberhalb des kritischen Wertes der Zusammendrückung &* beult der Streifen aus, und zwar 
in der Mitte stärker als in der Nähe der festgehaltenen Ränder; infolgedessen verteilt sich 
(der Längsdruck ungleichförmig: die mittleren Fasern tragen, wie man zu sagen pflegt, „nicht 
mehr voll mit“. 


Gesucht ist der Formänderungszustand im ausgebeulten Blech und insbesondere die neue 
Spannungsdehnungs-Kurve, d. h. der Zusammenhang zwischen mittlerem Druck und mittlerer 
Stauchung: 


Die Zusammendrückung #, (die für die gerade bleibenden Randprofile auch oberhalb der 


Beullast „Stauchung* im eigentlichen Sinne des Wortes ist) betrachten wir als die unabhängige 
Veränderliehe unseres Problems. 


| 


sie 


on 


en 
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Es ist zweekmäßig, die Verrückung u (#y) eines beliebigen Punktes der Plattenmittel- 
täche zu unterteilen in eine „mittlere“ Verschiebung — #, «= und eine zusätzliche Verschiebung 
y), also zu setzen: 

Amlgist (40), 


wenn wir unter + ®,5 die seitlichen Verschiebungen der (geradebleibenden) Längsträger 
9 


verstehen 
Für die Zusatzverschiebung # haben wir dann als Randbedingungen an den Rändern 
+ b]2 in beiden Fällen 


Für die Verschiebungen # fordern wir gleichfalls 


d.h. wir setzen voraus, daß die Längsträger hinreichend dehnungssteif sind, um durch die 
von der Platte ausgeübten Scherkräfte nicht deformiert zu werden. Diese Annahme ist übrigens, 
wie die Rechnung zeigt, auf das Ergebnis ohne merklichen Einfluß. — Entsprechende Rand- 
bedingungen an der Schmalseite einer Platte von endlicher Länge werden beim unendlich 
langen Streifen durch die Forderung ersetzt, daß “ und # endlich bleiben, was für unseren 
Fall darauf hinauskommt, daß und in periodisch werden. (# an den Knotenlinien 0.) 

Die Gleichgewichtsgestalt eines elastischen Systems bestimmt sich nach dem Prinzip 
der virtuellen Verrückungen aus der Forderung, daß die potentielle Energie (d. h. die Differenz 
von F. E. und Arbeit der äußeren Kräfte) ein Minimum wird. Da wir den Gleichgewichts- 
zustand bei vorgegebener Zusammendrückung, d. h. vorgegebenen Verschiebungen der Ränder 
untersuchen, so ist die Anwendung des Prinzips besonders bequem; denn indem wir bei der 
Variation des Formänderungszustandes die Ränder festhalten a = leisten die äußeren 
Kräfte?) bei der virtuellen Verrückung keine Arbeit, und die Minimalforderung für die poten- 
tielle Energie reduziert sich auf die Minimalforderung für die F. E. selbst. Die gesamte 
Formänderungsenergie A der deformierten Platte erhalten wir, indem wir den Ausdruck (36) 
für die F. E. der Flächeneinheit über die Platte integrieren. Mit den Ansätzen (40) für u 
und » wird | 


A- a,dedy=Gh = \ 28, +2, 
Gh? 1> 1). \ dady— (Wr 


Von der Verschiebung »w senkrecht zur Plattenmittelebene wissen wir, daß sie beim Erreichen 
der Beulgrenze in der Form (35) beginnt. Wir machen nun die Näherungsannahme (eine 
Erweiterung wird in Abschnitt IX durchgeführt werden), daß auch nach Überschreiten der 
Beullast nur die Amplitude f sich mit wachsender Last verändert, d. h. daß die Beulen sich 
ähnlich vergrößern. 

Dann ist A nur noch Funktion der Verschiebungen (x, y) und f(x, y) und der 
von x und % unabhängigen Größen f,*, und e,. Die Forderung de = Min liefert zunächst für 
die Verschiebungen # und # zwei simultane lineare Differentialgleichungen, mit deren Hilfe 
sie sich als Funktionen von f,*,unde, bestimmen lassen. Aus 


ergibt sich dann f als Funktion von x, und e,. 
Sind die (starren) Längsträger derart befestigt, daß sie als gerade Stäbe seitlich aus- 
weichen können, so daß sie im Mittel keine Querspannungen o, aufnehmen, dann ist 


oA 


In dem (später zu behandelnden) entgegengesetzten Extremfall vollkommen unverschieblicher 
Längsständer ist 


3) Quer zur Platte,(d. h. auf Ober- und Unterseite) sollen’ja keine Kräfte wirken. 
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KEliminiert man #, mit Hilfe von (45) oder (46), so ist der gesamte Formänderungszustand als 
Funktion von bestimmt. 


V. Berechnung der Verschiebungen. Aus der Forderung, daß der Ausdruck (43) für alle 
am Rande verschwindenden d # und 9% einen Minimalwert annehmen soll: 


erhalten wir mit der üblichen Schlußweise der Variationsrechnung für “# und # die partiellen 
Differentialgleichungen 


(da, /da, (da, /da, 
bzw. nach elementarer Umrechnung 
(45), 


mit den oben abgeleiteten Randbedingungen (41), (42): 
für y=+bl2. 


Dieselben Gleichungen erhalten wir, wenn wir das Gleichgewicht einer Platte bestimmen, die 
in ihrer Ebene von den Kräften 


Y xWyy yWxy 
COS sin 
—] h 
(49) 
— (rh wm — W, 
3 ») ») 
COS l|sın 
m—i 


beansprucht wird und den Randbedingungen (41), (42) genügt. Als Partikularlösung der voll- 
ständigen (unverkürzten) GIn. (48) finden wir 


- 


7 
165 


2ny m—! nf? . 2ny[| m—! 

Hierzu tritt zur Erfüllung der Randbedingungen (41), (42) eine Lösung der verkürzten Gleichungen 


(ie wır mit dem Ansatz 


gewinnen. 
Setzen wir zur Abkürzung 
) 

h 


so ergeben sieh für pop) und die Differentialgleichungen 


Die Lösungen dieser Gleichungen lauten: 


(DD). 


a+Bdm— Sniy— 2m 


e 
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wo a und ß Integrationskonstanten sind. Insgesamt erhalten wir also 


sinix +cos/iy+ + (3m - +2 (Im + DAySinäy!, | 
165 \ m J = 
sinäy + HP) SinAy 


Die Integrationskonstanten «a und 5 bestimmen sich aus den Randbedingungen (41), (42) zu 


Im Grenzfall m =» wird a0 und $m=0. Für den mittleren Wert m = 10/3 wird 
und damit 
+ 00849 0.0136 y+ 001264 Sind (9), 


Zur numerischen Durchführung der Rechnung schreiben wir die F. E. der Flächen- 
einheit in der Form . 


h 
2 m - 4 1 2 + x + 
0a, Od, oda, Od, , 


Integrieren wir, um die gesamte F. E. zu erhalten, so fallen die Integrale 


\ 


\ \( x I \ \ F x 73 dy 
0 —b]2 —bi2 


fort, wie man erkennt, wenn man in Gl. (47), die für beliebige, den Randbedingungen genügende 
Funktionen und gilt, speziell und setzt. Ferner wird 


h 


was man hier durch Einsetzen von w — fcos „ 608 2 feststellen kann; (63) gilt bekanntlich 
überhaupt für jede geradlinig begrenzte Platte mit festgehaltenen Rändern. Fine weitere Ver- 
einfachung erzielen wir, wenn wir beachten, daß die Glieder, welehe nur die Differential- 
uotienten von und enthalten: 


b b2 


.. . . r . A . 
die Formänderungsenergie darstellen, die bei den Verschiebungen # und ö in der Platte ent- 
halten ist. Nun sind aber, wie wir oben bemerkten, # und # die Verschiebungen, welche sich 


N 

| 
n | 

| 
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einstellen, wenn die am Rande festgehaltene Platte in ihrer Ebene pro Flächeneinheit von 
den Kräften (siehe Gl. 49) 
Im 
X: Nina 
beansprucht wird. Die in einem elastischen System gespeicherte Energie ist aber gleich der 
Arbeit, welehe die äußeren Kräfte bei allmählichem Aufbringen leisten. Das Integral (64) 


wird also 


cos/y-+| 


b 

0 — bR 


Wir erhalten schließlich die in einem Streifenstück der Periodenlänge 5b enthaltene F. E. aus 
den vier Integralen 


Gm h 
M 
h hir 


x (i hı ) 


(67) 
b 
h b]2 
Gm 
J, \ IP 
b]2 


in der Form J, +J,—- J,+J,: 
Setzt man hier die Werte von mw aus (35), von « und ® aus Gl. (56) und von X und Y 


aus (6») ein, so erhält man 


J, (i h 4 


(68), 
m 1 


umel damit 


m 
64 


Hier ist in dem unwichtigen Gliede VO, 5; welches von den Integrationskonstanten «a 
und > herrührt, m = 10/3 gesetzt. 
Aus der Minimalbedingung YA 0 und der Bedingung ge —0 (Gl. 45) erhalten wir 


(die Gleichungen: 
m m-— _ 
m-+130? m-+-1 sb? | 0.01) 
und 
| 


m h? m 
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also 


7.) 


Der von f unabhängige Anteil von &;: 


2 2 
m 
2 2 x . . . . . . . . . . . 
30 m?’ —1 
ist die aus der elementaren Theorie bekannte kritische Dehnung, bei welcher das Ausbeulen 
beginnt. Es wird also 


2 


VI. Berechnung der Spannungen im ausgebeulten Streifen. Überträgt man das Hookesche 
Gesetz auf die Elastizitätstheorie endlicher Deformationen, so lauten die Spannungs-Dehnungs- 
sleichungen ebenso wie in der klassischen Theorie (für den hier vorliegenden ebenen Spannungs- 
zustand): 
22 Y 

(i 


Die Spannungen im Blech setzen sich aus den Dehnungsspannungen und den Biegungs- 
spannungen zusammen. Die Dehnungsspannungen erhalten wir, wenn wir in Y,1, Yass Ya Ihre 
Werte für 2=0 (Gl. 24) einsetzen: 


+2ü, +20, tt, tn. . (76). 


Bei der Ausrechnung lassen wir in den Formeln für « und #® jetzt die kleinen, von den 
Integrationskonstanten « und 5 herrührenden Glieder weg und schreiben einfach 


| .„ af... m-1 
165 m 165 m 
Damit erhalten wir (2°-0)) 
n?f? . Im 1 . (9) 
‘ m 4b’ 


.. 
2 


cosiy 
m 1 m 


bzw., wenn wir für e, und &, ihre Werte aus (74) und (71) einsetzen und statt o, den Längs- 
druck p = — 0, einführen: 


Mer 


Wir haben also das einfache Ergebnis, daß die Schubspannungen verschwinden und die 
Normalspannungen in ihrer Wirkungsrichtung konstant sind. Der Mittelwert der Zug- 
spannungen o,„ über eine Periode ist null, in Übereinstimmung mit der Voraussetzung, daß 


die Träger an den Längsseiten frei verschieblich sein sollen. Den Längsdruck p = — 9, zer- 
legen wir in einen mittleren Druck 


m +1 


m 


2 


(79), 


cos/x 


cos 


und einen mit y veränderlichen Anteil 
7? f? 


| 
und daraus die Dehnungsspannungen 
| 

| 
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Bemerkenswert ist, daß die Amplitude dieses veränderlichen Anteils ebenso groß ist wie die 
Erhöhung des mittleren Druckes über den kritischen Druck 

Em’ a®h? 

— 

m’ — (83) 
In der Mitte, wo das Blech am leichtesten nachgeben kann, steigt der Druck also nicht über 
den kritischen Wert, während die den Rändern benachbarten Teile das Doppelte der mittleren 
Druckerhöhung aufnehmen. 


Definieren wir die mittragende Breite b,, des Bleches durch die Gleichung 


D 


| 


b| 


Pınax ‚pr 


wo nach Gl. (80): 


Es=p*+E7T 


der am Rande herrschende maximale Druck ist, so ergibt sich wegen p = 3 (Pax + pP”): 


)= 


I max 


145]. 


> 


Die mittragende Breite nimmt also mit wachsender Belastung ab. Würde unsere Näherungs- 
theorie bis zu beliebig hohen Lasten gelten, so würde in der Grenze #*/e, >0 gerade der halbe 
Streifen mittragen. 


Aus den Gl. (74), (SL), (83) erhalten wir 


d.h. der Zusammenhang zwischen Spannung und Dehnung ist auch oberhalb der Knicklast 
linear, mit einer „scheinbaren Dehnungssteifigkeit“, die unabhängig von m, halb so groß 
ist wie im eigentlich-elastischen Bereich unterhalb der kritischen Last. 


Vil. Der Fall festgehaltener Längsträger. Zum Vergleich sei noch das Ergebnis für den 
"all angegeben, daß die Längsträger von Anfang an unverschieblich festgehalten werden. 
(Geht man bei der Untersuchung dieses Falles wieder aus von der Beulform (35), so bleibt 
der Reehnungsgang bis Gl. (69) genau der gleiche; an die Stelle der Gl. (72) und (81) aber 
treten dann (70) und (79) mit &,—=V, also: 


m = m im? (88) 

Nun sind aber infolge des aus der verhinderten Querdehnung entstehenden Querdruckes die 
Beulenfelder im Augenblick des ersten Ausbeulens nicht mehr quadratisch. Man muß an 
Stelle von (35) vielmehr ansetzen: 

Wiederholt man, was hier nicht in extenso vorgeführt werden soll, mit diesem erweiterten 
Ansatz die bisherige Rechnung, so erhält man an Stelle von (69): 


f' N 41 (0), 
und statt (58): 
(1). 


- 


f?/b® | 
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Das in diesen beiden Gleichungen noch unbestimmte Verhältnis b/! bekommt man aus der 
‘übrigens aus der elementaren Theorie bekannten’) Beziehung: 


h? 


| 


- 


lie im Augenbliek des Ausbeulens (f=0) gilt, vermöge der Minimalforderung \bin 0: es 
( 


ergibt sich: 


N Mm — 2 (92) 


I ist also von der @Querkontraktionszahl sehr stark abhängig; im Grenzfall der Inkom- 
pressibilität (m = 2) würde die Wellenlänge in der Längsriehtung sogar über alle Grenzen 
ochen aber auch für mittlere m-Werte unterteilt sich der Streifen in ausgesprochen ungleich- 
seitige Rechtecke (b/1=0,4 für m 10/3). Geht man mit (92) in (OL) ein, so folgt: 


m ’+2 ,„, ”f3m—1+2im? 
& 272 3 =e"+5 (93), 
m sh m sh im 
Der kritische Wert liegt bei | 
(m — 10/3), 


also ungleich tiefer als im Falle verschieblieher Ränder. Dafür sinkt die scheinbare Dehnungs- 
steifigkeit oberhalb &** sehr viel weniger ab. Eliminiert man f aus (93) und (94), so findet 
man: 


Im+ 
m” 


m? — 1 > m? - 13m 


3m —1+ 


m“ 


Die scheinbare Dehnungssteifigkeit verhält sieh also zur elastischen Dehnungssteifigkeit wie 
2/3 Va); im Gegensatz zu dem Fall der beweglichen Ränder ändert sich dieses Verhältnis 
mit u, und zwar zwischen den Grenzen 0.52 für m 2 und 0.67 für m. Für die prak- 
tisch in Frage kommenden m-Werte liegt es etwa bei 0.75, also sehr viel höher als im 
anderen Fall. 

Es mag an dieser Stelle erwähnt werden, daß man eine spürbare Abweichung von den 
„im ersten Augenblick“ gültigen Geraden (87) und (95) (den Tangenten an die „wirklichen* 
Spannungs-Dehnungskurven oberhalb des kritischen Punktes (a) und te) ın Abb. 2 —) er- 
halten kann, wenn man das Verhältnis !/b im Ansatz (89) in Abhängigkeit von e, aus der 
Forderung < 10° O nach (90) bestimmt, anstatt, wie dies oben geschah, seinen von &, unab- 
hängigen Wert am kritischen Punkte einzuführen. — Die Erörterung dieser Verbesserungs- 
möglichkeit soll aber (zusammen mit anderen Erweiterungen) einer späteren Veröffentlichung 
vorbehalten bleiben. 


VI. Die Biegungsspannungen. Für die Biegungsspannungen erhalten wir mit ww: 
[eos ,, cos , (der Fall verschieblicher Ränder!) wie in der klassischen Biegungstheorie 


m’—1\ 79% m--l 2 2 
. m+1 b? 2 2 


°, 7. B.: Handbuch der Physik VI, S. 204 (1928). 
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Von Interesse ist das Verhältnis des Höchstwertes 
Em a’hf 


Max m h? (97 
zum Höchstwert des Druckes (Gl. (85): 
Em’ n’h” En’f? 
Wir erhalten 
Für f/h ergibt sich aus (SI) und (83) 
__ 8m? 
(4 p* (200). | 
Setzen wir m 10/3, so wird 
h p 
und «damit 
max _ (102). 


Man sieht, daß bei einiger Überschreitung der Beulgrenze die Biegungsspannungen von gleicher 
Größenordnung sind wie die Druckspannungen in der Längsrichtung. 
nz 

IX. Erweiterung des Ansatzes ı" = f cos „es, Die bisherigen — sehr einfachen und 
übersichtlichen — Ergebnisse wurden gewonnen unter der Voraussetzung, daß die Beulform (35) 
für w auch bis zu großen Amplituden f der Ausbeulung hin erhalten bleibt. Experimentelle 
Beobachtungen deuten aber darauf hin, daß schon kurz nach Überschreitung der Beullast 
sich die Beulen in der Weise verformen, daß die Querprofilform nach der Mitte zu sich ab- 
lacht. Diese Abweichung von der reinen Sinusform läßt sich qualitativ am bequemsten er- 
fassen, wenn man den Ansatz (35) erweitert in der Form: 


37 


Die Erweiterung durch ein Glied f, cos”, COS ni bietet sich in natürlicher Weise an, da 
(diese Funktion die nächste, den Randbedingungen # = m,,„=V für y == + b/2 genügende gerade 
Funktion ist. Das Minuszeichen wird zweckmäßig in den Einsatz hineingenommen, um die 
zu erwartende Abflachung in der Mitte (y 0) zum Ausdruck zu bringen. 

Die weitere Rechnung verläuft ganz wie bisher: Nach Bestimmung der # und # und 
Auswertung der Integrale wird die F. E. Funktion von #,,e,f, und f,: die Minimalforderungen 


04 0A 
.0 und die Bedineune 


erlauben die Größen f,. f, und e, und damit den Spannungs- und Formänderungszustand als 
Funktion von #, zu bestimmen, 

Die Rechnung läßt sich sehr vereinfachen, wenn wir, indem wir uns auf den Fall ver- 
schieblieber Längsränder beschränken, die beiden (durch die bisherigen Ergebnisse gerecht- 
fertigten ’) Hypothesen einführen: 


- 


2. r(y=+b2) 0, d.h. Gleiten der Platte an den Längsversteifungen (aus dieser 
Annahme ergibt sich im Gegensatz zu der anderen Grenzannahme (42) [r =0]- 
eine untere Schranke für die scheinbare Dehnungssteifigkeit der Platte, was für 
den vorliegenden Zweck einer Abschätzung des Einflusses des Zusatzgliedes von 


(105). 


) Man kann übrigens beweisen, daß im Falle —=0 und (x, =" die Querkontraktionsziffer herans- 


fallen muß. 


— 
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Der Ausdruck für die F. E. der Volumeinheit nimmt unter der Voraussetzung (105) die 


einfache Form an: 


worin 
Ya 2, +20, — Ya — 22m yy, 


Für die in den Verschiebungen quadratischen Anteile der y finden wir nach (103): 


f, cos 2 f,f, «cos 608 5 f,” cos 


und entsprechend 


/"sinzy  2f, f, 


Die Minimumforderung (#7) für « und # führt auf die Differentialgleichungen: 


mt 
und 
Y 
71° 


Partikularintegrale der Gl. (110) erhält man in der Form: 
\ 


7 


sind 
165 


165 


(107), 


(108). 


(109). 


(110) 


(112). 


(113). 


Man erkennt, daß diese Lösungen für y==+b/2 den Randbedingungen #=0 und 7-0 


la h. wegen # = 10.0, einfach 0) genügen. Wir können also unmittelbar anschreiben: 


| 


2 
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Die Schubspannung r (Gy, verschwindet diesmal also nicht, und entsprechend sind 
die Normalspannungen @y, und oo!) in ihrer Wirkungsrichtung ver- 
änderlich. Bei der Ausrechnung der F. E. zeigt sich allerdings, daß diese Anteile nur den 
Zahlenkoeftizienten von fj’. fs um etwa 10° beeinflussen, so daß die von anderen Autoren 
eingeführte vereinfachende Annahme 7-0 keine wesentlichen Fehler im Gefolge hat*). 


Die gesamte F. E. setzt sich, wie früher, aus der F. E. der Dehnung: 


b 


und der F. E. der Biegung (m = x!) 
b bj2 
— HR 
zusammen. Da Integrale über Produkte von Kreisfunktionen verschiedener Frequenz im 
Intervall ©... verschwinden, vereinfacht sich der Ausdruck für A sehr stark: 


\ 


| (117). 


S 


Der umrandete von den veränderlichen Anteilen der Spannungen herrührende Betrag 


| = ist in der Tat gegen 33- nicht groß. Läßt man zur Vereinfachung 
der Schreibweise das unerhebliche Glied 1/25 beiseite, so gilt: 


- 


> 


| (118). 
Für ergibt sieh: 
| 
Die Minimalforderungen (104) liefern: 
(120), 
| 1" 


- 


Setzen wir (121) in (120) ein und lösen nach #, auf, indem wir wieder 32 durch #* abkürzen, 
so ergibt sich: 


(122). 


") Vgl. die in der Schlußbemerkung angekündigte Arbeit. 


um | 

= | 


; 
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Nach (114) ist der mittlere Druck durch den von y und ‚= unabhängigen Anteil von 


E _ 
5 gegeben in der Form: 


T 
p= E €, — ) . . . . - - (123). 


An Stelle des Gleiehungssystems (122) können wir daher auch schreiben: 


(124). 


25 + 16 Is +82f; 


ll 


Durch (122) bzw. (124) sind die Amplituden f (d. h. Größe und Form der Beulen) als Funk- 
tionen der Zusammendrückung bzw. der Druckkraft grundsätzlich bekannt. 


Da es nieht möglich ist, die Systeme nach f, und f, einfach „aufzulösen“, soll ihre Ab- 
hängigkeit von &, und e* bzw. von p und p® in anderer Weise diskutiert werden. 
Zunächst bemerkt man, daß f, 0 nieht Lösung eines dieser Gleichungssysteme ist. 
Die neue Beulform (f, 5£0) entsteht daher nicht erst oberhalb einer neuen Instabilitätsstelle, 
sondern geht aus der Form (35) stetig hervor. 
Führt man in (122) an Stelle von f, den dimensionslosen Parameter 
ein, so läßt sich durch Elimination von f, dieser eine Parameter als Funktion von »,[e* ge- 
winnen. Man findet: 
45 


—e, 19. — 840° 


(126); 


(diese Gleichung erlaubt in der Tat z,/e* als Funktion von £ bequem zu bestimmen. In Abb. I 


Abb. 1. 


ist das Ergebnis: in Abhängigkeit von »,/s* aufgetragen. Man erkennt, daß in einem: 
e,/£ Bereich, in dem die entwickelte Theorie hier überhaupt gilt, sehr klein bleibt. Man 
kann übrigens auch leicht zeigen, daß die Kurve bei © = 0.0665 eine Asymptote hat, d.h. daß 
der Wert © = 1/15 (aueh für noch so hohe Überschreitung der Knieklast) nicht überschritten 
werden kann. 


Die Wellenform ist nach (103) gegeben durch: 


w(e,y)=f, cos c0o8 —- — [cos 


Die Abtlachung der Profile == const, ist also ganz unbedeutend’). 
Aus (122) und (124) ergibt sich nunmehr: 

3 


- 


21—-34C 2 
Da für #, = verschwindet, beginnt die neue Spannungs -Dehnungskurve wieder mit der 
Fangente E]2, und auch im weiteren Verlauf weicht sie selbst bei sehr hoher Überschreitung 
der Beullast von der Geraden (87) nur um wenige Prozent ab. 


| ‘) Eine „siehtbare“ Verformung, d.h. das Auftreten einer Gegenkrümmmung, ist, wie man analytisch leieht ein- 
Sieht, nur bei <> 1/9 möglich, kann also in unserem Falle nieht eintreten. 


| 
| 
| | 
| 904} 
| 
0 Ai 5 70 75 20 25 
j 
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Aus der ersten Gl. (122) läßt sich noch der Biegepfeil f, dureh &, und den Parameter £ 
ausdrücken: 
_ 


Die Maximalspannung p,,., an den händern ist natürlich (vgl. (114)) wieder gegeben 
dureh: 


Der Minimalwert in der Mitte aber steigt mit wachsender Zusammendrückung &, über die 
kritische Spannung p” hinaus: 
für y=V daher: 


2.) 


wofür sich wegen (131) auch schreiben läßt: 


Pin =P* + El - Ele, . (135). 


Für die mittragende Breite b,, ergibt sich: 


Die mittragende Breite wird daher kleiner als nach (86), 
praktisch äußert sich das aber nur darin, daß der 
Wert 1/2 etwas früher erreicht und dann sehr wenig 
unterschritten wird. 


Die Biegespannungen lassen sich nach (96) leicht 
errechnen. Die Größtwerte in der Mitte eines Feldes 
werden um ein weniges abgemindert; bezügl. 6, treten 
näher den Rändern noch einmal Extremwerte auf, die 
aber auch nieht höher liegen als das frühere Maximum 
Gh | | nach (97). Die höchsten Biegespannungen sind daher 
wieder von der gleiehen Größenordnung wie die höch- 
sten Druckspannungen. 


Ä In Abb. 2 sind die Spannungs-Dehnungskurven, 
/ wie sie sich nach (57), (95) und (125) ergeben, zu- 
sammen mit Versuchsergebnissen von €, „. (Kurve d), auf- 
ke ER N getragen. Man sieht, daß die Verbesserung, die man 
dureh die Erweiterung (105) erreicht, außerordentlich 
nr | geringfügig ist. (Vgl. aber die Ergebnisse der in dem 

unten genannten Aufsatz durchgeführten Rechnung.) 


X. Zusammenfassung, Bemerkung. Die Tragfähigkeit dünnwandiger versteifter Bleche 
ist bekanntlich beim Erreichen der Blech-Beul-Last keineswegs erschöpft. Eine theoretische 
Untersuchung des Verhaltens ausgebeulter Bleche kann mit Hilfe der gewöhnlichen (lineari- 
sierten) Elastizitätstheorie nicht geleistet werden. Für die praktische Anwendung der Elasti- 
zitätstheorie „endlicher” Deformationen bieten sieh in erster Linie die energetischen Näherungs- 
methoden (Ritz) an. Wegen der Unsicherheit, die diesen Methoden bei mehrdimensionalen 
Problemen anhaften können, ist es von Wichtigkeit, mit einem Minimum an willkürlichen 
Voraussetzungen auszukommen. 

Als ein Beispiel für eine in diesem Sinne möglichst strenge Rechnung wird der Spannungs: 
und Formänderungszustand im gelenkig gelagerten Plattenstreifen unter reinem Längsdruck 
berechnet. 

Kine etwas andere Darstellung der Ergebnisse dieser Arbeit sowie ein Vergleich mit 
den Arbeiten anderer Autoren erscheint im Märzheft der Zeitschrift Luftfahrtforschung (1937). 
Dort wird auch über eine andersartige Erweiterung des Ansatzes (35) für die Ausbeulform 
und über eine ‚Vereinfachung des Rechnungsganges berichtet. Gl 


.) 

| 
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Strenge Lösungen zur Prandtlschen 'T'heorie 
der tragenden Linie. 


Von Harry Schmidt in Köthen. 


Vorbemerkung. Die Aufgabe, bei vorgegebener Tiefenverteilung und vorgegebener An- 
stellwinkelverteilung die zugehörige Zirkulationsverteilung längs einer tragenden Linie zu 
ermitteln, hat eine überaus reichhaltige Literatur hervorgerufen, auf deren kritische Würdigung 
ich demnächst an anderer Stelle ausführlich zurückkommen werde. 

In der vorliegenden Arbeit handelt es sich lediglich darum, erstens aus den physikalischen 
Hilfsvorstellungen (Prandtlsche Theorie) unter Beschränkung auf die durch sie allein 
bedingten Voraussetzungen strenge, d.h. durch keinerlei mathematisch nieht voll geklärtes 
Vorgehen beeinflußte Folgerungen zu ziehen, zweitens im Anschluß hieran für spezielle 
Klassen von Tiefenverteilungen bei beliebiger Anstellwinkelverteilung die Lösung des Problems 
in geschlossener Form tatsächlich herzustellen, und drittens die allgemeine Erledigung 
‚ler Aufgabe auf ein unendliches lineares Gleichungssystem zu reduzieren, das der Erh. Schmidt- 
sehen Auflösungstheorie zugänglich ist. 

Ihren praktischen Nutzen haben die hier erhaltenen Resultate bereits insofern erkennen 
lassen, als sie die Möglichkeit bieten, irgendein technisch interessierendes Beispiel im Sinne 
einer Fehlerabschätzung in eine Manmnigfaltigkeit streng beherrschter Typen einzubetten. 

Einen besonderen Genuß bereitete mir die Feststellung, daß die mathematischen Konse- 
quenzen der P randtlschen Ansätze eine innere Folgerichtigkeit aufdeeken, wie sie wohl nur 
wenigen Theorien der mathematischen Physik eigentümlich ist. 


Inhaltsübersicht. Im ersten Abschnitt werden die allgemeinen Grundlagen der Theorie 
der tragenden Linie kurz zusammengestellt, die nach L. Prandtl zu einer hochgradig 
singulären Integro-Differentialgleichung für die gesuchte Zirkulationsverteilung führen. 


Im zweiten Abschnitt wird die von E. Trefftz angegebene potentialtheoretische 
Formulierung des Problems dargelegt. 


Im dritten Abschnitt wird gezeigt, daß durch einfache Fortsetzungsvorschriften für die 
Tiefen- und Anstellwinkelverteilungen eine Reduktion auf das dritte Randwertproblem der 
Potentialtheorie für das Außengebiet des Einheitskreises möglich ist. 


Im vierten Abschnitt wird auf Grund einschlägiger potentialtheoretischer Ergebnisse 
unter Voraussetzung stückweise stetiger Tiefen- und Anstellwinkelverteilungen die Existenz 
einer einzigen, längs der ganzen tragenden Linie stetigen Zirkulationsverteilung bewiesen, die 
sich nach Übertragung auf den Einheitskreis in eine daselbst absolut und gleichmäßig kon- 
vergente Fourierreihe entwickeln läßt. 


Im fünften Abschnitt werden einige in der älteren Literatur bereits vielfach benutzte, 
auf die beiden ersten Fourierkoeffizienten der Zirkulationsverteilung bezügliche Aussagen 
sichergestellt. 


Im sechsten Abschnitt wird die strenge Lösung des im dritten Abschnitt erhaltenen 
kandwertproblems für den beliebig verwundenen elliptischen Flügel auf potentialtheoretischem 
Wege in geschlossener Form gewonnen und über einige Anwendungsergebnisse kurz berichtet. 


Im siebenten Abschnitt wird auf funktionentheorischem Wege die strenge Lösung für 
eine spezielle Klasse nieht-elliptischer Tiefenverteilungen bei beliebiger Verwindung in 
geschlossener Form hergeleitet. Überdies wird erwähnt, daß die Entwiecklungskoeffizienten 
der zugehörigen Potentialfunktion einem unendlichen linearen Gleichungssystem genügen, das 
nach der von Erh. Schmidt angegebenen, auf dem Orthogonalisierungsprozeß beruhenden 
Methode explizit gelöst werden kann. 


Im achten Abschnitt wird die Übertragbarkeit der im siebenten Abschnitt benutzten 
[unktionentheoretischen Methode auf allgemeinere Grundrißfamilien hervorgehoben; überdies 
werden einige Bemerkungen über die Approximation praktisch in Betracht kommender Tiefen- 
verteilungen durch solche Grundrisse angeknüpft, die durch das dargelegte Lösungsverfahren 
exakt erfaßbar sind. 


Im neunten Abschnitt wird der Nachweis erbracht, daß sich für die Fourierkoeffizienten 
der Zirkulationsverteilung irgendeiner stück weise stetigen Tiefenverteilung bei beliebiger Ver- 
windung ein unendliches lineares Gleichungssystem ableiten läßt, das der Erh. Schmidtschen 
Auflösungstheorie zugänglich ist. 


| 
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1. Die Prandtlsche Integro-Differentialgleichuug. Die Grundlage der von L. Prandtl'; 
entwickelten Tragflügeltheorie bildet der Ersatz einer senkrecht zur (uerachse des Flugzeugs 
mit der zeitlich konstanten Geschwindigkeit angeblasenen 'Tragfläche von der Spannweite 
durch eine tragende Linie, d.h. durch ein geradliniges Wirbelfadenstück von veränderlicher. 
in den beiden Endpunkten auf den Betrag Null absinkender Wirbelstärke, von dem ein ebenes, 
vemäß Abb. 1 sieh einseitig ins Unendlichferne erstreckendes Wirbelband abgeht. Die 


- Wirbelband 
Abb. 1. 
Intensität dieses Wirbelsystems wird durch eine auf dem Intervall ==, der.-Achse 
definierte, den Bedingungen | 
| ji 


eenügende Zirkulationsverteilung 7'(@.r) charakterisiert. Als Repräsentanten der Tiefenverteilung 
der Tragfläche bzw. der gegen die Strömungsrichtung im Unendlichen gemessenen Verteilung 
der Anstellwinkel der Nullauftriebsriehtungen der einzelnen Flügelschnitte sind auf dem 
eenannten Intervall zwei Funktionen F= bzw. a== vorgegeben, während als einzige 
Profileigenschaft unter Voraussetzung eines linearen Zusammenhangs zwischen dem Auftriebsbei- 
wert ce, und dem Anstellwinkel a die Konstante 


da 


de, 


Un = 


in die Theorie eingeht’). 

Längs der tragenden Linie wird durch das Wirbelband eine gegenüber der Grund- 
strömungsgeschwindigkeit » als klein anzusehende Abwärtsgeschwindigkeit w (#) von der Größe 


b|2 
| ve ds 
-bj2 


hervorgerufen, falls man sich auf solehe Zirkulationsverteilungen 7'(.*) beschränkt, für die der 


Cauchysche Hauptwert des auf der rechten Seite von (1,2) auftretenden uneigentlichen 
Integrals existiert und einen positiven Wert besitzt. Dies hat zur Folge, daß die Anblase- 
h h 

richtung um den auf dem Intervall variierenden induzierten Anstellwinkel a; (#) 


gesenkt wird, für den wegen seiner Kleinheit mit praktisch ausreichender Annäherung 


mir) 


zesetzt werden darf. Die Verteilung des sogenannten effektiven Anstellwinkels «a. geht 
demnach aus der geometrischen Anstellwinkelverteilung a(#) durch Subtraktion der Funktion 
a;la) hervor. 

Nach der Kutta-Joukowskischen Theorie der beiderseits unbegrenzt breiten Trag- 
fläche ergibt sich die längs der ganzen Spannweite konstante Zirkulation 


wenn «t den Radius desjenigen Kreises bedeutet, auf dessen Außengebiet das Außengebiet der 
Protilkontur (bei geeigneter Normierung der Abbildungsfunktion im Unendlichfernen) umkehrbar 


I) L. Prandtl: Göttinger Nachr., math.-phys. Kl 1918, S. 451. 


=) Allgemeiner darf a’. als eine wesentlich positive, anf dem Intervall — ,„ <x< , vorgegebene Funktion 


„(x) aufgefaßt werden: alsdann behalten sämtliehe folgenden Darlegungen Gültigkeit, wenn (x) dureh (X) 
ersetzt wird. 
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eindeutig und konform ‘abgebildet werden kann. Überdies liefert die ebene Theorie einen 
konstanten Auftrieb 4 pro Einheit der Tragflächenbreite vom Betrage 


om 
für den in der Praxis die Darstellung 


A=r,ta)- 
mit o als Luftdiehte und mit Fals Profiltiefe (Maximalabstand zweier Punkte der Profilkontur) 
üblieh ist. Demnach gilt 


woraus weren (1.3) bei Beschränkung auf kleine Anstellwinkel «a 


Ca a, 
mithin 
de, 
da 
oder 
3 
| la . . . . . . . . . (1.53) 
folet. 


Die mit (1,31) aus (1,5) hervorgehende Beziehung 


[4 


wird nunmehr in der Theorie der tragenden Linie zu der Forderung verallgemeinert, daß auf 


dem Intervall — die Relation 


erfüllt sein soll. Führt man hierin (1,21) mit (1,2) ein, so erhält man die Integro-Differential- 
gleiehung 


—b12 


die den von L. Prandtl angegebenen Zusammenhang zwischen den drei Funktionen Fr), 
ale) und /'(#) zum Ausdruck bringt. Soll also nach Vorgabe einer technisch brauchbaren, 
Tiefenverteilung # (#*) sowie einer beliebigen, praktisch in Betracht kommenden Anstellwinkel- 
verteilunge ar) die zugehörige Zirkulationsverteilung 7'(@r) ermittelt werden, so ist die 
Funktionalgleichung (1,5) nach 7'(r) aufzulösen. Die sehr erheblichen mathematischen 
Schwierigkeiten, die der Bewältigung dieser Aufgabe entgegenstehen, liegen wesentlich in der 
Tatsache begründet, daß (1,5) eine hochgradig singuläre Gleichung darstellt, die den 
üblichen Methoden der Integralgleiehungstheorie keineswegs ohne weiteres zugänglich gemacht 
werden kann. 


2. Die Trefftzsche Formulierung als potentialtheoretische Randwertaufgabe. Zu einer 
ebenso eleganten wie prinzipiell bedeutsamen Formulierung der uns interessierenden Problem- 
stellung ist E. Trefftz*) durch die Erkenntnis gelangt, daß sich die Zirkulationsverteilung 
und die induzierte Abwärtsgeschwindigkeit durch eine Potentialfunktion zur 
Darstellung bringen lassen. Denkt man sich nämlich die @y- Ebene (vgl. Abb. 1) längs der 
die tragende Linie repräsentierenden Strecke 


>) E.Trefftz: Math. Ann. 82, S. 306, 1921: ZAMM 1, S. 206, 1921. Vgl. auch „Vorträge aus dem Gebiete der 

Iydro- und Aerodynamik (Innsbruck 1922)*, S. 34: Berlin 1924. Eine etwas abweichende Darstellung der einschlägigen 

Zusammenhänge habe ich im 7. Kapitel meiner „Aerodynamik des Fluges“ (Berlin 1929) gegeben. 
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aufgesehnitten, und kennzeichnet man (vgl. Abb. 2) den Übergang aus jener Ebene auf das 
obere Sehlitzufer 


b 
bzw. auf das untere Sehlitzufer 
h 


kurz durch Anhängen des Index o bzw. u, so gilt 


une 


wobei unter g Gr, y) eine im ganzen Schlitzgebiet einschließlich des unendlichfernen Punkts 
reguläre Potentialfunktion zu verstehen ist. hiermit nimmt aber die Forderung (1,4) die Form 


an, und so erweist sich die Bestimmung der bei vorgegebener Tiefenverteilung f(x) und vor- 
gegebener Anstellwinkelverteilung «a resultierenden Zirkulationsverteilung als auf das 
potent’altheoretische Randwertproblem reduzierbar, eine in der ganzen längs der 
Strecke (2,1) aufgesehlitzten «Ebene reguläre Potentialfunktion y) zu ermitteln, die 
beim Übergang auf die Schlitzufer (2,11) und (2,12) den beiden Bedingungen (231) und (2,32) 
genügt. 


(0) (0) 
ty 
[ 
(0) Ä (0) 
bo u) z 
| (u) (u) 
r=7 
Abb. 
Durch die Transformation 


wird das Detinitionsgebiet der Potentialfunktion 2 (#@, y) auf das Außengebiet r>1 des Ein- 
heitskreises #1 der Ebene der komplexen Variablen [== r-e‘” derart konform abgebildet, 
daß dem oberen Sehlitzufer (2,11) der obere Halbkreis 


und dem unteren Schlitzufer (2,12) der untere Halbkreis 

entspricht (vel. Abb. 3). Da gemäß (2,4) 


wird, so ıst längs des Schlitzes 


j| für y-+0, \ 
fü 


mit \ 
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zu setzen; ferner ergibt sich nach einfacher Rechnung 
5) 


„ b-sınd 


or 
so daß (2,31) unter Beibehaltung der mit Berücksichtigung von (2,43) zu verstehenden Be- 
zeiehnungen f(x) und a (x) in 


sowie (2,32) In 


übergeht; das Argument 9 bleibt dabei der Einschränkung 


unterworfen. Auf diese Weise ist also unsere Aufgabe auf das Problem zurückgeführt, eine 
im Außengebiet r>1 des Einheitskreises r = 1 überall reguläre Potentialfunktion y (r, d) zu 
bestimmen, die längs des oberen Halbkreises (2,41) die Bedingung (2,51) und längs des 
unteren Halbkreises (2,42) die Bedingung (2,52) erfüllt. 


3. Zurückführung auf das dritte Randwertproblem der Potentialtheorie. Schreibt man 
(2,52) in der Form 


so muß auch 
werden. Da nun überdies 
[sin = — [sin 9]. . (3,12) 


eilt, so läßt sich offenbar für das zweite auf der linken Seite von (2,51) auftretende Glied 
das Bestehen der zu (3,1) analogen Relation 


sin® Or|, 


(ur) 
sind) Or 


höchst einfach dadurch erreichen, daß die Festsetzung 


getroffen, mithin die gemäß (2,43) zunächst nur auf dem oberen Halbkreis (2,41) vorgegebene 
Tiefenverteilung zu einer auf dem ganzen Einheitskreis r 1 definierten ungeraden 
Funktion erweitert wird. Um alsdann auch für die rechte Seite von (2,51) die durch 


a = — a 
zum Ausdruck gelangende Eigenschaft zu erzwingen, brauchen wir lediglich 


zu wählen, d.h. also die Anstellwinkelverteilung im Gegensatz zur Tiefenverteilung als 
eine gerade Funktion auf den unteren Halbkreis fortzusetzen. Dies alles besagt aber nichts 
anderes, als daß die Einführung der Vorschriften (3,21) und (3,22) die beiden von der ge- 
suchten Potentialfunktion 9 (r, 9) zu verlangenden Forderungen (2,51) und (2,52) zu der einen 
längs des ganzen Einheitskreises #1 zu erfüllenden Randbedingung 


l 
mit 


zusammenzufassen gestattet. Damit ist jetzt das Fundamentalproblem der Theorie der tragen- 
den Linie auf das dritte Randwertproblem der Potentialtheorie für das Außen- 
gebiet des Einheitskreises zurückgeführt. 


| 
j 
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4. Existenz einer und nur einer stetigen, in eine absolut und gleichmäßig konvergente 
Fourierreihe entwickelbaren Zirkulationsverteilung. Sınd «die zur Flügelmitte z#—=0 symme- 
trische Tiefenverteilung f(x) und die Anstellwinkelverteilung «a(.“) als auf dem Intervall 

h 
9 #0, stückweise stetige Funktionen von . vorgegeben, so gehen sie mit der 


Substitution (3,31) unter Beachtung von (3,21) und (3,22) in längs des Einheitskreises r — | 
stückweise stetige Funktionen von # über. Da ferner auf dem offenen Intervall 


stets >O vorausgesetzt werden darf, so stellt auch 
sin) 
) — 
T(9) (4,1) 


eine längs r = | stück weise stetige, wegen (3,12) und (3,21) beständig positive Funktion von d 
dar, falls als Funktion von für und d.h. also an den beiden Flügelenden, 
je eine einfache Nullstelle besitzt. Wird überdies g (r, d) als eine beschränkte, für 
reguläre und mit beschränkter Ableitung versehene Potentialfunktion voraus- 
zesetzt, die für r>1, von den Unstetigkeitsstellen der Funktionen a(x) und T(#) abee- 
sehen, der Randbedingung (3,3) genügt, so existiert nach bekannten potentialtheoretischen 
| 
r>1 auf jedem von Unstetigkeiten der Funktionen ar) und T(#) freien Bogen des Einheits- 
kreises eleichmäßig. Fassen wir nun die somit längs r=1 überall beschränkte und, von 
den Unstetigkeitsstellen der Funktionen ale) und T(9) möglicherweise abgesehen, stetige 


Ergebnissen?) stets eine und nur eine derartige Funktion, und es konvergiert für 


Funktion als Vorgabe’) des zweiten Randwertproblems für das Außengebiet des 


Kinheitskreises r—=1 auf, so erscheint die Funktion @(r,d) als Lösung dieses Randwert- 
problems in der Form des logarithmischen Potentials einer einfachen Belegung, deren Diehte u (#) 
unbedingt integrierbar und, möglicherweise wiederum von den Unstetigkeitsstellen der Funktionen 
ala) und T(9) abgesehen, längs = 1 stetig ist: die Lösung (r, selbst ist längs r=|1 
überall stetig"). Damit folgt aber aus unserer Randbedingung (3,3) die stückweise Stetig- 
Orlr>i 

einer einfachen Belegung mit der jetzt als längs des ganzen Einheitskreises stück weise 
stetig’) erkannten Dichte «(9 Demnach gilt mit beliebig kleinem Z>0 und mit 
M gleichmäßig in auf dem Intervall — a = die Abschätzung‘) 


keit von und damit wiederum die Darstellbarkeit der Funktion g (r, 9) als Potential 


Per, 9 >ı - Mr >ı) <M- Al): 


dureh die schließlich die Entwickelbarkeit der stetigen Funktion [| (r, #)Ir>ı in eine auf dem 

Intervall absolut und gleichmäßig konvergente Fourierreihe sichergestellt 

wird"). Wegen (2,21) ist daher für alle den genannten Bedingungen genügen- 

den Tiefenverteilungen f(r) und Anstellwinkelverteilungen a(r) die Existenz 

einer einzigen Zirkulationsverteilung /'(r) bewiesen, die auf ihrem ganzen 
h h 


Definitionsintervall 


> 5: „ stetig ist, und die sieh als Funktion von 


in eine absolut und gleichmäßig konvergente Fourierreihe entwickeln läßt. 

5. Die Fourierkoeffizienten der Zirkulationsverteilung. Da die durch (3,3) eindeutig 
festgelerte, für > I reguläre Potentialfunktion (r, eine ungerade Funktion von sein 
muß, so kann sie in eine für r > I absolut und gleichmäßig konvergente Reihe von der Form 


A 


L.. Liehtenstein: Neuere Entwicklung der Potentialtheorie. Konforme Abbildung. Ene. d. math. 


Wiss. 11, 3, I, Artikel II C3, Nr. 28, S, 279/280. 
"T 
Die Erfüllung der Bedingung | ist dabei dureh (3,11) gewährleistet. 
Vel. L. Liehbtenstein:]. Nr. 24, S. 202. 
Wird 7 7) als eine stetige Funktion von * vorausgesetzt, so folgt dieses Ergebnis einfacher aus den Dar 
gungen von T.Carleman, Diss. Upsala 1916 ($ 3). 
>, Vgl. L. viehtenstein:l. e Nr. 6, S. 201, insbesondere Fußnote 8). 
Vel. hierzu hinsiehtlieh der Konvergenz W. Rogosinski: Fouriersche Reihen 
(Sa ne Goschen Bd. 122), S. 59: hinsiehtlich der absoluten Konvergenz vel.S. Bernstein: C. R. 158, S. 1661, 1914. 
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ontwiekelt werden. Da andererseits die Fourierreihe des Randwerts [9 (r, A)]->ı konver- 
viert. so darf auf (irund des Abelschen Grenzwertsatzes 
| 
reschrieben werden, und die in (5,1) auftretenden Koeffizienten A, sind daher als die 
Fourierkoeffizienten 
A, = (ir, > ı sin (v (2,12) 
der Funktion (r, anzusprechen. 
Der der Zirkulationsverteilung /'(#*) zugeordnete Gesamtauftrieb A des Tragflügels ist 
bei Beachtung von (2,21) dureh die Beziehung 


b|2 
Az= \ I h: \ (r, > ı 
oder gemäß (5,12) dureh 
orba w 


vereben, und wenn man den in der Technik üblichen Ansatz 


- 4 


benutzt, in dem F den Flächeninhalt der größten Flügelprojektion mißt, so gilt 


Ebenso findet man für das Rollmoment (Moment der Auftriebsverteilung um die Längsachse 
des Flugzeugs) Mo, wegen 


b/2 
oob? 
M,=or: x-I(e)de= € 4 (r, sin (20) dd 
bj? 
einfache Relation 
Ma="—— (5,3) 
die mit dem Ansatz 
l 
zu 
A 
führt, unter 
A= 


das Seitenverhältnis der Tragtläche verstanden. 


6. Der elliptische Flügel. Die tatsächliche Herstellung der Lösung g(r,d) des 
kandwertproblems (3,3) gestaltet sich besonders bequem für den elliptischen Flügel, bei 
dem im Einklang mit (3,21) 


zu setzen ist. Alsdann kann nämlich (3,3) mit Benutzung der Abkürzung 


(6,11) 


| 
| 
| _ 
| 
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auf die einfache Form 


. 
gebracht werden, und da mit y (r, 9) auch von i eine im Unendlichfernen reguläre Potential- 


funktion ist, so läuft die Randwertaufgabe (6,12) auf ein Dirichletsches Problem für das 
Außengebiet des Einheitskreises hinaus, das mit Hilfe des Poissonschen Integrals 


h 
I ob —1)-siny”-a | „ cos 
\ d 


2r-cos(d — y*) (6,2 
«dureh 


zelöst wird"). Nun hat die gewöhnliche homogene Differentialgleichung 
1 du | 
k dr 

als partikuläres Integral die Funktion 

u(r)=rk, 
so daß, wenn F in (6,21) als Parameter aufgefaßt wird, mit Anwendung der Methode der 
Variation der Konstanten als einzige der Differentialgleichung (6,21) genügende, im Unendliech- 
fernen reguläre Potentialfunktion der Ausdruck 


k-rk. 


o(r,d)= 


r 
resultiert. Ihiermit ist demnach die Lösung des Problems der tragenden Linie für den 
beliebig verwundenen elliptischen Flügel in geschlossener Form aufgefunden. 

Ist «die Anstellwinkelverteilung a als eine längs der ganzen Spannweite - 


eanz-rationale Funktion von vorgegeben, gilt also 


m 
; ! [7 
) 
adlır) \ \ E (m 0), 


so ist die Bedingung (3,22) automatisch erfüllt, und die Auswertung von (6,2) läßt sieh mit 
Benutzung «der bekannten trigonometrischen Identitäten 


| 


de 


une 

[2.00 


auf «ie mit Hilfe des Residuensatzes leicht zu bereechnenden Integrale 


und > 1) 
\ 2 fi d TE 


10) Vgl. hierzu MH. Schmidt: Math. Annalen MH2, S. 322, 1936, wo auf analogem Wege die Lösung eines statischen 
Problems für die gelenkig gestützte Kreisplatte auf das entsprechende Problem für die eingespannte Kreisplatte zu- 
rüekgeführt wird. 
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reduzieren. Für den unverwundenen Flügel 


wird auf diese Weise 


sın 

womit (6,5) sofort 


ergibt, während man bei der linearen V erwindung 


zunächst er sin (29) 
und damit dann sin @9) 
y(r,d)= 1: (k +2) 


erhält. 

Auch für stückweise ganz-rationale Anstellwinkelverteilungen a (.“) können unter 
der Voraussetzung rationaler Werte der Konstanten % geschlossene, durch elementare 
Funktionen ausgedrückte Lösungen erhalten werden. Als Beispiel führe ich hier einige von 
Herrn K. Bauseh DVL in einer noch nieht veröffentlichten Arbeit erhaltene Resultate an, 
die sich auf einen an einem Ende um den festen Winkel 5 gemäß Abb. 4 verwundenen 


| 
d b 
"7 %o 0 | 7 
| N 2 
| 


Abb. 4. 


elliptischen Flügel beziehen, der überdies zur Erzielung eines verschwindenden Gesamtauf- 
triebs A längs der ganzen Spannweite um den festen Winkel «, gedreht ist. Mit 


ergibt sieh dabei die vom Seitenverhältnis A unabhängige (!) Relation 


ß > (6,61) 
die für wandernde Sprungstelle in Abb. 5 dargestellt ist. Abb. 6 veranschaulicht bis auf 
konstante, durch technische Gesichtspunkte bedingte Faktoren für 

2b. 

die der Z-Verteilung proportionale Größe die e„-Verteilung, die a;-Verteilung, die 
Verteilung und die dem örtlich induzierten Widerstand proportionale Funktion 
während in Abb. 7 der Rollmomentenbeiwert eo bei beliebigem Seitenverhältnis „I sowie 
in Abb. S bzw. Abb. 9 die Beiwerte (e,.;)ea-0 bzw. (Oyr)aa=0 des bei verschwindendem Ge- 
samtauftrieb infolge der Verwindung induzierten Gesamtwiderstands bzw. des zugleich auf- 
tretenden Wendemoments für =3 zur Wiedergabe gelangen. 

Für beliebige Anstellwinkelverteilungen a(x)'bei beliebigem % führt die für r > 1 

absolut und gleichmäßig konvergente Reihenentwieklung 


\ cos [r (9 
= 
v=] 
unter Beachtung von (3,22) zu 
A, sin(r 


v=] 
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>: h- \ sın cos y*)-sin . (6,72) 
0 


so daß gemäß (6,5) 


\’ 4A, sind) 


( 
v—] 
wird. 
da 
0 0,5 b/2 10 
Abb. ». 
Ofr 
+ 
= 
0 0,5 da 7,0 
b/2 
7. 
0.057 
| 
| 
ala 3 
0.025} 
In 
0 ba. | 
i 
| 
% 
7 
| 
D 
0,005 w | 
| 
Abb. ®. | Abh. 6. 


7. Eine spezielle Klasse nicht-elliptischer Tiefenverteilungen. Um nunmehr zu nicht- 
elliptischen, zur Flügelmitte symmetrischen Tiefenverteilungen überzugehen, wählen wir 
im Kinklang mit (3,21) zunächst den Ansatz 


mit 7,,< 1, für den wir bei beliebiger Verwindung «(.“) das durch (3,3) gekennzeichnete 
kandwertproblem in geschlossener Form auf einem Wege lösen wollen, der auch unter 
Zugrundelegung allgemeinerer Klassen von Tiefenverteilungen (vgl. Abschnitt S) olıne weiteres 
anwendbar bleibt. 
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Setzen wir zur Abkürzung 


» 


so geht (3,3) mit (7,1) unter Beibehaltung von (6,11) in 


| 
über. Mit Benutzung der komplexen Variablen 
führen wir sodann durch 
v=] 


diejenige bis auf die willkürlich wählbare, hier der Einfachheit halber als Null angenommene 
(imaginäre) Konstante F(x) eindeutig festgelegte, im Außengebiet 2) = r>1 des Einheits- 
kreises r— 1 überall reguläre analytische Funktion von 2 ein, deren Realteil mit der gesuchten 
Potentialfunktion y(r, 9) übereinstimmt. Dann wird wegen 


v=] 
und 
>| 


die linke Seite der Gl. (7,2) der längs r=1 angenommene Randwert des Realteils R {f(z2)} der 
Funktion 


mit der Potenzreihenentwicklung | 
Demnach stellt 
it, 


eine auch im Unendlichfernen reguläre Funktion von z dar, für die 


it 


-A,-2 


it, _ - 


gilt. Addiert man also auf beiden Seiten der Beziehung (7,2) die Funktion 


T 


2 
2k 


A, -sın 


so erhält man als Randbedingung für die für !z!>1 überall reguläre Funktion (7.34) 
die Forderung 


-A,-sınd+ (29:40), 


aus der bei Beachtung von (7,35) mit Benutzung des Poissonschen Integrals sofort 


it l i 
mit 


| 
8% 
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folet. Führen wir (7,33) in (7,36) ein, so ergibt sich für die Funktion F(z) die lineare Differen- 
tialgleichung 


Du 


+3 [2 + ;) .2: F’(e) = (741 


(deren der Bedingung 


eenügendes Integral 
mit 
2 
und mit 
h 


lautet. Diese Funktion F(z) muß aber im allgemeinen (d. h. für irgendwie gewählte Werte 
von |, und 4,) in den beiden dem Außengebiet |z!>1 des Einheitskreises angehörenden 
Stellen 


Singularitäten aufweisen, deren Charakter dureh die in (7,52) verzeichnete Zahl x bedingt 
wird. Demnach sind die Konstanten A, und A, dureh die Forderung festzulegen, daß sich 
die Funktion F(z) auch in den Punkten z2=2z, und z=z, ihres Definitionsgebiets regulär 
verhalten soll. 


Indem wir zunächst die Stelle 2=ia ins Auge fassen, schreiben wir (7,5) in der Form 


) A, f, 


woraus sofort ersichtlich ist, daß notwendig 


a? 


sein muß. Diese Bedingung ist aber auch hinreichend für reguläres Verhalten von F(z) 


in 2 7a, «da für eine passend gewählte Umgebung jener Stelle 
| 


gesetzt werden darf, mithin dort die reguläre Entwicklung 


ia) 


gilt. Nach entsprechendem Vorgehen bei Vertauschung von # mit  i bekommt man gemäß (7,6) 
(die beiden linearen Gleichungen 
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2 1 T, (3? + a?)* 
ia “ia 
und 
-ia +) ia +) 
da 
wird. 


Mit den als Lösung der Gl. (7,61) und (7,62) ermittelten Konstanten A, und A, wird 
dureh (7,5) die der Tiefenverteilung (7,1) bei beliebiger Verwindung a (.*) zugeordnete Funktion 
F(z) und damit gemäß (7,3) auch die zugehörige Potentialfunktion y (r, 9) in geschlossener 
Form dargestellt. Die Auswertung einschlägiger praktisch interessierender Spezialfälle hat 
Herr J. Gehlen in Angriff genommen. 


Wird die in (7,4) verzeichnete Funktion @ (z) gemäß 


in eine für >1 absolut und gleiehmäßig konvergente Potenzreihe entwickelt, so liefert 
die Differentialgleichung (7,41) nach Eintragung der in (7,3) angegebenen Entwicklung von (2) 
dureh Koeffizientenvergleich das unendliche Gleichungssystem 


(vl) 
(v„— 2)r, k (v+2)r, 


das als einzige mit konvergenter Quadratsumme versehene Lösung die Koeffizienten A, von 
F(z) ergibt. Hierüber sowie über die explizite Berechnung der A, aus (7,71) nach der von 
Erh. Sechmidt'') angegebenen, auf dem Orthogonalisierungsprozeß beruhenden Methode wird 
Herr H. Schubert demnächst berichten. 


8. Allgemeinere Grundrißfamilien. Die Übertragung der im vorhergehenden Abschnitt 
dargelegten Lösungsmethode auf den Fall von solchen Tiefenverteilungen f(x), die in Ver- 
allgemeinerung von (7,1) durch den der Bedingung (3,21) genügenden Ansatz 

. Y fi 
mit erfaßt werden können”), bereitet selbstverständlieh keinerlei prinzipielle 
Schwierigkeiten, da ja in Analogie zu (7,32) 


Y f \ 


eilt. 


11) Erh. Schmidt: Rend. Palermo 35, 8. 53, 1908. Vgl. hierzu E. Hellinger und ©. Toeplitz: Integral- 
sleichungen und Gleichungen mit unendlich vielen Unbekannten. Ene. d. math, Wiss. 11, 3, 11, Artikel II C 13, Nr. 19. 


n 
1+ > 
1?) Noch allgemeiner darf, wie man leicht sieht, f(x) fo-sin"» genommen werden. 
m 
1-+ 10 COS (2u%) 
u=1 
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Um jedoch die zur Integration der alsdann an Stelle von (7,41) auftretenden linearen 
Differentialgleichung benötigten Integrale 


d; 
| 
) 


bequem auswerten zu können, sind der Anwendung von (S,1) gewisse praktische Grenzen 
gesetzt'’). Nehmen wir beispielsweise an, daß man einen Rechteckflügel durch Tiefen- 
verteilungen von der Form (8,1) derart zu approximieren wünscht, daß im Einklang mit Ge- 
pflogenheiten der Flugteehnik und neuerdings auch bei Windkanalversuchen zwecks Ver- 
meidung örtlicher Ablösungsgefahr eine geringfügige Abrundung der beiden äußeren Ecken 
resultiert. Auf Grund geläufiger trigonometrischer Identitäten kann nun (8,1) dureh 


. Y 
I’ 
oder wegen 
2x 
cos 
auch durch 
„ 


ersetzt werden, so daß t(x) eine lineare Superposition von Funktionen der Form 


darstellt. Diese Funktionen sind für n=1bisn =5 in Abb. 10 über dem Intervall 0=&2<1 
graphisch wiedergegeben: man erhält so eine nützliche Handhabe für ein zwecekmäßiges Vor- 
gehen bei der Durchführung der Aufgabe, eine bereits erreichte Approximation einer an- 
eestrebten Umrißform dureh Hinzufügung eines neuen Summanden zu verbessern. Man wird 
sich dabei aber aus den zuvor erwähnten Gründen die Einschränkung aufzuerlegen haben, 
daß die Nullstellen des in (5,2) auftretenden Polynoms in £ mit mehr oder minder großem 
Rechenaufwand explizit zu ermitteln sein müssen: ein Beispiel, für das dies zutrifft, ist in 
Abh. unter der Annahme veranschaulicht. 


m) 
| _7\ n=4 


| 


05 70 


Abb. 10, Abb. 11. 


'3) Von dieser Schwierigkeit unberührt bleibt die Auflösbarkeit des jetzt an Stelle von (7,71) sich ergebenden 
Gleichungssystems nach der Erh. Sehmidtschen Methode. 


H 
n=7 
n=2 
05 | 10% 
| 
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Sehr viel zweekmäßiger erscheint es daher, zum Aufbau von Tiefenverteilungen #(.r), 
uf die unsere Methode übertragbar ist, neben der in (8,2) ausschließlich benutzten additiven 
Verknüpfung auch multiplikative Bildungen heranzuziehen. So hat Herr M. Zoll einen 
Rechteekflügel durch den Ansatz 

(140333). +0,25) 0,835 5%) (+1,11 


so weit approximieren können, daß längs zirka 95"/o der Spannweite die Abweichungen von 
‚ler Wurzeltiefe f, den Betrag 0,005 -#, nicht überschreiten (vgl. Abb. 12). 


Abb. 


9. Beliebige stückweise stetige Tiefenverteilungen. Schreiben wir (3,3) mit Benutzung 
von (4,1), (6,11) und (7,11) in der Form 


| 
und führen neben (7,3) durch 


diejenige bis auf eine imaginäre (hier gleich Null gesetzte) Konstante festgelegte, für 2, >| 
überall reguläre Funktion 


der komplexen Veränderlichen z=r-e'” ein, deren Realteil U(r, d) der Randbedingung | 

genügt, so kann (9,1) nach Addition von kt,-[V (r, d)-y(r, Od), >ı auf beiden Seiten wegen 


dureh 


ersetzt werden, woraus aber sofort die lineare Differentialgleichung 


mit 


— 


(081) 


| 
| 
| 
| | 
| 
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folet. Da nun T(#) wegen (3,21) eine gerade Funktion von # ist, so wird gemäß (9,11) 


— 


mit 
| 


u 
ferner gilt weren (7,3). d. h. 


\’iA, 


die absolut und gleichmäßig’) konvergente Entwicklung 
_ 
mit der man unter Beachtung von 
sowie von 
(51) die Beziehung 
It, 
2 
und 
2 


erhält. Trägt man die Entwieklungen (9,5), (9,4) und (9,6) in (93) ein, so ergibt sich durch 


Vergleich der Koeffizienten von —, ein unendliches lineares Gleichungssystem für die Unbe- 


kannten ‚I, das «den folgenden Bedingungen genügt: 


I. Für die gesuchte Lösung A. ist eine konvergente Quadratsumme zu fordern; 
>. in der Koeftizientenmatrix sind die Elemente jeder Zeile reell und besitzen eine 


konvergente (Juadratsumme. 


Überdies hat Herr H. Schubert im Zusammenhang mit seinen im vorhergehenden Abschnitt 
erwähnten Untersuchungen gezeigt, daß auch 


>. je endlich viele Zeilen der Koeffizientenmatrix des Gleichungssystems linear 
voneinander unabhängig sind. 


Dammt ıst dann das resultierende Gleichungssvstem als ein solches vom Erh. Schmidtschen 
Typus erkannt, für das eine vollständige Auflösungstheorie fertig vorliegt '’) '®). 

Mit aufriehtiger Freude hebe ich am Schluß dieser Arbeit hervor, daß mir die Möglich- 
keit zu einer eingehenden Beschäftigung mit den hier behandelten Fragen durch die Deutsche 
Versuchsanstalt für Luftfahrt in Berlin-Adlershof geboten wurde. 656 


14, Vgl. hierzu S. Bernstein: EC. R. 158, S. 1661, 1014. 


15) Verl. hierzu wieder E. Hellinger und O. Toeplitz, Integralgleichungen und Gleiehungen mit unendlich- 


vielen Unbekannten. Ene. d. math. Wiss. II, 3, II, Artikel II C 13, Nr. 19. 


16) Auf die praktischen Auswirkungen dieser Tatsache wird Herr H. Schubert zurückkommen. 


| 
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/um Übergang von Unterschall- in Überschallströmungen. 
Von W. Tollmien in Göttingen. 


(Aus dem Kaiser Wilhelm-Institut für Strömungsforschung:.) 


adiabatisch verhält. Die Strömungsgeschwindigkeiten können dann aus einem Potential 
abzeleitet werden. Kompressible Potentialströmungen sind aus diesem Grunde genau durch- 
forseht worden. Besonders leistungsfähige Rechenverfahren liegen für den Fall vor, daß die 
Geschwindigkeit der stationären Strömung überall die Schallgeschwindigkeit übersteigt. Aber 
auch dann hat man kompressible Strömungen berechnen können, wenn die Schallgeschwindig- 
keit nirgends erreieht wird. Diese Unterschallströmungen sind qualitativ den inkompressiblen 
Strömungen recht ähnlich. In der Tat hat man dureh schrittweise Näherung von der volumen- 
beständigen Strömung her den Einfluß der Kompressibilität bei den Unterschallströmungen 
berücksichtigen können. Große Schwierigkeiten treten aber auf, wenn weder eine aus- 
gesprochene Überschall- noch Unterschallströmung vorliegt, sondern die Schallgeschwindigkeit 
von der Strömungsgeschwindigkeit durehschritten wird. Mathematisch haben diese Schwierig- 
keiten ihren Grund darin, daß die Differentialgleichung für das Strömungspotential beim 
Übergang aus dem Unterschall- in das Überschallgebiet vollständig ihren Charakter ändert. 
Im Unterschallgebiet ist die Differentialgleichung elliptisch, im Überschallgebiet hyperboliseh. 
Außerdem liegt die Grenze zwischen beiden Gebieten nieht etwa von vornherein fest, sondern 
muß bei der Lösung der Differentialgleichung erst mitbestimmt werden. Auch die Sehall- 
geschwindigkeit selbst, die für diese Grenze maßgebend ist, hängt von den Strömungs- 
geschwindigkeiten ab. 


A uch eine kompressible, reibungslose Strömung bleibt wirbelfrei, solange das Gas sich 


Trotzdem ist schon frühzeitig ein Versuch gemacht worden, den Durchgang der Strömung 
durch die Schallgesehwindigkeit an einem Beispiel zu studieren. In der Göttinger Dissertation 
von Th. Meyer') aus dem Jahre 1907 ist das Potential einer zweidimensionalen Strömung 
beim Durchgang dureh die Schallgesehwindigkeit in eine Potenzreihe nach kartesischen 
Koordinaten entwickelt. Freilich kann man über die Konvergenz dieser Entwickelung von 
vornherein nichts aussagen. Aber Th. Meyer erhielt, wenn er Glieder bis zum 6. Grade 
einschließlich berücksichtigte, gute Übereinstimmung mit seinen Experimenten in einem Kanal 
mit einer Einschnürung (Lavaldüse). Die Abb. I zeigt ein etwas aufgefrischtes Bild der 
Ergebnisse Th. Meyers. Die „-Achse steht senkrecht zur eingezeichneten .„-Achse, der 
Koordinatenursprung ist durch einen Nullenkreis markiert. Aus den Experimenten ist zu 
schließen, daß die Rechnungen etwa innerhalb eines Kreises mit dem Radius 0,3 um den 
Koordinatenanfangspunkt gültig sind. 


/ / \y 
[A679.1] 


Abb, I. Stromlinien und Machsche Wellen einer ebenen Strömung dureh einen Kanal mit einer Einsehnürung (Laval- 
dise) beim Durchgang dureh die Schallgesehwindigkeit nach Th. Meyer (1908). Die Schallgrenze ist gestrichelt. 


Viele Jahre später hat G. J. Taylor’) dasselbe Entwickelungsverfahren angewandt, 
ebenfalls um Aufschluß über das Verhalten einer Gasströmung beim Durchgang durch die 
Schallgeschwindigkeit bei der engsten Stelle eines Strömungskanals zu erhalten. S.G. Hooker‘) 


') Th. Meyer: Über zweidimensionale Bewegungsvorgänge in einem Gas, das mit Überschallgesehwindigkeit 
stromt. Mitteilungen über Forschungsarbeiten a. d. Geb. d. Ing.-Wesens, H. 62, Berlin 1908, Abschnitt ©. 
°),G.J. Taylor: The flow of air at high speeds past eurved surfaces. Reports and Memoranda Nr. 1381 of 
"he Aeronautieal Research Committee, 1930 - 1931, London (1952). 
») S.G. Hooker: The flow of a compressible liquid in the neighbourhood of the threat of a eorstrietion in a 
eireular windehannel. Proe. of the Royal Soc. A. Vol. 135, p. 498, 1932. 
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benutzte diesen Ansatz für den rotationssymmetrischen Fall und erhielt befriedigende Über- 
einstimmung mit Experimenten von T. E. Stanton’"). 

In einer anderen Abhandlung’) gibt G. J. Taylor beiläufig eine exakte Lösung mit 
Überschall- und Unterschallgeschwindigkeiten an. Es ist dies der Potentialwirbel mit genau 
derselben Geschwindigkeitsverteilung, wie sie in den Elementen der Hydrodynamik ohne 
Berücksichtigung der Kompressibilität gelehrt wird. Um einen kreisförmigen Kern mit Über- 
. schallgesehwindigkeit ist also das Unterschällgebiet gelagert. Taylor schlägt zwei verschiedene 
Wege ein, um aus «dieser einfachen Grundströmung zu anderen Strömungen zu gelangen, die 
auch Überschall- und Unterschallgebiete enthalten. Einmal wird eine rein radiale Geschwindig- 
keitsverteilung dem Potentialwirbel überlagert, so daß als Stromlinien eine Art Spiralen auf- 
treten. Wir werden uns mit dieser Strömung noch näher zu befassen haben®). An zweiter 
Stelle betrachtet Taylor eine schwache wellenförmige Störung des Potentialwirbels, was uns 
hier weniger interessiert. 

Damit ist die Aufzählung aller bekannten Lösungen unseres Problems erledigt. Es 
besteht demnach ein dringender Bedarf an exakten Lösungen des Potentialproblems kom- 
pressibler Strömungen, die Unterschall- und Überschallgeschwindigkeiten aufweisen. Solche 
sollen im folgenden gegeben werden. 

Wir beginnen mit der Aufstellung der Grundgleiechungen kompressibler stationärer 
Potentialströmungen. 


1. Gleichungen der kompressiblen stationären Potentialströmungen. In einer stationären 
wirbelfreien Strömung mit dem Geschwindigkeitsvektor d, dem Druck p und der Dichte o gilt 
die Bernoullische Gleiehung, erstens in ihrer differentiellen Form: 


zweitens in ihrer integrierten Form: 
(dp 
+\ (2). 


In einem Gase, das sich gemäß der Adiabate 


p 


verhält, kann man in (2) 
dp “ 
\ o A) 
setzen. Da andererseits die Schallgesehwindigkeit 
p 
ist, kann man für (2) auch 
— — (6) 


schreiben. wo der Wert a, für die Schallgeschwindigkeit zum Geschwindigkeitsvektor vd, 
vehört. 


Die Kontinuitätsgleichung lautet: 


oder odivv+pgradoe=0 ....:.:. 
Nun ist 
do l 
erado: geradp= „grad 


Indem noch grad p aus (1) entnommen wird, erhält man schließlich für die Kontinuitäts- 
gleichung (7): 


v- grad 


4 T. E. Stanton: The variation of veloeity in the neiehbourhood of the throat of a eonstrietion in a wind 
channel. Teehnieal Report of the Aeronautical Research Commitee 1930-31, Vol. II, p. 728 (London 1932), Reports 
x Memoranda No. 1388, 

5), G..). Taylor: Recent work on the flow of eompressible fluids. Journal of London Mathem. Soe. Vol. >. 
p. 224, 1950, Deutsche Übersetzung in ZAMM, Bd. 10, S. 354, 1930 Vgl. Gl. (15) bzw. (30) ). 

“, Diese Strömung ist von H. Bateman, Proc. of the Nation. Acad. of Seiences, Vol. 16, p. 816, 1930, mit 
anderen mathematischen Hilfsmitteln als bei Taylor behandelt worden. 
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Man könnte nun den Geschwindigkeitsvektor d durch ein Potential ausdrücken, wir wollen 
statt dessen meist die Bedingung der Wirbelfreiheit 


unmittelbar hinschreiben, so daß wir in (8) und (9) die Grundgleiehungen der kompressiblen 
stationären Potentialströmungen haben. 


Die Differentialgleichung für das Geschwindigkeitspotential 9 würde für eine ebene 
Strömung in kartesischen Koordinaten lauten: 


ER 2 
Diese Differentialgleichung ist elliptisch, wenn 


04 - 2 2 


ist, also für Unterschallströmungen; sie ist hyperbolisch, wenn dv?’ > a” ist, also für Überschall- 
strömungen. 


x: 


120). 


Es sei kurz an das Kennzeichen hyperbolischer Differentialgleichungen erinnert. Kleine 
Störungen (die in unserem Fall etwa von Kratzern an den Wänden eines Strömungskanals 
ausgehen) pflanzen sich längs gegebener Bahnen, der Charakteristiken, in den Lösungsbereich 
hinein fort. Für die Überschallströmung heißen die Charakteristiken Machsche Linien und 
sind leicht nach der Töplerschen Sehlierenmethode beobachtbar. Bei elliptischen Differential- 
gleichungen (Unterschallströmungen) kann man dagegen von vornherein nichts über die Fort- 
pflanzungsbahnen von kleinen Störungen aussagen. 


2. Ebene Strömungen, die von einer Koordinate eines krummlinigen Koordinatennetzes ab- 
hängen. Wir wollen ebene Geschwindigkeitsfelder, die den Gl. ($) und (9) genügen, nach 
folgenden Richtlinien aufbauen. Es werde statt der kartesischen .-y-Koordinaten ein krumm- 
liniges &-n-Koordinatennetz eingeführt. Die Geschwindigkeitskomponente tangentiell zu den 
Kurven = const sei r,, die Normalkomponente sei Sowohl », als auch sollen nur von. 
abhängen, also auf den Linien 7 =const unveränderlich sein. Mit anderen Worten, es wird 
eine die Ebene oder ein Ebenenstück lückenlos und einfach erfüllende Schar von Kurven 
7 =const gesucht, auf denen die Tangential- und Normalkomponente der Strömungs- 
geschwindigkeit konstant bleibt. Ohne Beschränkung der Allgemeinheit können wir danach 
in der Ebene eine zweite Schar von Kurven £ == const annehmen, die zu den Kurven 7) = const 
orthogonal ist. 


Ein ähnliches Verfahren zum Aufsuchen von Strömungen hat G. Hamel’) bei zähen, 
inkompressiblen Flüssigkeiten mit großem Erfolg angewandt. Der Hamelsche Beweisgang 
benutzt aber wesentlich die zusätzliche Voraussetzung, daß das krummlinige Koordinatennetz 
quadratische Maschen im Kleinen bildet. Daher müssen für unsere Zwecke andere Beweis- 
mittel herangezogen werden. 


Wir haben die Gl. (8) und (9) auf das krummlinige orthogonale &7-Koordinatennetz um- 
zurechnen. Das Linienelement ds sei durch 


bestimmt. 9, und g, sind von Null verschieden und positiv. Erinnern wir uns daran, daß; 
die Geschwindigkeitskomponente in Richtung der &-Koordinate »,, in Richtung der -Koor- 
ddinate », sein soll, so ergibt sich: 


l 


7 G. Hamel: Spiralförmige Bewegungen zäher Flüssigkeiten. Jahresberieht der Deutschen Mathematiker- 
Vereinigung 1916, 85. 34. Den Stromlinien wird in dieser Arbeit die Bedingung auferlegt, daß sie mit einer Sehar von 
Koordinatenlinien zusammenfallen. 


| 


120 Follmien,. Zum Überganz von Unterschall- in Überschallströmungen 


wobei i, der Einheitsvektor in der Richtung, t, der in der „-Riehtung ist. Führen wir nun 
(die Forderung ein, daß », und ®, nur von ») abhängen sollen, so schreiben sich die Gl. (8) und (9): 


Or, v, 109, 109 


109; 
07 9, 05 ° 
Nach (6) ist 


Damit », und », nur von 7) abhängen, ist nach diesem Gleichungssystem hinreichend, 
daß die Ausdrücke 
109, 199, 


16 

9, 08 
"unktionen von », allein sind. Dann tritt nämlich in den Differentialgleichungen (13) und (14) 
S überhaupt nicht auf. Diese beiden Bedingungen sind auch notwendig; denn wir können 
diese Größen aus (15) und (14) ausrechnen und erhalten: 


7 + | | | | 
‚dm 


ıst von Null verschieden. Die rechten Seiten dieser Gleichungen sollen nur von 
abhängen, also ist die Bedindung für (16) auch notwendig. 


Aus «diesen Bedingungen folgt durch Quadraturen 
wo die p und 4 stetig differenzierbare, sonst aber bis auf weiteres willkürliche Funk- 
tionen sind. 


Außerdem müssen g, und 9, wegen ihrer geometrischen Bedeutung noch gewissen Be- 
ziehungen genügen. Es ist ja unsere Aufgabe, aus den Maßstäben 9, und g, das &-7-Kurven- 
netz selbst zu bestimmen. Wir gehen zweckmäßig dabei so vor, daß wir x, y als Funktionen 


von £,n ermitteln. Nach der Gleichung für das Linienelement (11) ist 


Da in dem Ausdruck für ds’ kein Glied mit d&d auftritt, also wegen der Orthogonalität 
der Kurven, muß 


| = 23 
sein. 
Wir fassen .r und y zu einer komplexen Koordinate zusammen: 
Dann ist dureh (21) und (22) \. und vegeben. Folglich ist 
02 ik; 
0: 
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wo k, und k, reelle Funktionen von £ und 7 sind. Mit diesen Ansätzen schreibt sieh die 


Bedingung (23) | | 
9, 9, (cos k, cos k, + sin k, sin k,) = 0 


oder 
cos (k,—k)=d. 


Daraus folgt, daß 
sin(k,—k)= + 


sein muß. Wir rechnen mit 

ei lvo N ei 8). 
Die Durchreehnung des Problems mit der zweiten in (27) vorgegebenen möglichen Form 
von eiks vollzieht sich natürlich genau so wie mit (28), liefert aber keine neuen Lösungen, 


so daß wir uns mit (25) begnügen können. 
Die Gl. (25) und (26) sind dann und nur dann miteinander verträglich, wenn die 


Integrabilitätsbedingung 


Dies ergibt zusammen mit (28): 
09, 
on Be 8% 
oder, nach Trennung in Real- und Imaginärteil, 
) k 
01 
k 109, 
on 9, 05 
Aus diesen Gleichungen läßt sich #, dann und nur dann ermitteln, wenn 
0 /1® 
on \y, OS, 


ist. Diese Beziehung zwischen 9, und 9, ist grundlegend. Ist sie erfüllt, so wird nach dem 
Beweisgang k, aus den Gl]. (30), (31), sodann z aus den Gl. (25), (26) bestimmt. 

Die Grundgleichung (32) für y, und 9, läßt eine unmittelbare differential-geometrische 
Deutung zu. Das Linienelement nach (11) könnte zunächst auch auf einer gekrümmten Fläche 
liegen, wobei die &- und n-Kurven Gaußsche Parameterkurven auf der Fläche wären. Damit 
die Fläche eine Ebene wird, was sie in unserer Anwendung sein soll, muß die Gaußsche 
Krümmung überall verschwinden. Die Gaußsche Krümmung läßt sich aber nach dem 
Theorema egregium von Gauß allein durch den Ausdruck für das Linienelement ausdrücken. 
Durch einfaches Einsetzen können wir bestätigen, daß unsere Grundgleichung (32) für q, und g, 
in der Tat mit der Bedingung identisch ist, daß die Ganßsche Krümmung überall verschwindet. 

In die geometrische Grundgleichung (32) führen wir die eingangs aus dem dynamischen 
Problem gefolgerten Ausdrücke (19), (20) für q, und q, ein: 


Da p, (5) wegen 9, >O von Null verschieden ist, verschwindet der Klammerausdrueck in (33) 


reziproke Wert des Klammerausdruckes, nämlich 


nur dann, wenn —=(0 ist. Schließen wir diesen Fall zunächst aus, so muß auch der 


von unabhängig sein. Dies ist auch hinreichend. Da (3) nicht (&) proportional 
) 
1 \> 


und Konstante sein, um (3b unabhängig von zu machen. 


Falls weder 4, so muß mit zwei reellen Konstanten e, und e, 


sein kann, müssen 


dd (m) (1) (35) 


sein, 


| 
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Die Bedingung (35) ist nur erfüllbar, wenn 


C, 
ist. Es wird also i 


ce, 


Nach (35) ıst 


wobei das Integral ein unbestimmtes ist. Setzen wir: 


er 


so wird 


dı, 
also sehließlich, mit Unterdrückung eines jetzt belanglosen Index bei der Konstanten 
AP (2) 


Damit haben wir das erste Lösungssystem gewonnen, das allen unseren Bedingungen genügt. 
Nehmen wir q,, das bisher von Null verschieden angenommen wurde, jetzt gleich O0 an, 
so muß man, um den Ausdruck (34) von »; unabhängig zu machen, mit einer Konstanten e 
setzen, also 
nehmen, wonach dann 


wird. Wenn wir noch 
= Pi), 


setzen, was weiterhin zweckmäßig ist, erhalten wir damit als zweites Lösungssysten 


dPi(2) 
Te 


(il). 


Die Annahme, daß q,=0 ist, während dann q, wegen 9,=-0 von Null verschieden sein muß, 
liefert eine Lösung, die bereits im ersten Lösungssystem enthalten ist. 
Bisher hatten wir stets den Fall ausgeschlossen, daß 
dq, 
dn (4) 
ist. Die Bedingung (33) ist aber auch in diesem Falle befriedigt, so daß wir nach (19) und 
(20) mit 


(das dritte Lösungssyvstem 
dr (2) 


erhalten. Die Konstante ist dabei in die ohnehin willkürliche Funktion 4, (1) einbezogen. 
2 
Durch die drei Lösungssysteme (1), (ID), sind sämtliche Möglichkeiten, und 9, 
entsprechend den Bedingungen (19), (20) und (33) zu wählen, erschöpft. 
Aus g, und 4, berechnen wir k, nach den Gl. (30), 81) und z nach den Gl. (25), (26). 
Nach dem ersten Lösungssystem (I) erhalten wir 


d d 
Ok, dE 1 dn 


5 
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Daraus folgt bis auf eine belanglose Integrationskonstante 
| | 


Nach (25), (26) wird 


02 


on dı 
Setzt man In P(&)= P*(£), In (n) = (17), so Ist 


()ll—iele 


de 
Bis auf eine unwesentliche Integrationskonstante ist wegen are ; danach 
— 
5 
i] 
re (1,2). 
Drücken wir z durch Polarkoordinaten aus, also 2=re'”, so wird 
Durch Auflösung nach und P*(£) wird 
1 2 
elnr + Inr -ced=(l+e) . . (13). 


Demnach sind die Kurven 7=const logarithmische Spiralen. Der Tangens des Winkels 
zwischen hadiusvektor und Tangente der Spirale ist 2 Die Kurven £= eonst sind eben- 


falls logarithmische Spiralen. 

Unabhängig von der hier dargelegten allgemeinen Theorie hat L. Prandtl bereits im 
Jahre 1955 die Vermutung geäußert, daß die logarithmischen Spiralen Linien konstanter 
Tangential- und Normalkomponente einer kompressiblen Potentialströmung sein könnten. 
Diese Vermutung war durch die schon bekannte Tatsache nahegelegt, daß die Eigenschaft 
jedenfalls für zwei Grenzfälle der logarithmischen Spiralen gilt, für die konzentrischen Kreise 
nach der eingangs erwähnten Lösung von G. J. Taylor, für das Strahlenbüschel nach der 
von Th. Meyer‘) konstruierten Überschallströmung um eine Eeke. Der Beweis für die 
Prandtlsche Vermutung folgte aus unserer allgemeinen Theorie. | 

Nach dem zweiten Lösungssystem (II) wird mittels (30), (31): 


ok, __ Ok, 
On 
Man erhält also 


Nach (25), (26) ist 


e-ier® do (N) e-ieP 


07 
Demnach wird bis auf eine unwesentliche Integrationskonstante 


Mit Polarkoordinaten erhält man: 


Th. Meyer, |]. e. Abschnitt A. 
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Die Kurven 7 = const sind demnach konzentrische Kreise. Diese Möglichkeit ist von G.J. Taylor 
benutzt worden, wie in der Einleitung bemerkt wurde. 


Das dritte Lösungssystem (III) liefert nach (30), (31) die Gleichungen für k;: 


ok, ok, 
Daher wird 


Wo 


ist. Nach (25), (26) wird 


- 


Daher wird 


Bei y, VO erhält man durch Einführung von Polarkoordinaten 9 = const bei 1) = const, also 
des Strahlenbüschel, das sehon in der Arbeit von Th. Mever benutzt ist. Im allgemeinen 
"all zerlegen wir in (111,2) das Integral in Real- und Imaginärteil: 


3,4), und 9, sind durch die Beziehung 
do, 
tg 10, 
d 7 


verknüpft, die durch Differentiation von (1113). folgt. Wir zerlegen (II1,2) in Real- und 
Imaginärteil: 


Dieses System stellt die Gleichungen der Kurven ») = const in Parameterform dar, als Parameter 
ist Sanzusehen. Durch Elimination von P(£) bekommen wir 


Die Kurven »)  const sind demnach Geraden, die mit der Neigung tg 0 durch die Punkte 
0,07), 9= 9, gehen. Denken wir uns durch die Punkte (1). 49=%,(n) eine 


Kurve gelegt. so ist der Riehtungstangens dieser Kurve was nach (IIl,4) gleich — 
% 


ist. Demnach steht die Geradenschar (111,6) auf der Kurve 4, (n). 9 9, (1) senkrecht. 
Wegen der Willkür in 0, und 0, ist das keine besondere Beschränkung der Allgemeinheit 
der Geradensehar. Natürlich darf die Schar der Geraden in dem Ebenenstück, wo sie als 
Kurven benutzt werden, sich nicht überschneiden. 


Auch IH. Bateman hat diese Eigenschaft von Geradenscharen im 3. Absehnitt der 
angeführten Abhandlung nachgewiesen und dargetan, daß die Strömung dann nur eine Über- 
schallströmung sein kann. Solche Geradenscharen haben im übrigen schon längst Verwendung 
zefunden als Träger konstanter Tangential- und Normalgeschwindigkeit. Sie treten auf bei 
Überschallströmung um Tragflügel in unendlich ausgedehnten Gasen®). Die Geraden sind 
dann Machsche Linien oder, mathematisch gesprochen, Charakteristiken der Differentialgleichung. 


Wir fassen zusammen. Als mögliche Kurven konstanter Tangential- und Normal- 
geschwindigkeit einer kompressiblen Potentialströmung stellten wir fest, erstens logarithmische 
Spiralen mit den beiden Grenzfällen der konzentrischen Kreise und des Strahlenbüschels, 
zweitens beliebige Geradenscharen in solchen Gebieten der Ebene, die sie einfach überdecken. 
Wesentlich neue Strömungen lassen sieh nur bei Benutzung logarithmischer Spiralen erwarten. 

"», Vel. etwa J. Ackeret: Beitrag „Gasdynamik“* im Handbuch der Physik (Geiger-Scheel) 7. Bd., Berlin 1927, 


Abb. 18, 19. X. Busemann: Beitrag „Gasdynamik* im Handbuch der Experimentalphysik (Wien-Harms) 4. Bd., 
I. Teil, Leipzig 1931, Abb. 54, 64. 
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3. Logarithmische Spiralen als Koordinatenlinien. Als krummlinige Koordinatenlinien 
ihren wir die logarithmischen Spiralen 
E=alnr „=iIinr+tad . . . (88) 
»in. Die Bezeichnung ist gegen (1,3), wo die logarithmischen Spiralen als tauglich für unseren 
/ week erkannt wurden, etwas geändert. «a ist der Tangens des Winkels zwischen der Tan- 
vente an die Spirale £=const und dem Radiusvektor. Die zur Umreehnung der Differential- 
oleichungen der Strömung notwendigen Größen 9, und 9, sind beide gleich: 


(39. 


Die Geschwindigkeit werde auf die Schallgeschwindigkeit bei Ruhe, «a, bei vd, 0, be- 
zoeen, die man auch als Schallgeschwindigkeit im Kessel bezeichnet, wenn man an die 
experimentelle Verwirklichung denkt. Es werde also 


oesetzt. Die Gl. (6) geht demnach über in 


Nach (5) und (3) folgt daraus 


wenn p, der Druck im Kessel (bei der Geschwindigkeit d, 0) ist. 

Es empfiehlt sieh nieht, die Geschwindigkeitskomponenten durch ein Potential aus: 
zudrücken, für das man eine nichtlineare Differentialgleichung 2. Ordnung bekommen würde. 
Vielmehr führen wir den Betrag von w und den Winkel der Geschwindigkeit gegen die 
Linie 7 = const ein. Wir setzen also die auf a, bezogenen &Geschwindigkeitskomponenten | 


w,=—=Weos9, .... 0... 


N und © sind nach unseren Ausgangsvoraussetzungen Funktionen von 7 allein. 


Die Kontinuitätsgleichung (5) schreibt sich dann: 


Die Bedingung der Wirbelfreiheit (9) wird 
asınJ) cos daW do 


Da in beiden Differentialgleichungen (44), (45) die unabhängige Variable », nicht explizit auf- 
tritt, so läßt sich die Ordnung des Differentialgleichungssystems herabsetzen. © wird als 
Funktion von W betrachtet: 


d; dWdn 
/ wird aus (45) ermittelt: 
7 
dW MW asın cos an 
44). 
AN 
In (44) eingesetzt, ergibt dies die Differentialgleichung für (WW): 
do 7 cos a sin? |. 


_ 
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Wie aus dieser Differentialgleichung O(\W) bestimmt wird, soll später auseinandergesetzt 
werden. Ist erst (WW) bekannt, so sind nur noch Quadraturen zu erledigen. Nach (47) und 
(45) erhält man nämlich: 

dn a@+1 (@+1l) 

r + . . . . . . . . . (49) 
dW x 
1— —;—- 

so daß bei bekanntem (WW) durch Quadratur 7 (W) bestimmt werden kann. Dann sind 
und bekannt, also auch die Geschwindigkeitskomponenten und als Funktion 
der Ortskoordinate 7. 


Die Differentialgleichung der Stromlinien ist 


v, 

5 
Bei bekanntem 9d,) läßt sieh für die Stromlinien (67) nach (50) durch Quadratur bestimmen, 
Die im UÜberschallgebiet auftretenden Machschen Wellen bilden mit den Stromlinien den 
„Machschen Winkel“ y, dessen Sinus gleich dem Verhältnis der lokalen Schallgeschwindigkeit 


5 


um Geschwindigkeitsbetrag a, W ist. Der Tangens dieses Winkels ist demnach 


- 


y ist nach dem Vorhergehenden als Funktion von », bekannt. Die Differentialgleichung der 
Machschen Wellen ist unter Benutzung von (50) und (51) 


= (17) läßt sich nach dieser Gleichung für die beiden Scharen der Machschen Wellen durch 
Quadratur bestimmen. & und ,, sind dann nur noch dureh die ursprünglichen kartesischen 
Koordinaten auszudrücken, um Stromlinien und Machsche Wellen zeichnen zu können. Der 
Druckverlauf läßt sich aus (42) berechnen. 


max 


4 


—- N 


Abb. 2. P(W)für@e=1. Die Nummern der Integralkurven 
sind inder Figur in Kreise gesetzt, während sie im Text und 
in den Tabellen dureh fetten Druck hervorgehoben sind. 


Damit ist ein vollständiger Rechnungsplan aufgestellt bis auf die Integration der 
Differentialgleichung (48), der wir uns jetzt zuwenden. 


Zunächst stellen wir fest, daß wegen der Periodizität der rechten Seite von (48) die 
Differentialgleiehung nur im Bereich O< © 7 diskutiert zu werden braucht. Sobald eine 
Integralkurve 9-0 oder = 7 erreicht, setzt man sie bei demselben W und einem um + 7 


_ | 
/ 
ig) 
w’—1 
| | | | | | 
| | | 7 | | | 
| 3) | | 
. N - + + + + + + + 
2a 
+ + + + + + + + + — 
| 
3 | | "eo | | 
m" | | 
08 12 1 20° | 
| 
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versetzten © fort (vgl. etwa die Integralkurve 1 und 3 im O(W)-Bild, Abb. 2). Auch W ist 
beschränkt, da wegen (42) 
H / r -. 
sein muß. W,,.. gibt gerade das Verhältnis der Ausströmungsgeschwindigkeit ins Vakuum 
zur Schallgeschwindigkeit a, im Kessel an. Die Differentialgleichung (45) braucht demnach 
nur in dem rechteckigen Bereich 
x 
betrachtet werden. | 
Der Durehgang durch die Schallgeschwindigkeit erfolgt nach (41) bei 


a\? 


oder. wenn wir das betreffende W mit einem Stern versehen, bei 


/ 


Diese Scheidelinie zwischen Unterschall- und UÜberschallgebiet ist in das 9 (W)-Bild (Abb. 2) 
eingetragen. | 
Der einzige singuläre Punkt der Differentialgleichung (48) in dem bezeichneten Bereich 
do 0 do 


ct W— W ie unbestimmte Form -. annımmt. Da ferner nur 
ist W=W W di Pur am 


Rande W = W,,., unendlich werden kann, laufen die Integralkurven 9(W) durch den ganzen 

Bereich von W. Um eine doppelte Überdeekung der &-7;-Ebene zu vermeiden, muß deshalb n, 

das nach (49) zu bestimmen ist, eine monotone Funktion von W sein. Nun ist 7 5 ==) bei 
W=W,, wo 

Kommt man aus dem Unterschallgebiet, so kann man die Strömung nicht nach größeren W 

hin fortsetzen, als nach (55) bestimmt sind, da man sonst keine eindeutige Strömung mehr 

bekommen würde. W, bestimmt daher eine Grenzlinie der Strömung. Nach (43) und (41) ist 

auf der Grenzlinie W, nach (55) gerade die Geschwindigkeitskomponente senkrecht zu den 

Spiralen ») = const gleich der örtlichen Schallgeschwindigkeit. Die Kurve W, berührt die 


Schallgrenze W* bei © = und erreicht W,,., bei 9 ==. Die physikalische Bedeutung der 


Grenze W, wird später ausführlich erörtert werden. 
Die Integralkurven (ll) sollen für a=1, z= 1,405 konstruiert werden. Die logarith- 


mischen Spiralen des &-17-Koordinatensystems bilden dann mit den Strahlen durch den 
Ursprung Winkel von 4° Zur graphischen Lösung der Differentialgleichung (48) bedienen 
wir uns der Isoklinenmethode, deren Ergebnisse wir nach dem Picardschen Iterations- 
verfahren mittels graphischer Quadraturen schrittweise verbessern. Nur in der Nähe des 
singulären Punktes WW... 9=7 versagt dies Verfahren. 

Zur Diskussion der Integrale bei diesem singulären Punkt setzen wir 


Dann wird durch Entwiekelung nach den als klein angenommenen Größen W, und ©, 
d W, 1 . . . . . . . . . . . 


Dabei sind auf der rechten Seite Glieder weggelassen, die mindestens quadratisch in W, und ©, 
sind. Nach allgemeinen Kriterien ist der singuläre Punkt W, 0, ©, -0 ein Sattelpunkt. 
/,wei Integrale gehen durch den singulären Punkt, von denen das eine W, 0 oder WW 
beliebigem ©, also für uns unbrauehbar ist, während das andere mit 


+ 


MAN bei 


oder 


I\ı/z —1 
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beginnt. Indem dieser Ausdruck für © in die rechte Seite von (48) eingesetzt wird, gelingt 
es, durch Quadratur eine zweite Näherung für (W) formelmäßig zu finden'®). Mit dieser 
rechnerischen Näherung kommt man schon genügend weit vom singulären Punkt ab, um die 
weiteren Verbesserungen wie bei den anderen Integralen graphisch ausführen zu können, 


Aus (49) folgt unter Benutzung von (57), daß 7 sich in der Nähe von W,,., wie 
n=V -W) tn » 2 2 2020020. (88) 


verhält, 7, ist dabei eine Integrationskonstante, für die in unserem Zahlenbeispiel 1,345 
gewählt ist. Eine andere Wahl von „7, würde nur eine Drehung des gesamten Strömungs- 
bildes bewirken, also physikalisch belanglos sein. Durch Kombination von (57) und (58) 
ergibt sich 


(-) == (7 — Nu) . . . . . . (59) 


in der Nähe von W=W a, I(=7,9=%- 


Obwohl im äußersten Unterschallgebiet bei W=0 keine Singularität auftritt, sind einige 
Bemerkungen über das Verhalten der Integralkurven in dieser Gegend geboten. Die einzige 
Integralkurve, die durch 90, W=0 geht, ist W 0. Im übrigen ist nach (48) bei sehr 
kleinem W 


Man vergleiche dazu das letzte im Unterschallgebiet gezeichnete Stück der Integralkurve 1. Mit 
derselben Näherung folgt aus (49) 


und aus (60) und (61) 


Das Strömungspotential y, das zu der entsprechenden Strömung gehört, berechnet sich aus 


109 
7 
Ya 
zu 


wo CU und e zwei Konstante sind. Nach (63) ist p das Potential einer inkompressiblen 
Strömung, wie man von einer Strömung im äußersten Unterschallgebiet nicht anders er- 
warten kann. 


Damit ist die Diskussion der Integralkurven 9(W) abgeschlossen. Die Ermittelung von 


(U) nach (49) vollzieht sich durch Quadratur ohne Besonderheiten, abgesehen von der 


Zeichenumkehr von I bei W-W,, die wir schon erwähnten. Die bei dieser Quadratur 


auftretende Integrationskonstante ist nach einer früheren Bemerkung belanglos. In der Tabelle 
sind die Ergebnisse für die drei Integralkurven 9(W) der Abb. 2 mit den zugehörigen ;, 
wiedergegeben. Die Lösung 1 geht durch den singulären Punkt W= Wax = 2222, O=n, 
die Lösung 2 findet ihre Fortsetzung jenseits W,, das hier gleich 1,187 ist, in a. Von der 
Lösung 3 ist nur ein Stück bis zu W,, das in diesem Falle 1,124 ist, gezeichnet und tabuliert, 
da uns reine Überschallströmungen, wie sie jenseits W, existieren, weniger interessieren. Die 
Werte von © und 7 beim Durchgang durch die Schallgeschwindigkeit W* — 0,912 sind in die 
Tabelle eingetragen. Für die Fortsetzung der $(W') über O0 oder 7 hinaus ist die am Eingang 
dieses Absehnittes gemachte Bemerkung ausgenutzt worden. 


4. Diskussion der Strömungen im Koordinatennetz der logarithmischen Spiralen mit a1. 
Bei der Lösung 1 fällt die Grenze W, nach (55) mit der natürlichen Grenze jeder Überschall- 
strömung |, nach (583) zusammen. Wir beginnen daher mit der Diskussion der zur Lösung 1 
gehörenden Strömung. 


10), Vgl. O0. Perron, Math. Annalen Bd. 75, S. 256, 1914, wo Differentialgleiehungen des Typus 4 (x) Ur) 


in der Umgebung einer Nullstelle von g und f erörtert werden. 


| 
f 
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Die Abb. 3 zeigt das Stromlinienbild mitsamt einigen Machschen Wellen in der Nähe 
les Durchgangs durch die Schallgesehwindigkeit. Die logarithmische Spirale, wo die Schall- 
seschwindindigkeit erreicht wird, ist in dem benötigten Umfange gestrichelt eingezeichnet. 
Die Machschen Wellen stehen auf den Stromlinien an der Schallgrenze senkrecht. Der 
Abstand zwischen zwei Stromlinien verringert sich, wenn man aus dem Unterschallgebiet 
sich der Schallgrenze nähert; oberhalb der Schallgrenze tritt eine Erweiterung der Stromfäden 
ein. Der Durchgang durch die Schallgeschwindigkeit vollzieht sich entgegen manchen Ver- 
mutungen ohne weitere Besonderheiten der Stromlinien. Man kann sich zwei der gezeichneten 
Stromlinien verfestigt denken und erhielte dann die Strömung durch einen gekrümmten Kanal. 


Abb. 3. Stromlinien und Machsche Wellen der Strömung I beim Durchgang dureh die Schallgeschwindigkeit. 


Tabellen zur Strömunefl. 


2.222 | 312 1.345 85 0.999 
2,2 3139 1.365 8.0 0,998 
2,1 3,126 1.454 7.0 0.992 
2,0 3.111 1.548 6,0 0.981 
1.9 3.097 1.649 5.0 0.936 
1.8 3.076 1,753 4.6 0.902 
1.7 3,055 1.865 4.2 | 0.852 
1.6 3.033 1.980 3.8 | 0,780 
1.5 | 3.010 | 2.098 3,4 | 0.675 
1,4 2988| 2,228 3.15 | 0.588 
1.3 2,956 2.357 2,087 | 0.527 
1.2 2,918 2.5083 2,8 0.442. 
1,1 2,878 | 2.656 2.6 0,349 
1.0 2830 2,825 24 | 0.254 
0,912 2,781 | 2.987 2,3 | 0.206 
0.9 2,776 3.005 2,2 | 0,161 
0.8 2708 3.230 2,1 | 0,121 
0,7 2,6300 3.472 2,0 | 
0.6 2,531 | 3,748 1.9 0.0558 
0,5 2398 A,0SS 1.85 0.0439 
0.4 220 4,518 1.8 | 0.0331 
0.3 1.971 5.035 1.75 0.0242 
0.2 1.608 5.796 1.72 0.0200 
0.1 1.017 6,972 1.345 | 0 
0.05 0.321 8.367 | 
0.036 9,009 
0,020 — 0,584 10,177 
0.010 | 11,573 
0,005 -- 1,973 12.955 
0,00156 — 3,142 15.293 | 


| 


Abb. 4 
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getragen, Das Strömungspotential im äußersten Unterschallgebiet war in Gl. (63) gegeben. 


Abb. 5 gibt das Überschallgebiet wieder. Bei der Ausströmung ins Vakuum, wo W=W,, 
ist, wird 97, so daß die Strömungsrichtung mit der zu W,,.. gehörenden logarithmischen 
Spirale zusammenfällt. Diese ist in unserem Zahlenbeispiel durch 7, = 1,345 gegeben. (Vgl. die 
an Gl. (58) anknüpfenden Erörterungen.) Diese Spirale ist in Abb. 5 gestrichelt eingetragen. 


Die benachbarten Stromlinien, wo also 7 — 7, klein ist, sind nach den Gl. (50) und (59) in 
erster Näherung durch 


2x 2 —1 


sezeben. Die Integrationskonstante hängt von der Wahl der Stromlinie ab. Fügt man die 


Abb. 4 und 5 zusammen, so sieht man, daß die Ebene mehrfach überdeckt wird, entsprechend 
der Tatsache, daß sich „7 um mehr als 27 geändert hat. (Vgl. Tabelle). Man kann demnach 
nur ein Teilgebiet der berechneten Strömung verwirklichen. 

Der Druck p. bezogen auf den Kesseldruck p,. ist in der Tabelle als Funktion von „7 aul- 
geschrieben. Ferner ist in den Abb. 6, 7 und 5 der Druck senkrecht zu einer Stromlinie 
aufgetragen. Es ist die Stromlinie gewählt, die in den Abb. 5, 4 und 5 die zweitnächste am 
Ursprung des angedeuteten Koordinatensystems ist. Abb. 6 gibt den Druckverlauf beim Durch. 


gang dureh die Schallgeschwindigkeit:; die Schallgrenze (2 0,327) ist durch Strichelung her- 
0 


-50 


/ 


Stromlinien der Strömung im Unterschallgebiet. 


Abb. 5. Stromlinien der Strömung # im Überschallgebiet. 
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In Abb. 4 ist ein größeres Gebiet der Unterschallströmung in kleinerem Maßstabe auf- 
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Abb. 9. Stromlinien und Machsche Wellen Abb. 10. Druckverlauf der Strömung 2. 
der Strömung 2. 


Tabellen zur Strömung 2. 


) 

0.2 1,309 5.300 3.62 0.834 
0.3 1.778 4.501 8,5 
0.4 2,009 3.996 3,4 0,784 
0. 2,177 3,611 3. 
0,6 2.284 3.399 3.2 0,750 
0.7 2,359 3.099 3.1 0,699 
2.102 2.8090 3.0 0,666 
0,0 2.488 2.739 2,9 0.624 
0.012 2.439 2,720 2,8 0,573 
1.0 2,439 2,616 2.12 0,527 
1,3 2,402 2.551 2,69 0,482 
1.15 2.379 2.534 2.60 0,440 
1.187 2.348 2.530 58 0.416 

2,55 0,374 
2.53 0.312 


Abb. 11. Stromlinien der Strömung 2a. 


vorgehoben. Das Gebiet um die Schallgeschwindigkeit, das in Abb. 6 wiedergegeben wird, 
ist größer als das in dem Stromlinienbild 3 gezeichnete, wie man an der Bezifferung der 
Koordinatenachsen erkennt. Abb. 7 zeigt den Druckverlauf im Unterschallgebiet entsprechend 
den Stromlinien in Abb. 4. Abb. 8 bezieht sich auf den Druckverlauf im Überschallgebiet 
und ist den Stromlinien in Abb. 5 zugeordnet. 


Die Lösung 2% weist eine Grenzlinie der Strömung bei W,= 1,187 <W,,, auf. Abb. 9 
zeigt vier Stromlinien nebst einigen Machschen Wellen. Die Schallgrenze ist wiederum 


\\ \\9972 
\ \ 
y \ % 
an 
[7 
8] Go (A6760) 
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gestrichelt. Zu W, 1,157 gehört eine Spirale, die zahlenmäßig durch n,,;, = 2,530 definiert 
ist. Die vier Stromlinien der Abb. 9 nach der Lösung 2% sind durch solche nach der Lösung 2a 
in Abb. 11 fortgesetzt. Die Stromlinien haben also bei dieser Art Fortsetzung für 7 = 2,530 
eine Rückkehrkante, während im 9 (W)-Bild bei W, keine Singularität auftritt. Die Strömung 
nach Lösung 2a und Abb. 11 ist nicht besonders interessant, da sie sehr bald in eine radiale 
Ausströmung ins Vakuum übergeht. Es ist leicht nachzurechnen, daß die Stromlinien an der 
Grenzlinie unendliehe Krümmung haben. In Abb. 10 ist der Druckverlauf der Strömung von 
ldy 


bei 


Abb. 9 senkrecht über der untersten Stromlinie aufgetragen. Das Druckgefälle 


der Spirale 7),,;, unendlich, wie sofort aus 7 =0 folgt. Den größten Aufschluß über die 


eigentümliche Begrenzung der Strömung 2% durch die Spirale 7 = ,,;„n gewinnen wir aus der 
Betrachtung der Machschen Wellen. Es ist ohne weiteres festzustellen, daß die Grenz- 
spirale 7, selbst eine Machsche Welle ist. |, war ja nach (55) dadurch bestimmt, daß 
die Geschwindigkeitskomponente senkrecht zu der Spirale ,,;. gleich der örtlichen Schall- 
geschwindigkeit ist. Diese Spirale muß daher nach der Definition des Machschen Winkels y 


Tabelle zur Strömung 2a. Tabelle zur Strömung 3. 
1.187 2.348 2.53 0,2 2.283 
1.2 2.332 2.556 0.3 
1.3 2.500 0.4 2942 
1.4 2.031 2,758 0,5 3.160 
1.5 1.708 3.074 0,6 3.345 
1.6 1.541 3.558 0,7 3514 
1.7 1.303 1.117 3.660 
1.5 1.111 4,833 0,0 3.792 
1.0 5,624 1.012 3.804 
2,0 6.615 1.0 3.902 
1.1 3.995 
1.124 4.015 


vel. (S1)) mit einer Machschen Welle zusammenfallen. Die Schar der Machschen Wellen, 
zu der die Grenzspirale gehört. ıst durch 


d 7 1 + Ciy () ig 


bestimmt. Dabei ist von den beiden Differentialgleichungen (52) für die beiden Scharen der 


\Machschen Wellen natürlich die ausgewählt, die für „= 7— Obein= n,in unendlich macht, 


weil dann », = const = 7,,,;, zu der Schar gehören muß. Die Grenzspirale ist aber nicht nur selbst 
eine Machsche Welle, sondern eine Enveloppe ihrer Wellenschar. Freilich kann man die 
Wellensehar nieht explizit darstellen, aber man kann die ausgesprochene Tatsache dadurch 
zeigen, daß man die Differentialgleichung (65) in der Nähe von 7 = ],„in genähert löst. Die 
rechte Seite von (65) ist eine Funktion von „, allein, die sich bei ,,;n bis auf konstante oder 
0,426 
verschwindende Bestandteile wie ' verhält, wie eine elementare Entwickelung des 
in 
Zählers und Nenners der rechten Seite von (65) zunächst nach W— W,,, das schließlich durch 
7 N, ersetzt wird, nachweisen kann. Die Machschen Wellen dieser Schar sind demnach 


in der Nähe von „=, =2,583 durch 


bestimmt, wo Sp der Scharparameter ist, der gerade den Wert von £ auf der individuellen 
Machschen Welle im Berührungspunkt mit der Einhüllenden 7 = 2,53 ist. 
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5. Ebene Strömung mit Rotationssymmetrie. Wir wollen noch an einem anderen Beispiel 
die Verhältnisse an der Grenzlinie studieren. Die Konstante a in der Formel (35) für das 
Koordinatennetz der logarithmischen Spiralen werde gleich O0 gewählt. Die logarithmischen 
Spiralen 7 = const gehen in konzentrische Kreise über, während die Koordinatenlinien & = const 
in radiale Strahlen ausarten. Man hat also ein Polarkoordinatennetz '') (r, 9) und die radialen 
und tangentialen Geschwindigkeitskomponenten sind auf konzentrischen Kreisen konstant. 
Dies ist gerade das von Taylor behandelte Beispiel. Im Unterschied zu der Taylorschen 
Untersuchung wird hier eine geschlossene Lösung des Problems gegeben, die sich einfach 
dureh Spezialisierung der Theorie für «=-0 ergibt. Dadurch gelingt es, die Diskussion der 
Strömung wesentlich weiterzuführen. 


Bezeichnen wir die radiale und die zirkulare Geschwindigkeitskomponente mit r, und 
"y, so setzen wir nach (43): 


v,=a,Wsin®, . 


Der Zusammenhang zwischen W und © ist nach (48): 


do cos sın 9) 


(' ist die Integrationskonstante. Aus Gl. (49) könnte man sodann In r =, als Funktion von W 
berechnen: 


1— w°’ 


sin’ ist aus (69) zu entnehmen. 


Etwas einfacher gelangt man zum Ziele, wenn man berücksichtigt, das r, und », aus 
einem Potential 9 nach der Formel 


bestimmt werden. Da ”, nur von r abhängen sollen, ist 


= also vg: 
wo A eine Konstante ist. Demnach wird 


—- 


Normieren wir den radialen Abstand r so, daß die tangentielle Geschwindigkeitskomponente 
für r=1 gleich der Schallgeschwindigkeit im Kessel'?) wird, so wird also 


1 
und mit Benutzung von (69) 
\ ] 
l 
l 2 War 
— -,) 
11) Genan genommen wird für nach <= — n=Inr, so daß gegenüber dem üblichen Polarkoordinaten- 


system der Durehlaufungssinn der Winkelkoordinate umgedreht ist. Zur Bequemlichkeit des Lesers ist im folgenden 
gleich der Übergang zu dem üblichen Polarkoordinatennetz vollzogen, was nur geringfügige Änderungen an den 
ursprünglichen Formeln bedingte. 

1?) Diese Normierung ist nur für rein radiale Strömung unmöglich, die uns niehts angeht. 


| 
Diese Differentialgleiehung läßt sich durch Trennung der Veränderlichen sofort integrieren: 
f 
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Durch (69%) und (70) ist der Zusammenhang von r, W und © in geschlossener Form gegeben, 
also die Strömung vollständig bestimmt. Auf etwas anderem Wege gelangte H. Bateman 
(l. e.) zu einer geschlossenen Darstellung dieser Strömung. 


In der Abb. 12 ist nach Formel (72) W über = für die fünf Werte der Integrations- 


konstante (0:0,3: 1:2; ©) aufgetragen. Für =0 erhält man den Potentialwirbel, für 
die rein radiale Strömung. Die Stelle, wo die Strömungsgeschwindigkeit gleich der lokalen 
Schallgeschwindigkeit wird, ist durch die Linie W* gegeben; die maximale Geschwindigkeit, 
welche der Ausströmung ins Vakuum entspricht, ıst mit W,,, bezeichnet. Man erkennt aus 
der Abbildung, daß zu jedem Wert von r im allgemeinen zwei Werte von W gehören. Die 
Grenzlinie, bis zu der man die Strömung aus dem Unterschallgebiet kommend fortführen kann, 


ist ın der Abb. 12 mit bezeichnet und durch bestimmt. Würde man die Strömung 


formal nach den Formeln (69) und (72) über die Grenzlinie W, hinaus fortsetzen, so würde 


man doppelte Überdeekung der Strömungsebene, also eine Rückkehrkante der Stromlinien 
bekommen. 


Die Bedinzung IW 0 für die Grenzlinie W, ist nach der Differentialgleichung (70) dann 


befriedigt, wenn 
»—1l_.-. 
Wie wir in der Formel (55) schon für beliebiges a feststellten, ist danach die Geschwindig- 
keitskomponente normal zur Grenzlinie gleich der lokalen Schallgeschwindigkeit. Auf diese 
Tatsache hat bei dem vorliegenden Beispiel schon G. J. Taylor auf Grund einer Bemerkung 
von W. @. Biekley hingewiesen. Unter Benutzung von (72) gilt diesmal weiter auf der 
Grenzlinie, welche allgemein den Index 0 zugeordnet erhält: 


Danach ist eine Hyperbel, welche die Sehallgrenze W bei berührt und die maxi- 


male Geschwindigkeit W,,., bei —: erreicht. 


In der Nähe von W=0, also im äußersten Unterschallgebiet ist nach (72) angenähert 


une nach (60) 


| | | | 
15 
| | | 
| | 
| | | | | 
| | | | 
$ 05 19 15 20 25 
Abb. 2. 
») 
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Daraus ergibt sich in erster Näherung für die Stromfunktion 
v=-Cd—Inr 
und für das Geschwindigkeitspotential 
y=9+Chr. 
Daß dies das Potential einer volumbeständigen Strömung ist, war zu erwarten. 
Für wird 
also der Potentialwirbel erhalten. Für =» wird tg © auch unendlich, = 1, cos) 


also wird die radiale Strömung erhalten, die jeweils wegen des Zı ısammenfallens von W, mit W u 
als reine Unterschallströmung oder als Überschallströmung realisierbar ist. 


Abb. 13. Stromlinien und Machsche Wellen der Strömung (€ 0,5) 
beim Durehgang dureh die Schallgeschwindigkeit. 


Die Stromlinien wurden für nach der Formel r er = berechnet. 
‘ dr 
Abb. 13 zeigt fünf Stromlinien, die aus dem Unterschallgebiet kommen. Die. Sch: llgrenze 
ist wieder gestrichelt, im Überschallgebiet sind einige Machsche Wellen eingezeichnet. Die 
Machschen Wellen sind nach (51) und (52) und (71) dureh die Formeln: 


d etg (-) tg 
unter Benutzung 


bestimmt. Die Grenzlinie liegt bei r, 1,007: W, = 1,283. Denkt man sich diese Stromlinien 
der Formeln (69) und (72) oder der Abb. 13 formal weiter ins Überse "hallgebiet hinein fort- 
gesetzt, so bekommt man das Stromlinienbild 15, das sich wenig von einer rein radialen 
Expansion ins Vakuum unterscheidet. 


Abb. 14. Druckverlauf der Strömung von Abb. 13. Abb. 15. Reine Überschall- 
strömung bei 


| 0507 | 
024 \\ | 
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In Abb. 14 ist der Druckverlauf, der zum Stromlinienbild 13 gehört, über der mittleren 
dr 
dw. 
sofort folet. Aus dem Gleichgewicht des Druckgefälles normal zu den Stromlinien und der 
Zentrifugalkraft kann man «dann auch leieht erschließen, daß die Krümmung der Stromlinien 
an der Grenzlinie unendlieh ist. 


Stromlinie aufgetragen. Das Druckgefälle ist an der Grenzlinie unendlich, wie aus 0 
( 


Die Grenzlinie r = r, ist eine Machsche Welle, da nach (55) die Grenzlinie dadurch 
bestimmt werden kann, daß die radiale Geschwindigkeitskomponente gleich der lokalen Schall- 
geschwindigkeit ist. Der Grenzkreis gehört zu der Schar von Machschen Wellen, welche 
der Differentialgleiehung 


genügen, day © für r=r, wird. Eine Lösung dieser Differentialgleichung ist dann r  r,. 
wie es verlangt war. Es läßt sich auch hier nachweisen, daß die Grenzlinie eine Enveloppe 
der Machschen Wellen ihrer Schar ist. Ohne die Wellenschar explizit zu kennen, kann 
man diese Tatsache durch näherungsweise Lösung der Differentialgleichung (74a) bei W=W,, 
r 7, feststellen. Die rechte Seite von (74a) wird in der Nähe von W = W, unter Benutzung 
von (69) nach WW, entwickelt, das dann mittels (72) nach r — r, entwickelt wird. Man 
vyr—r, 
Wellen dieser Schar sind demnach in der Nähe von r = r, 1,007 durch 


erhält dann für die rechte Seite von (74a) näherungsweise — Die Machschen 


1,007 


bestimmt, wo 9, der Scharparameter ist, der gerade der Wert von 9 auf der individuellen 
Machschen Welle im Berührungspunkt mit der Einhüllenden r 1,007 ist. 


Grundsätzliche Bemerkungen zur Frage der Grenzlinien folgen in der Zusammenfassung, 


Zusammenfassung. In «dieser Untersuchung wurden exakte Potentialströmungen eines 
adiabatisch zusammendrückbaren Gases für den wichtigen Fall berechnet, daß die Strömungs- 
geschwindigkeit die Schallgesehwindigkeit durehschreitet. Zur Konstruktion ebener und 
stationärer Strömungen dieser Art wurden Kurvenscharen ermittelt, die als Träger konstanter 
Tangential- und Normalgeschwindigkeit möglich sind. Es erwiesen sich Scharen von loga- 
rithmischen Spiralen und von Geraden für diesen Zweck geeignet. Diese Kurvenscharen 
wurden dureh Hinzunahme ihrer orthogonalen Trajektorien zu einem krummlinigen Koor- 
dinatennetz ergänzt. Mit anderen Worten wurden also in einem krummlinigen Koordinatennetz 
ebene Strömungen berechnet, deren Geschwindigkeitskomponenten nur von einer Koordinate 
abhängen. 


Unter Verwendung logarıthmischer Spiralen als Träger konstanter Tangential- und 
Normalgeschwindigkeit wurden verschiedene Strömungen konstruiert. Unter anderen lassen 
sich dadureh Strömungen durch gekrümmte Kanäle mit einem Übergang von Unter- zu 
Überschallgesehwindigkeit darstellen. Auch der von G.J. Taylor und H. Bateman schon 
behandelte Grenzfall, daß die logarıthmischen Spiralen in konzentrische Kreise ausarten, 
wurde einer erneuten Diskussion unterzogen. 

Das auffallendste Ergebnis war bei fast allen Strömungen das Auftreten von Grenzlinien 
im Überschallgebiet, jenseits derer es nieht möglich war, die Strömungen in der gewählten 
Weise, nämlich in alleiniger Abhängigkeit von einer krummlinigen Koordinate fortzusetzen. 
Diese Grenzlinien hatten verschiedene bemerkenswerte Eigenschaften. Der Geschwindiekeits- 
betrag und damit der Druck war auf der Grenzlinie konstant. Das Druckgefälle senkrecht 
zu der Grenzlinie war unendlich, woraus die unendliche Krümmung der Stromlinien an 
demselben Orte folgte. Ferner war die Geschwindigkeitskomponente senkrecht zur Grenzlinie 
eleich der lokalen Schallgeschwindigkeit, so daß die Grenzlinie eine Machsche Welle der 
Strömung darstellte. Darüber hinaus hüllte die Grenzlinie sogar die Machschen Wellen 
ihrer Schar ein. Hierin hat man wohl den Hauptgrund dafür zu erblieken, daß die Strömung 
nicht mehr in der bisherigen Weise fortsetzbar ist. Es ist sogar zweifelhaft, ob die Strömung 
dureh eine Emveloppe von Machschen Wellen hindureh überhaupt noch adiabatisch und 
wirbelfrei weitergeführt werden kann. Die Behandlung dieses neuen Problems bleibe einer 


späteren Untersuchung vorbehalten. 679 
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